Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,218,574 T→C 14.2% G167G (GGA→GGG)  aer ← fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,218,5740TC14.2% 255.4 / 3.3 96G167G (GGA→GGG) aerfused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (64/15);  new base C (6/11);  total (70/26)
Fisher's exact test for biased strand distribution p-value = 3.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCA  >  NC_000913/3218341‑3218808
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
gCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑CTAAACATCAGGGTCATCACTCCAc                                                                                                                                                                                                                                          >  2:438379/1‑238 (MQ=255)
           cAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCCACAATACCCCTAACGCACAGAGGATATAC‑CTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATAAGGGTCACCACTCAACGCC‑CCCGCCacc                                                                                                                                                                                                                             >  1:285912/1‑237 (MQ=255)
           cAGCCCAAGCTCATCGCTGCGATTCAGATGCCCAACACTATTACGTTCCCC‑GGTCCCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGc                                                                                                                                                                                                                             <  2:961058/239‑1 (MQ=255)
            aGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAAACAAAGC‑ATGGCCGCCAGCAAG‑CTAAAAATCAGGGCCATCACTCAACGCG‑CCCGCCAgag                                                                                                                                                                                                                            >  2:1170699/1‑237 (MQ=255)
             gCCCCAGCTCATCGCTGCTATTCAGATGCTCAACACTATTTCGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑AAATTTCCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCCCCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGa                                                                                                                                                                                                                           <  2:968324/239‑1 (MQ=255)
              ccccAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑CTAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGaa                                                                                                                                                                                                                          >  2:1244427/1‑239 (MQ=255)
              ccccAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACCCCACGCG‑CCCGCCAGCGaa                                                                                                                                                                                                                          >  2:1156758/1‑239 (MQ=255)
               ccccccctctccgtcggtttCAGACGCCCACCACT‑TTTCCTTCTCCGGGTCGCCACCTTCAGGCCCGGATGGGCA‑ACATTTTCTATCGGGCACAC‑AATCTGCAATTAAAAACAGGCCTTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCCCCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAg                                                                                                                                                                                                                         <  1:1257515/221‑1 (MQ=255)
               cccattccATCGCTGCTATTCAGATGCTCAACACTTTTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCCCAGAGGATATCC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAg                                                                                                                                                                                                                         <  2:909568/232‑1 (MQ=255)
                ccAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCCCAAAGGATATAC‑GTCACCACCGGGGCAGCAACAACCAAAACC‑ATGGCCCCCACCAAG‑CTAAACACCCGGGGCACCACCCCACCCC‑CCCCCCCACGaacc                                                                                                                                                                                                                        >  1:354476/1‑237 (MQ=255)
                  aGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCCCTAAAGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAGGC‑ATGGCCGCCAGCAAA‑AAAAAAATCAGGTGCATAACCCCAAACG‑CCCCCAACCGaaaggg                                                                                                                                                                                                                      >  1:648624/1‑235 (MQ=255)
                        tCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCACCAACAAACCAAAGC‑ATGGCCGCCCGCAAG‑ATAAACATCGGGGCACTCACCCCAcccc‑cccccccacaaaacggaaagaa                                                                                                                                                                                                                >  1:119952/1‑223 (MQ=255)
                         cGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGAAATGaa                                                                                                                                                                                                               >  2:221754/1‑239 (MQ=255)
                         cGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGCCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGaa                                                                                                                                                                                                               >  2:50121/1‑239 (MQ=255)
                                 ttCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGACATCACCCCACGCG‑CCCCCCACCGCAGCGGAAAACAAGGCAGttt                                                                                                                                                                                                       >  1:1289921/1‑239 (MQ=255)
                                      aTGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCAGGCGCCC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGAACAGAGGATATCC‑GTCACCACCGGGGCCGCAACCACCCAAAGC‑ATGGCCGCCAGCAAG‑ATCAACATACGGGTCATCACTCCACCCC‑ACCGCCACGCAAGCGGCAAAGAAAGCAATTTAcac                                                                                                                                                                                                   >  1:3573/1‑236 (MQ=255)
                                           ccccgcttttcttgtcccc‑GGTCCCCACCTTCAGTCCCTGGTGGGCA‑AAATTTTCCTGCGGGAGCCA‑AATCTCCCATACAAAACAGGCGCTTGCCAACAATTCCACAAACCCACAGAGGATATAC‑GTCCCCCCCGGGGCAGCAACAACCAAAAGC‑CGGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACccagcc                                                                                                                                                                                               <  1:974991/220‑1 (MQ=255)
                                           cAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGACCGCCAGCAAG‑ATAAACATCAGGGCCATCACTCCACCCG‑CCCGCCACCGAACCGGTAATGAAAGCAGTTTAcccaccca                                                                                                                                                                                              >  1:883667/1‑238 (MQ=255)
                                                 tATTACGTTCCCC‑GGTCGCCCCCTTCAGTGCCTGATGGGCA‑ACATTTCCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATCCCACGAACGCACAGAGGATATAC‑GTCACCCCCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTa                                                                                                                                                                                       <  2:979217/239‑1 (MQ=255)
                                                   tttggTTCTCG‑CGTCGCCCCCTTCGTTCCATGATGGCCA‑ACTTTTTCGACCGGCCACAA‑AATCTACCATCCAAAATGGGCCCTGCCCAACAATCCCTATACCGCACGAAGATTACAC‑GTCCCCACCGGGGCAGAACCAAACAAAAGC‑CTGCCCCCCAGCAAG‑ATAAACATCGGGTTCATCCCCCCACGCT‑CCCGCCAGCGAATCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACg                                                                                                                                                                                     <  1:1170699/235‑1 (MQ=255)
                                                   ttACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAATC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAACCGGTAATGACGGCAGTTTACCCAGCCACCCTTTAc                                                                                                                                                                                      >  2:707561/1‑238 (MQ=255)
                                                   ttACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTAc                                                                                                                                                                                      >  2:348408/1‑238 (MQ=255)
                                                        ttCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCACG‑GTAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGAAATGACGGCAGTTTACCCAGCCAGCCTTTACGcacc                                                                                                                                                                                 >  2:1238879/1‑238 (MQ=255)
                                                        ttCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGcacc                                                                                                                                                                                 >  2:11220/1‑238 (MQ=255)
                                                           tcC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCCGCGAACCGGTAATGCAGGCAGTTTACCCAGCCCGCCTTTCCGCCccacca                                                                                                                                                                             >  2:863715/1‑239 (MQ=255)
                                                           tcC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGaa                                                                                                                                                                                                               <  2:859847/205‑1 (MQ=255)
                                                           tcC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGaa                                                                                                                                                                                                               >  1:859847/1‑205 (MQ=255)
                                                           tcC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGcctgtc                                                                                                                                                                               >  2:974991/1‑232 (MQ=255)
                                                           tcC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCAccacca                                                                                                                                                                             >  2:637713/1‑239 (MQ=255)
                                                             c‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAACCATCCGGGTCATCACCCCACGCG‑CCCCCCAGCGACACGGAAATAAAGGCAGTTTACCCACCCAGCCCTTACGCAcccccaa                                                                                                                                                                            >  1:1040574/1‑237 (MQ=255)
                                                                        atttctttcctcgtGGGCACACTTTTTCTATGCGGCCCAA‑AATCCCCCATTCAAAACAGGCGCTTGCCAAAAAACCCACTAACGCAGAGAGGATATAC‑GTCACCACCGGGACAGCAACAAACCAAAGC‑ATGGCCCCCAGCAAG‑ATAAACATCAGGGTCATCACCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGa                                                                                                                                                                  <  1:690536/226‑1 (MQ=255)
                                                                         ttCGGTGCCTGTGGGGCA‑A‑ATTTTTCAACGGGGGCACAAACCTGCAATTCAAAACGGGCCCTGGCAAAAAACCCCACAAACGCCAAGAGGTAAAAG‑GTCACCACCGGGGAACAAACAACCAAAAGC‑AGGGCCGCCAGCAAG‑AAAAACATCAGGTTCATCACCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAAt                                                                                                                                                                <  1:1156758/239‑1 (MQ=255)
                                                                         ttCAGTTCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAACCAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAAt                                                                                                                                                                <  2:1040574/239‑1 (MQ=255)
                                                                                  tGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCCCCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCCAG‑ATAAACATCCGGGTCATCAATCCACGCG‑CCCGCCACGGAACCGGTAATGAAGGAAGTTTACCCCGCAACCCTTTACGCACCACCAGGCCTTATTAAAAACGCtt                                                                                                                                                          >  1:378509/1‑236 (MQ=255)
                                                                                   gATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCttcc                                                                                                                                                        >  2:690536/1‑235 (MQ=255)
                                                                                   gATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTg                                                                                                                                                      >  2:211703/1‑239 (MQ=255)
                                                                                        gCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGAAATGAAGGCAGGTTACCCACCCACCCTTTACGCCCCACCCGGCCTTTAGGAAAACCCtacacggaac                                                                                                                                                  >  1:762225/1‑229 (MQ=255)
                                                                                        gCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCAAGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTAc                                                                                                                                                  >  2:841456/1‑238 (MQ=255)
                                                                                               ttttCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTCCGCCCCCCCCGGCCTTTATGAAAACGCTTCCGGGAACGGCCggg                                                                                                                                           >  1:738867/1‑238 (MQ=255)
                                                                                                    taccGGGCCCAC‑AATCTGCCATTCAAAACGGGCCGTTGCCAAAAATACCACTAACGCACAGAGGATATCC‑GTCCCCACCGGGGCAGCAACAAACCAAAGC‑AGGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCa                                                                                                                                      <  2:802813/236‑1 (MQ=255)
                                                                                                    caccGGGCACAC‑AATCTGCAATTCAAAACAGGCGCTTGCCAACCACACCACAAACGCACAGGGGATTGAC‑GTCACCACCGGGGCAGCAACAAACCAAGTC‑ATGGCCGCAAGCAAG‑ATAACCATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCa                                                                                                                                      <  2:707137/236‑1 (MQ=255)
                                                                                                       cGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTCCTGGAACGCCCGGCGTTCAAc                                                                                                                                    >  1:838175/1‑238 (MQ=255)
                                                                                                         ggCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCCCCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGACATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGAACGTCCGGccgtcacac                                                                                                                                   >  2:487793/1‑229 (MQ=255)
                                                                                                                           tCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCCCTCCACGCG‑CCCGCCAGCGAAGCGGTAATTAAGCACATTTACCCAACAAACCTTTCAGCACCCCCAGGCCCTTATTAAACCGCCTAACGGGAAGTCCGGCGTTCAAACCTTCGTAA‑agcagccaca                                                                                                                >  1:307294/1‑234 (MQ=255)
                                                                                                                           tCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCCCCAGGCCTTTATGAATACGCTTAATGGTACGTCCGGCGTTCAACGCTTTGGAA‑AGCGGCCCCa                                                                                                                >  1:1126630/1‑238 (MQ=255)
                                                                                                                             aaaaCAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTCCGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCAcc                                                                                                              >  2:192349/1‑238 (MQ=255)
                                                                                                                                  cggCCCTTGCCAACAAAACACGCAACGCCGAGGGGATAACC‑GCCACCACCGGGCCAGCAACAAACAAAGGC‑AGGCCCGCCAGCAAG‑AAAAACAGCAGGGTCACCCCCCCCCGCC‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCcgcg                                                                                                        <  1:211703/238‑1 (MQ=255)
                                                                                                                                  ccctcctggcaaaaacacccACGAACCCAC‑GAGGGTATACGGCCCCCCCCGGGGAAACAACAACCAAAACC‑ATGGCC‑CCAGGAAGCATAACCTGCGGGTTGACCTCCCCCCGCC‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTAGGAATTCGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCcgcg                                                                                                        <  1:637713/221‑1 (MQ=255)
                                                                                                                                    gcgcTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCCGCGACGCGGTAATGAAGGACGTTTACCCATCCACCCTTTCCGCACCACCAGGCCTTTATGACTACGCTTCCTGCTACCTCCGACGTTCAACGCTTCGTAC‑AGCAGCCCCACCGCCCCGa                                                                                                       >  1:798532/1‑238 (MQ=255)
                                                                                                                                         tGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGCCCGCCAGCAAG‑ATAAACATCAGGGTAATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAAGCAGTTTACCCAGCCAGCCTTTACGCCCCACCAGGCCTTTATGAAAACCTTAACTGGTACGCCCGGCGTTAAAAGCCTTTTAC‑AACGTCTCCACCCCCCCGGTCCCtt                                                                                                 >  1:1182148/1‑239 (MQ=255)
                                                                                                                                                      aCTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGAACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGCTCAACCCTTTGTCC‑ACCGGCCCCCCACCCGGGACCTCTtcctcccccccccg                                                                                    >  1:1080042/1‑231 (MQ=255)
                                                                                                                                                        tAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGGCCGGCGTTCAACGCTTTGTAC‑AGAGGCACCCCCCCCGCGATCTCTTATACCGCCGCCCggg                                                                                  >  1:1275991/1‑239 (MQ=255)
                                                                                                                                                          aCGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTAt                                                                                                                                                                    >  2:690291/1‑156 (MQ=255)
                                                                                                                                                          aCGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGCATACGCTTACTGGTACCTCCGGCGTTCAACGCTTTGTAC‑ACAGGCTCCACCGCCGCGACTCCTTCATCCCTCCCCCggtt                                                                                 >  1:169827/1‑236 (MQ=255)
                                                                                                                                                          aCGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATAc                                                                                                                                                              >  2:253478/1‑162 (MQ=255)
                                                                                                                                                          aCGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGt                                                                                 >  2:115004/1‑238 (MQ=255)
                                                                                                                                                          aCGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCCTTACGCACCACCAGGCCTTTCTGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTAC‑AACGGCTCCACCGCCCCGCACTCTTCATCCGTCGCCCgggg                                                                                 >  1:554400/1‑237 (MQ=255)
                                                                                                                                                               caGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCaccacc                                                                                                                                                                              >  1:567923/1‑141 (MQ=255)
                                                                                                                                                               caGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCaccacc                                                                                                                                                                              <  2:567923/141‑1 (MQ=255)
                                                                                                                                                               acagggaataacctcccccccggggcCAACAACAACCCAAAGC‑AGGCCCCACAGCAAG‑AAAAACATCGGGGTCATCCCTCCCCGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGaa                                                                            <  2:354476/217‑1 (MQ=255)
                                                                                                                                                                    gATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCAACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTCCGCACCACCAGGCCTTTATGAATACGCTTACTGGTCCGGCCGGCGGCCACCCCTTT‑TACAAACGGTCCCACCCCCCCCATCCCTTCTACCGGCCCCCGGGTACGaaacacca                                                                      >  1:1226716/1‑233 (MQ=255)
                                                                                                                                                                       ataC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGc                                                                                                                                                                          >  2:465690/1‑137 (MQ=255)
                                                                                                                                                                             tCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATg                                                                                                                                                                   >  2:487592/1‑139 (MQ=255)
                                                                                                                                                                             tCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGACCTCTTCATCCGTCGCCCGGGTACGAATCGCCTAATAGcc                                                               >  2:557414/1‑238 (MQ=255)
                                                                                                                                                                             ccaccacGGGGGCAGCAACAACCCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCCCTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCa                                                              <  2:1289921/238‑1 (MQ=255)
                                                                                                                                                                               accaccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTCCGTCCGGCGGTCAACGCTTTGTAC‑AACGGCCCCCCCCCCGCGATCCCGTCATCCCTCCCCCGGGTAcacacccacaaaaaaccccc                                                            >  1:438316/1‑220 (MQ=255)
                                                                                                                                                                               accaccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGaa                                                                            >  2:197318/1‑223 (MQ=255)
                                                                                                                                                                               accaccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGCACGAATCGACATATAGCCACt                                                            >  2:487797/1‑239 (MQ=255)
                                                                                                                                                                               accaccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCCCTTCATCCGTCGCCCGGGTACGCATCGACATATAGcc                                                               >  2:119170/1‑236 (MQ=255)
                                                                                                                                                                               accaccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGCCGTTCAACGCTTTGTCA‑AACGGCCCCACCGCCGCGGCCCCCTCATCCCACCCCCGGGGACGaaacagacaaaaacc                                                               >  1:960596/1‑223 (MQ=255)
                                                                                                                                                                                 caccGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTTTAC‑AGCGGCCCCACCGCCGCGCCCTCCTCACCCGTCCCCCCGGTACGAAACGACataaaacccccga                                                          >  1:388034/1‑229 (MQ=255)
                                                                                                                                                                                       ggaagcaACAACCCAAACC‑ATGGCCGCCAGCAAA‑ATAAACATCAGGGTCATCCTCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGAttttt                                                     <  2:169827/235‑1 (MQ=255)
                                                                                                                                                                                        gcagcaACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCAACCGCCGCGATCTCTTCATCCGTCGCCCTGGTACGAATCGACATATAGCCACCGATTTTTc                                                    >  2:37983/1‑238 (MQ=255)
                                                                                                                                                                                           gcaACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCCTCCGTCGCCCGGGTACGATTCGACATATAACCACTGCTTTTTCCCt                                                 >  2:144303/1‑238 (MQ=255)
                                                                                                                                                                                           gcaACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCAt                                                                                              >  2:849896/1‑193 (MQ=255)
                                                                                                                                                                                                aaccccaagcccctccccacAGCAAT‑ATAAACATCGGGGT‑TTCCCCCCCCGCGCCCCGCCGGCGAAGCGGTAATGAAGGAATTTTCCCCAGCCAGCCTTGACGCACCCCCCGTCCTTTATGAATGCGCTTACTGGTACGTCCGGCGTTCACCGCTTTGTAC‑AGAGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCgcgc                                            <  1:557414/220‑1 (MQ=255)
                                                                                                                                                                                                aaaaCAAAGCAAGGGCCGCAGACAAG‑TAAAACATCGGGGTCATCCCCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAGGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATCCGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCgcgc                                            <  2:1240116/239‑1 (MQ=255)
                                                                                                                                                                                                  cccAAAGC‑AGGGCCGCCAGCAAG‑AAAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACg                                                                              <  2:1193567/201‑1 (MQ=255)
                                                                                                                                                                                                      aaGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCg                                     >  2:201217/1‑239 (MQ=255)
                                                                                                                                                                                                             ggCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCTTCCGCCGCCCGGGGACGAACCgcccat                                                                    >  1:1193567/1‑197 (MQ=255)
                                                                                                                                                                                                              gccgccAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATGCACAAATAGCCCCTGTTTTTTCCCTTCCGCACCATCGGTACCGc                              >  2:578808/1‑239 (MQ=255)
                                                                                                                                                                                                                  ccAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGATTCGACATATAGCCACCACTTTTTCCCTCGCGCACCATTGGTACCCCATTg                          >  2:986211/1‑239 (MQ=255)
                                                                                                                                                                                                                      aaaG‑ATAAACAGCAGGGCCACCCCTCCCCGCC‑CCCGCAAGCGAAGCGGTAATGAAGGCAGTTTCCCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGccc                      <  2:1126630/238‑1 (MQ=255)
                                                                                                                                                                                                                       aaaaaTAAACATCAGGGTCATCACCCCAGCCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGGTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGccc                      <  2:648624/235‑1 (MQ=255)
                                                                                                                                                                                                                        aG‑AAAAACATCAGGGTCACCCCCCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTCCGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGc                    <  2:119952/239‑1 (MQ=255)
                                                                                                                                                                                                                                     gggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGTTCCCTTCATCCGTCGCCCGGGTACGAATCGCCATTAAGCCCCTGATTTTTCCCCCGCGCCCCATCGTGACCCCATTTGCCCCCACCCCATAATg         >  1:703889/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAGTCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCCTTGGCCCCCACCCcata            >  2:1244539/1‑235 (MQ=255)
                                                                                                                                                                                                                                     gggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATgg        >  2:702198/1‑239 (MQ=255)
                                                                                                                                                                                                                                      ggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACg                                                                                                                                   <  2:786386/116‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACg                                                                                                                                   >  1:786386/1‑116 (MQ=255)
                                                                                                                                                                                                                                      ggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTCC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTTGCGCCCCATCCGTACCCTATTGGCCCGAACCCCATCATgg        >  2:720484/1‑238 (MQ=255)
                                                                                                                                                                                                                                      ggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑ATCGGCTCCACCGCCGCGATCCCTTCCTCCGTCGCCCGGGTCCGAACCGGCATATACCCCCTGATTTTTCCCTCGCCCCCCCACGGTACCCCATTGGCCC‑CACCCaaaaaaagg       >  1:258518/1‑234 (MQ=255)
                                                                                                                                                                                                                                      ggTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCCCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATgt        >  2:130751/1‑237 (MQ=255)
                                                                                                                                                                                                                                          atcaCTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCAcc   >  2:175458/1‑239 (MQ=255)
                                                                                                                                                                                                                                          atcaCTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGCCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTCCGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCCCATTGGCCCGCACCCACTAATGGTcc     >  2:733350/1‑236 (MQ=255)
                                                                                                                                                                                                                                           tcaCTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCCCTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCa  >  2:1165664/1‑239 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
GCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCC‑GGTCGCCACCTTCAGTGCCTGATGGGCA‑ACATTTTCTATCGGGCGCAC‑AATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACAGAGGATATAC‑GTCACCACCGGGGCAGCAACAAACCAAAGC‑ATGGCCGCCAGCAAG‑ATAAACATCAGGGTCATCACTCCACGCG‑CCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTAC‑AGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCA  >  NC_000913/3218341‑3218808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑