Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,438 T→G 9.9% pseudogene (1047/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,4380TG 9.9% 173.7 / 2.7 66pseudogene (1047/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base G (3/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAG  >  NC_000913/1466141‑1466738
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                  
aGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                                                                                                                                                                                                                                                                                 <  2:617642/301‑1 (MQ=255)
 gTGATTCAGGTCGGCGATCTGGAGTTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGTGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCAGCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                                                                                                                                                                                                <  2:611271/301‑1 (MQ=255)
  tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGCTAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTAATTCa                                                                                                                                                                                                                                                                                                               >  2:58396/1‑301 (MQ=255)
  tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCa                                                                                                                                                                                                                                                                                                               <  2:579487/301‑1 (MQ=255)
                           cAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAGCACGGCCGTGCGCCG‑TGCGGGGACAGTGTGTATCCAGATGTACGTCGCCATGGCCGCCGTCTAAAACGAAGACGAGAGCACCATCCCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACAGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGATATTCAGAACGACGAACTGGATGTCAgcgag                                                                                                                                                                                                                                                                                      >  1:104439/1‑299 (MQ=17)
                                gcgggggTCACGGTATTGATATCACCGGTGACAGTGCGACGGGGGATAACAAGGGCACCATGACCGT‑CACCGATTCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                                                                                                                                                                                                 <  2:473978/301‑1 (MQ=255)
                                    cggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CAc                                                                                                                                                                                                                                                                            >  2:143226/1‑301 (MQ=11)
                                          gggtattgtttgtcccGGTGAAAGCGTGGCGTTGGATAACAAGGGCACCATGCCCGT‑GACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATGTCACCGG‑TACGGAAAGTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGTGGCGGT‑CACGGTAtt                                                                                                                                                                                                                                                                      <  1:143100/287‑1 (MQ=11)
                                           ggtattgttattatcgtgtactgcggatcGGTGGTGAACAGGGGCACCATGATCGT‑CATCTAGTCGGAGTCCATGGGTATGCAGATTGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGGTCACCATCACCAACGGTGGCACCGGCACCCAGATTATGGGTGATGACGCTACGGCATACATCATCGGCAGGCGCACCGTGGACGGCACAGCTCCCACCGC‑TCCGTCCCTACACGGCGATATTCCGAACGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTg                                                                                                                                                                                                                                                                     <  1:58328/273‑1 (MQ=9)
                                              aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGata                                                                                                                                                                                                                                                                  >  2:639416/1‑301 (MQ=255)
                                                      ctcccctacacggctgtgttccggtgttaaGGGACCCGTGACCGT‑CACCTGTCCGGAGTCAAGGGGGTTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGc                                                                                                                                                                                                                                                          <  1:12072/273‑1 (MQ=17)
                                                             gACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCACCGGTAATAACGGAAAGGTGATTCAGGACGACGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGaacggg                                                                                                                                                                                                                                                   >  1:408611/1‑296 (MQ=255)
                                                               cAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCa                                                                                                                                                                                                                                                 >  2:500692/1‑301 (MQ=255)
                                                                       ggttttatgaCAAGAGTAGCATGACCAT‑CTCCGCTCCTGCGTTCTTGGGTATTCAGATCGATGGTGACAAGGCCGTCGTCAAAAACGAAGGCGAGAGCACCAACACCAACGGCGACACCGGGATCCATATTACGGGTGATAACGCCACGGTAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGa                                                                                                                                                                                                                                         <  2:194732/292‑1 (MQ=17)
                                                                         ggggAGAAGAAGGGCACCATGACCGT‑GACCGAGCCGGAGTCCATGGGTATCCAGATCGCCGGTGCCAAGGCCATCGTCAATAACGAAGGCGCGAGCACCTTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATa                                                                                                                                                                                                                                       <  1:500045/299‑1 (MQ=255)
                                                                          tGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATaa                                                                                                                                                                                                                                      >  2:137016/1‑301 (MQ=255)
                                                                             aTAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                                                                                                                   >  2:290565/1‑301 (MQ=255)
                                                                                aCAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                                                                                                                                <  2:641139/301‑1 (MQ=255)
                                                                                 cAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGcagatcggaagagcacacgtctgaactc                                                                                                                                                                                                                               >  2:12072/1‑274 (MQ=255)
                                                                                 aaaGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGc                                                                                                                                                                                                                               <  2:104563/300‑1 (MQ=255)
                                                                                     ttggagttcagacgtgtgctcttccgatctcATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                                                                                                                            <  2:211075/271‑1 (MQ=255)
                                                                                     ggCACCATGACCGT‑CACCGAGTCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATAACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                                                                                                                                                                                           <  1:638539/301‑1 (MQ=255)
                                                                                               ccGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑Ac                                                                                                                                                                                                             <  1:290163/301‑1 (MQ=255)
                                                                                                cGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATCTCACCGGCGACAGCGCAACGGTTGATAACAAGGGCACCATGACCTTC‑‑‑‑Acc                                                                                                                                                                                                            >  2:565640/1‑301 (MQ=255)
                                                                                                        gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACAGGCGACAGCGCAACGGTGAATAACAAGGGAACCATGACCATC‑‑‑‑ACCGATACgg                                                                                                                                                                                                     >  2:188029/1‑301 (MQ=255)
                                                                                                                   cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCagatcggaagagcgtagtgtagggaaagag                                                                                                                                                                                              >  1:210832/1‑272 (MQ=255)
                                                                                                                           tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGt                                                                                                                                                                                  >  2:426710/1‑301 (MQ=255)
                                                                                                                              tcagacgtgtgctcttccgatctgTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGa                                                                                                                                                                               <  2:11302/278‑1 (MQ=255)
                                                                                                                                gttgACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATTGCTGGCCATAACGGGAAGGTGATTCAGGACGGCGACCTGGATGTCAGCGACGGCGGT‑CACGGTATTGATATCACCGGCCACAGCGCACCGGCGGAACACAAGGGACCCATGCCACAC‑‑‑‑ACCACCCCCGAGTCCACCGGGCTCCAAGTTCAcc                                                                                                                                                                             >  1:333885/4‑300 (MQ=17)
                                                                                                                                  cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCCACACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGCACTCAACGGTACTAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGc                                                                                                                                                                           >  2:61943/1‑301 (MQ=255)
                                                                                                                                              cgCCTTCGTTATTGACGATGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCAACGGTATCCAGATTGACGGAGACCAGGCGgtg                                                                                                                                                               >  2:501089/2‑300 (MQ=255)
                                                                                                                                               gCCATCGTCATTAACGACGGCGCGAGCTGCATCTCCAACCGTGGCACCGGCACTCTTATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttg                                                                                                                                                              <  2:433677/301‑1 (MQ=14)
                                                                                                                                               gCCATCGTAAATAACGAAGGCGATAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTTCTGACGCAACGCCACACAACAACGGCAAACCCACCGGTGACGGCAAGGATTCCTCCGG‑TACTGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttg                                                                                                                                                              <  2:129666/301‑1 (MQ=14)
                                                                                                                                                     gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGACAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAATCCAGCGGGACCCAGGtagcagaacggaaaagccgcctgtagg                                                                                                                                                        >  1:11302/1‑275 (MQ=14)
                                                                                                                                                      cggtctaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGGTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                                                                                                       <  1:29042/295‑1 (MQ=255)
                                                                                                                                                       cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAATTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                                                      <  2:250214/301‑1 (MQ=255)
                                                                                                                                                       cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAAGCGGTTGTTaacaaa                                                                                                                                                      >  2:138879/1‑300 (MQ=255)
                                                                                                                                                            aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAGCGGCAAAACCACCGTTGGCGTCAAGAATGCCACCGA‑TACCGAAATCAACGGTAACAACGGAAAAGTCATTCACGCCGGCGATCTGGATGTCAGCGCCCGCGGT‑CACGGTATAGACATCACCGGCGACACCCCCACGGTGGATAACACGCGCACCCTGACCGTC‑‑‑‑ACCCAACCCGACTCCACCGGCACCCAAGTCGCCGGccaccaccacctcgcccaccacccccc                                                                                                                                                 >  1:384307/1‑280 (MQ=255)
                                                                                                                                                            aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAAAGatcg                                                                                                                                                 >  2:29093/1‑298 (MQ=255)
                                                                                                                                                                   ggagaGCATAATCACCAACTGTGGCATCGTGACCGAGATTAATGGGGATGACGCTACGGCAAACAACAACGGAAAAACAACCGTGGACGGCAAGGATTCCAACGG‑GACGGAAATCAACGGTAATAACGGAAAAGGTATTCAGGACGGCGATCTGGATGTAAGCGGCGGCGGT‑CACGGGATTGATATCACCGGAGACAGCGCAACGGTGGATAAAAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAgc                                                                                                                                         <  1:564926/300‑1 (MQ=17)
                                                                                                                                                                             aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGTTGTCAGCGACGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCa                                                                                                                               >  2:61348/1‑301 (MQ=255)
                                                                                                                                                                                      ggtggatacggcaacatgttgaaTGGTGATGTGGCTACGGCAATATACAACGGCAAAACCATCGTTGACGGCAAGGATGCTACCGG‑TACGGAAATCAACGGTAATAACTGAAATGTTTTACAGAACGGCCGACT‑GCCGGGATCTACACCTCTATACACTCTTTCCCTACACGACGCTCTTCCGATGTTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCa                                                                                                                      <  1:136891/280‑1 (MQ=9)
                                                                                                                                                                                                 aagccAGGTTAATGGGGGTGCGGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCAACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGTTCTGGAGGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAg                                                                                                           <  1:138754/298‑1 (MQ=255)
                                                                                                                                                                                                    taggATTAATGGTGATGACGCAACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAAGCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                        <  1:44675/298‑1 (MQ=255)
                                                                                                                                                                                                    ccAGATTAATGGTGATGACTCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCAGAta                                                                                                        >  2:608931/1‑300 (MQ=255)
                                                                                                                                                                                                      aGATTAATGGTGATGAAGCCAGGGAAAACAATAACGGGAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCGGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGCCCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGCATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                      <  1:426146/301‑1 (MQ=255)
                                                                                                                                                                                                       gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaga                                                                                                     >  2:44728/1‑299 (MQ=255)
                                                                                                                                                                                                       aaTTAATGGTGATGACGCTACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                     <  2:320786/300‑1 (MQ=255)
                                                                                                                                                                                                                       gCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATTACACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACgg                                                                                     <  2:384803/301‑1 (MQ=255)
                                                                                                                                                                                                                       acgctggCAAAGAACAGCGGCAAATACGCATGTGACGGCAAGGAGTACGACGG‑TAGGGAAAGCAAGGGGAAGAATGTAAAAGGTAGTGAGGAGGGCGATCGGGAGGTCAGCGGCGGCGGT‑CACGGTATTGATATAACTGTGGACGGCACAACGGTGGCCAAGAAGGGCACCACGACAGAC‑‑‑‑ACTGACCCGAAGGGAAGCGGTAGCGAGGTGGACGGTGACAAGGCGGTTGTCAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACgg                                                                                     <  1:61875/296‑1 (MQ=255)
                                                                                                                                                                                                                               aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCCGGTTGACGGCGACCAGGCCATTGTCCACAACGAAG‑CCGAGAGCACTATCACCAACGGCGGCACCGGCACACACATCAAcagctcagcagccccgactacacac                                                                             >  1:588242/1‑277 (MQ=11)
                                                                                                                                                                                                                                gccgtcaacggtggttttgccGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGCCAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATCACCCCACGGCAaacaaca                                                                            >  2:247360/20‑301 (MQ=14)
                                                                                                                                                                                                                                             aaaCCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGCGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACACCACGGCAAACAACACAGGCAAAaccac                                                               >  1:445476/1‑301 (MQ=255)
                                                                                                                                                                                                                                               ctcagtgTTGAGGGCATGGAGTCCGTCGG‑TACGTAAATCAACGGGAATAACGGAAAAGTTAGTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATAACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                                                             <  1:509857/295‑1 (MQ=17)
                                                                                                                                                                                                                                                ccaccGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCga                                                            >  2:510529/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      ttGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                      <  2:446040/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttggCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACggtggttttgccgttgttgtttg                                                   <  2:473965/299‑22 (MQ=255)
                                                                                                                                                                                                                                                          cGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaaa                                                  >  2:279181/1‑300 (MQ=255)
                                                                                                                                                                                                                                                           ggCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAAGCCAAGGCAAACAACAACggtggatttgccggtgttgtttgag                                                 >  1:473364/1‑278 (MQ=255)
                                                                                                                                                                                                                                                              aaaggggtgaCCGGTTACGGAATTCATAGGTAATTACTGAAAAGTTATTAAGGACGTGGATCTGTATGTCAGCGGCGGCGGT‑CACGTGTTTGATTTCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGTCTTTCTTTTATTCATACGGCGACCACCGCGATC‑‑‑‑TACACCTCTA‑TACGCTCTTTCCCTACACGACGCTCTTCCGATCTATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                               <  1:61280/292‑1 (MQ=2)
                                                                                                                                                                                                                                                                           g‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                 <  2:150506/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tcGGAATTAAACGTTAATAACGGAAAGGTGTTTCAGGACGGCGATCTGGATTTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTCGACGGCAAGGATTCCACCGGTACGGaaa                               <  2:16763/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGATAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGTCAAGGATTCCACCGTTACGGAAAttc                            >  2:240533/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         cGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCaaa                    >  2:631812/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                               aaaCGGAAAAGTTATTCAGGACGGGGTTCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACggg              <  1:630935/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  cGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaaa           >  2:214906/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                           tATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATCTCACGAGCGACAGCGAGGCGATGGGTGACTAGAGCAGCATGTCCGGC‑‑‑‑TCCGCGCCGTAGTCCATCGGAATCCAGGTTGACGGCGACCAGGCCGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCCCGCAGATTAACGGTGATGACGCCCCGGCACACACCCACGGCACAACCACCCTAAACGGCAACGATTCCACAGGTCCGAAAATAGCTACCAACAAAGCGAAggaaaagcac  >  1:558762/1‑293 (MQ=9)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                  
AGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGT‑CACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGG‑TACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGT‑CACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC‑‑‑‑ACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAG‑GCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAG  >  NC_000913/1466141‑1466738

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