Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C93.1% 451.9 / 24.1 164intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base C (66/83);  total (73/91)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
*NC_0009134,296,3802.G93.2% 463.4 / 23.0 164intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base G (69/83);  total (73/91)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 164 total aligned reads displayed.

GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTT  >  NC_000913/4296193‑4296617
                                                                                                                                                                                              ||                                                                                                                                                                                                                                                   
gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGAAGCTTGCCGCCTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                                   >  1:1172538/1‑251 (MQ=2)
     ttctgaggaccaCGCCAGGCAGCGGAATAGCCTGAGTTAGCTGGAGTTTGTGGGTCGGTTAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGGCAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                              <  1:1910638/240‑1 (MQ=11)
     ttatacggCCTACCCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                              <  1:886023/245‑1 (MQ=255)
     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                              <  2:420435/251‑1 (MQ=17)
     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                              <  2:148164/251‑1 (MQ=17)
     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                              <  2:1255538/251‑1 (MQ=17)
         cAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATTTTACCCGTTGGTTac                                                                                                                                                                          >  2:1083879/1‑250 (MQ=9)
               tacaccacccAGCACAATAGCCTGATTTCCCGTTATTTTGGATGTCGGATAAGGCGTTTATACCGCATCCGACGTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCGTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                                   <  2:983339/242‑1 (MQ=11)
                accccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGTCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                  <  2:679072/248‑1 (MQ=255)
                accccACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                  <  2:1031557/248‑1 (MQ=255)
                aCGCCAGACAGCGCAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                  <  1:1362633/251‑1 (MQ=255)
                aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                                                                                                                             <  2:163391/241‑1 (MQ=17)
                aCGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                  <  2:1321344/251‑1 (MQ=255)
                  gCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                                                                                                <  2:1172550/251‑1 (MQ=255)
                    cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTTGGTCGGACAAGGCGTTTA‑‑CGTCGCATCCGACATCAAGGCCCGATGCGACGCTTGCCGCGTTTTACCAGGCCTATCTTAACCGGTGGTTAATTTTTCaaaac                                                                                                                                                              >  1:679065/1‑249 (MQ=11)
                         aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATAAATGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTATCTTACCGCTTCGTTAATTTTTCAACCTCTTca                                                                                                                                                         >  1:1280514/1‑251 (MQ=11)
                          gcgcgtaaGCCTGATTTAACGTAATTTTGTAGTTCGGATAAGGCGGTTAAGCCGCATCCGTCATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  1:1200149/244‑1 (MQ=14)
                          gcgcTATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  1:1805022/251‑1 (MQ=255)
                          gcgcGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  2:279737/251‑1 (MQ=255)
                          gcgcAATAGCCTGTTTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCTCTTCCCACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  2:950682/251‑1 (MQ=14)
                          gcgcAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAGGCGACGCTTAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  2:613898/251‑1 (MQ=255)
                          gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAGGCTTAATGCGTCTTATCAGGCCGAAGTCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTCTAGGTCGGAGCAGGCGTTTACGC‑‑CGCATACGACCTCAATCCCGGCTGCTACGCTTCCCGAATCTTACAAGGCCTATCTTACCCGTGGGTCAATTTTTCAAACACTTcac                                                                                                                                                        >  1:139063/1‑251 (MQ=11)
                          gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  2:1874792/251‑1 (MQ=255)
                          gcccaatagctgtatgtggCGTGATTTTGTAGGTCGGATAAGGCATTTACGCCGCATCCGGCATCAACGCCTGGTGCGACGCTTACCGCGTCTTATCAGGCCTACGCCATACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  1:652219/233‑1 (MQ=14)
                              aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGC‑‑CGCATACGACACCAATGCCTGATGCGACGCTTGCCTCGCCTTACCAGCCCTATCTTAACCGTGGGTTAATTTTTCAAACTCTTCTCAccc                                                                                                                                                    >  1:97283/1‑251 (MQ=14)
                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGCATTAGCGTGATTTTGTCGGTCGGATCAGGCGTTTACGC‑‑CGCATTCGACCTCAATCACTGATTCGACGCTTTCCGCGTCTTAACAAGCCTCTCTTACCCGTTGGTAAATTTTTCACACTCTTCACACGCGGTAt                                                                                                                                               >  1:100420/1‑251 (MQ=14)
                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                                                                                                                               >  2:1129744/1‑251 (MQ=255)
                                    cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                              <  2:866626/251‑1 (MQ=255)
                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCGGCGCAATAGCCTGCATTAGCGTGATTTTGTAGGTTGGATAATGCGTTAA‑‑CGCCCAAGCCGACATCCAGGCCTGACGCGCGCTTTGCCCACTCTTAACAGGCTTATCTTAACCGTCGGTTAATTTTCCAAACTCTTTACAACCGCTATCCAAc                                                                                                                                          >  1:1775284/1‑251 (MQ=14)
                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGTCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTAAATTTTTCACACTGTTCACACCTGGTATCAAAc                                                                                                                                          >  2:861147/1‑251 (MQ=14)
                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATCCCTGATGCGACGCTTGCCCTGTCTTATCAGCCCTATCTTCATCGTTGGTTAATTTTTCACACCCTTCACACCCGGTATCAAAc                                                                                                                                          >  1:574646/1‑251 (MQ=14)
                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACCCTTCACACCCGGTATCAAAc                                                                                                                                          >  2:15332/1‑251 (MQ=255)
                                         ttggcgggttttcGTATGTCGGATAACGCGTTTTTGCCGCTTCCGACTCCACCGCTTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                         <  1:377147/242‑1 (MQ=17)
                                              tcagatgtgtataaggaacagggCGTTTGTGCCGCATCCGACATCAACGGCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                    <  1:1208432/230‑1 (MQ=17)
                                              gTGATTTTGTATGTCGGATAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                    <  2:1730400/251‑1 (MQ=255)
                                              gTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                    <  2:507844/251‑1 (MQ=255)
                                                  ttttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATGCCAGAAAGCGCAATAGCCTGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCCGATGCGACGCTTGGCGCGGCGTATCAGGTCTATCTTACCAGTTGGTTAATTTTTCAAAATCTTCCAAACCGGTTGCAAAcccatacaaac                                                                                                                                >  2:1253605/1‑243 (MQ=11)
                                                   tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                                                                                               <  2:527720/251‑1 (MQ=255)
                                                    ttGTAGGTCGGAAAAGGCGTTTAGGCCCCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAGACTCTTCACACCCGGTATCAAACCCTTCCATACAg                                                                                                                              <  1:310986/251‑1 (MQ=255)
                                                     tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                                     >  2:1829711/1‑227 (MQ=255)
                                                         ggTCGGATAAGGCGTTTATGCCGCAACCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGTGATTTTGTACGCCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATCCCTGAGGAGAGCCTTGCCGCCTCTTCTCAGGCCTATCTTAACCGTGGGTTAATTTTCCAAACTCTTCCCCCCCGGTAACAAATCTTTCCATACACCTccg                                                                                                                         >  1:1028301/1‑249 (MQ=14)
                                                          gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CTCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCACACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                                                                                                                        >  1:1219159/1‑251 (MQ=255)
                                                          gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTCGTAGGTCGGATAAGTCGCTTACGC‑‑CGCATCAGACATCAATGCCTGATGCGAAGCTTGCCGTGTCTTACCAGGCCTATCTTAACCGTTGGATAATTTTGTCACATCTTCACAATCGGCATCAAACCCTTACATACCGCTCCgg                                                                                                                        >  1:1702861/1‑251 (MQ=14)
                                                              gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTGTCAAACTCTTCACACCCGGTATAAAACCCTTCCATCCAGCTCAGGtttc                                                                                                                    >  2:86997/1‑249 (MQ=255)
                                                                   ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                        <  1:1620630/210‑1 (MQ=255)
                                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc                                                                                                         >  1:1448685/1‑251 (MQ=255)
                                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCCGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc                                                                                                         >  2:488095/1‑251 (MQ=255)
                                                                            gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc                                                                                                      >  2:1163981/1‑251 (MQ=255)
                                                                             ccacATCCGCCATCACCGCCTAATGCTACGCTTAACGCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCGTACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                                     <  2:1002169/248‑1 (MQ=18)
                                                                             ccGCCTCCGACATCAACGCCTGATGGGACGCTTAACGCGTTTTCTCAGGCCTATGCCAGACAGCGGAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                                     <  1:1074615/251‑1 (MQ=18)
                                                                                 aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTg                                                                                                 >  2:245290/1‑251 (MQ=255)
                                                                                      ccagcaacccctgagcctccGCTTAGCGCCGCTTACCAGACCGACGCATGAAAGCGAAAGTGCCTAATTTAGCGTTATTTTGTAGGTACGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATTCCTGAGTCGACGCTGTCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAGCTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTGCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                            <  1:850931/232‑1 (MQ=16)
                                                                                      acagtcaggCCTAATGCGACGCTTAACGCGTCTTATCAGACCGACGCCAGACTGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                            <  2:53997/243‑1 (MQ=21)
                                                                                      aCATGAACGCCCGATGCAAAGCTTAACGCGTCTTTTCAGGGTTCCGCACGACGACTCAGTAGCCTGATTTAGCGTGATTTGTTAGGTCGGATAAGGCGTTTACGC‑‑CGCAACCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                            <  1:1864275/251‑1 (MQ=21)
                                                                                      aCATCAACGCCTGGTGGCACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                            <  2:1280529/251‑1 (MQ=255)
                                                                                           aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                >  2:1620667/1‑210 (MQ=255)
                                                                                           aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                                                                            >  2:1208446/1‑230 (MQ=255)
                                                                                              gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGC‑‑CGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                    <  1:152316/251‑1 (MQ=255)
                                                                                              gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                    <  1:245289/251‑1 (MQ=255)
                                                                                              gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                    <  1:1621852/251‑1 (MQ=255)
                                                                                               ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGA‑‑CGCATCCGAAATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCTGGTATCAAACCCTTCCATACAGCTCAGGTTCAAACAGTGCAGAGCATTTAGCATGttcgctacga                                                                                   >  2:383871/1‑244 (MQ=18)
                                                                                                         cccttcacgctgattaccagtcatccGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCTTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCTGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                         <  1:370172/226‑1 (MQ=25)
                                                                                                          gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                                     >  2:241154/1‑222 (MQ=255)
                                                                                                          gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCCCCTTTTGCTTGTTCTCTTCGATAAAACcttct                                                                        >  1:1005495/1‑251 (MQ=255)
                                                                                                                cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCg                                                                  >  2:1065093/1‑251 (MQ=255)
                                                                                                                   gtgttACCAGGCCTACGCCCGACAGTGTAGTCGCCTGAGTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTA‑‑CGTCGGCCACCGCAGTAATGGCTGATGAGACAGTGGTCGCGTCTTATCAGGCCTATCATATCCGTTGGTTAATTTTTCACACTCGTCACGCCCGGTATCAAAACCTTCCATACAGCTAAAGTTCGACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTACACCGttt                                                               <  1:1515235/248‑1 (MQ=255)
                                                                                                                     cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACTGTTGGTTAATTTTTCAACCTCTTCACAcc                                                                                                                                                     >  2:846774/1‑163 (MQ=255)
                                                                                                                        agcAGGCGTAGGCGAGACAGCGCAATAGCCTGATTTAGCTGGATTTTGTAGGTCGGGTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                          <  2:1775325/249‑1 (MQ=37)
                                                                                                                        aTCAGGCTGACTCCAGACAGCGAAAGAGCCGTATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTGTCG‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCACAAGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                          <  1:628561/251‑1 (MQ=255)
                                                                                                                        aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                          <  2:661831/251‑1 (MQ=255)
                                                                                                                          cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAca                                                        <  2:1517922/251‑1 (MQ=255)
                                                                                                                              ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                    <  2:1615417/251‑1 (MQ=255)
                                                                                                                               ttACGCCAGGAAGCGCAATAGCATGATTTAGCGTGAGTTTTTAGGTAGGAGAATGCGTTTACGC‑‑CGGAACCGGAATCAATGCCTGGTGCGGCGATTGCCGCGTCTTATCAGGCTTATCTTAACGTTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTTAACCAGGGCAGCACCTTTTGCGTGGTCTTCTCGATAAAACCTTCTGCACCGTTTAAGAAcatcatc                                                   <  1:1437027/250‑2 (MQ=255)
                                                                                                                               cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                                   <  1:240563/251‑1 (MQ=255)
                                                                                                                                   gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTTGCAGCGTATTATCAGGCCCATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                >  2:1369081/1‑170 (MQ=18)
                                                                                                                                   gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCCTCTTCACCGTTTAAGAACATCATTCCCg                                               >  1:1615380/1‑251 (MQ=255)
                                                                                                                                       gACATCGCAGTAACCGTATTTAGCGTGGTTTTGTAGCTCGGATAACGCGTTTACGC‑‑CGCATTCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCCACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCa                                           <  1:488091/251‑1 (MQ=37)
                                                                                                                                        aCAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTTATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTACTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAg                                          <  1:1923165/251‑1 (MQ=255)
                                                                                                                                        aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAg                                          <  2:557393/251‑1 (MQ=255)
                                                                                                                                          aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTa                                        >  2:1018799/1‑251 (MQ=255)
                                                                                                                                              cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                    <  2:185683/219‑1 (MQ=255)
                                                                                                                                                  aGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt                                <  1:601487/251‑1 (MQ=255)
                                                                                                                                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTCCTTCACAGTTTAAGAACATCATTCCCGCCTAGTACTCGGAATaa                              >  1:1245595/1‑251 (MQ=255)
                                                                                                                                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCCTCCATACAGCCCAGGTTCAACCAGTGCAGCGCCTTTTGCTTTTTCTTCTCGATAAAACCTTCTTCACCGTTCAAGACCATAATTCACGCCCACTACTTGCAATac                              >  1:1321734/1‑250 (MQ=25)
                                                                                                                                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTGCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATaa                              <  2:453682/251‑1 (MQ=255)
                                                                                                                                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTc                                     <  1:1058728/241‑1 (MQ=255)
                                                                                                                                                         tttAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTCCACCGTTTAAGAACATCATTCCCGCCCAGTCCTCGGAATAACCGGt                         >  2:1864061/1‑251 (MQ=255)
                                                                                                                                                         tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGTCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTATTTTTTCAAACTCTTCATACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACTTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGt                         >  2:1034239/1‑251 (MQ=255)
                                                                                                                                                         tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAAATCATTCCCGCCCAGTACTCGGAATAACCgtt                         >  1:1832476/1‑249 (MQ=255)
                                                                                                                                                         tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGt                         >  2:1908273/1‑251 (MQ=255)
                                                                                                                                                         tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTCACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTAAAACTATTCACACCCGGTATCAAACCATTCAATACAGCCCAGGTTCAACCAGTGCAGCGCCTTTTCCTTGTTCTTATCGATCAAACCTTCTTCACCGTTTGAGAACATCCTTCCCGGCCAGTACTCGGAATacacgct                         >  2:951892/1‑245 (MQ=25)
                                                                                                                                                         tttAGCGTGATCTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGTGTCTTATCAGGCCTATCTTAATCGTTGGTTAATTTTTCAAACTATTAACACCCGGTATCAAAACCTTCCATACAGCCCAGGTTCAACCAGTGCAACGGCATTTGCTTGATCTTTTTGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                    >  1:1562012/1‑224 (MQ=25)
                                                                                                                                                            aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg                      <  1:1319240/251‑1 (MQ=255)
                                                                                                                                                            aGCGTGATTTTGTAGGTCAGATAAGCCGTTCA‑‑CACAGCACCTGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCTTTGGTTAATTTTTCATACTCTTCACACCCGTTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCTCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg                      <  2:1289278/251‑1 (MQ=38)
                                                                                                                                                                 gATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCCTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTCCCATACAGATCAGGTTCCATCAGTGCAGCGCCTTTCGCTTGTTCTTCCCGATCAAACCTTCTTCACCGTTTACGAACACCATTCCCGCCCATTACTCGGCATAACCGGTTGTGGaaa                 >  1:1838601/1‑251 (MQ=37)
                                                                                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAAtt               >  2:12988/1‑251 (MQ=255)
                                                                                                                                                                         atggctgATAAGTCGTTTACGC‑‑CGCATCCGACATCAATGCCTTATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTACTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCagag         <  1:954983/245‑1 (MQ=255)
                                                                                                                                                                               gATAAGGCGTTTA‑‑CGGCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTACAGAGCTGCGt   <  1:1320265/251‑1 (MQ=255)
                                                                                                                                                                                aTAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGtt  <  2:1005502/251‑1 (MQ=255)
                                                                                                                                                                                 tAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                    <  2:1414074/200‑1 (MQ=255)
                                                                                                                                                                                              ||                                                                                                                                                                                                                                                   
GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTT  >  NC_000913/4296193‑4296617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: