Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,216 T→C 6.9% Q73R (CAG→CGG) ‡ mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,2160TC 6.9% 887.0 / 36.0 274Q73R (CAG→CGG) ‡mobBmolybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base T (133/122);  new base C (11/8);  total (144/130)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 292 total aligned reads displayed.

GCCAGCCCCTCAACATCATTAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040946‑4041481
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                          
                catTAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgt                                                                                                                                                                                                                                                        >  2:440283/1‑275 (MQ=255)
                       atCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAg                                                                                                                                                                                                                                              >  1:89247/1‑278 (MQ=255)
                             tAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCCTGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGCTGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                                                                                                                                        >  1:208267/1‑278 (MQ=255)
                                tGCGACATCAAGCTTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAACCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCCGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCCGCTCTTCTTCGTCTAGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGACAACGATGGTTTGGGCCGACCCAgac                                                                                                                                                                                                                                            >  4:183774/1‑269 (MQ=255)
                                       tCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGcc                                                                                                                                                                                                                                                     >  4:511157/1‑255 (MQ=255)
                                       tCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGCTATCCATTCGACTCGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCACCCTTGCGCACCTccac                                                                                                                                                                                                                              >  1:296107/1‑275 (MQ=255)
                                                                 tACAGCAATACCAATCCTGTTTATAACTAATTTTTCAGGTCGATGACCGGCTCCATCGCGAAACCGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAAGTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCaa                                                                                                                                                                                                       <  4:471929/275‑1 (MQ=255)
                                                                         tAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                             >  2:640175/1‑261 (MQ=255)
                                                                             aTGCCTTGCTATCACTAATTATTCAGTTCAATGTCCGGCTCAATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTCTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGGCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                         <  4:640170/261‑1 (MQ=255)
                                                                                                  ttcAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTgctggc                                                                                                                                                                                                                                                                         <  2:407122/176‑1 (MQ=255)
                                                                                                                          gcgAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTgct                                                                                                                                                                                                                                                                            <  5:892442/149‑1 (MQ=255)
                                                                                                                                   cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCaa                                                                                                                                                                                                                                                                   <  8:181099/149‑1 (MQ=255)
                                                                                                                                   cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCaa                                                                                                                                                                                                                                                                   >  5:298446/1‑149 (MQ=255)
                                                                                                                                    acaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAAc                                                                                                                                                                                                                                                                  >  8:964232/1‑149 (MQ=255)
                                                                                                                                     caATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTgctggct                                                                                                                                                                                                                                                                        <  5:963635/142‑1 (MQ=255)
                                                                                                                                     caATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTgctggct                                                                                                                                                                                                                                                                        >  6:963635/1‑142 (MQ=255)
                                                                                                                                      aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                                                                                                                                                                <  7:58433/149‑1 (MQ=255)
                                                                                                                                        tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                              <  6:1466654/149‑1 (MQ=255)
                                                                                                                                                        aTGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgt                                                                                                                                                                                                                                                        <  5:737820/139‑1 (MQ=255)
                                                                                                                                                         tGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGc                                                                                                                                                                                                                                             >  6:211151/1‑149 (MQ=255)
                                                                                                                                                            ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCCAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                          <  5:462888/149‑1 (MQ=255)
                                                                                                                                                            ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                          >  8:711500/1‑149 (MQ=255)
                                                                                                                                                            ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                          >  7:743145/1‑149 (MQ=255)
                                                                                                                                                            ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                          >  6:1260314/1‑149 (MQ=255)
                                                                                                                                                            ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTATTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                                                                                                                                             <  7:496302/130‑1 (MQ=255)
                                                                                                                                                             tAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                             >  4:762122/1‑177 (MQ=255)
                                                                                                                                                              aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                                                                                                                                        <  7:1082419/149‑1 (MQ=255)
                                                                                                                                                               aaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcg                                                                                                                                                                                                                                       >  7:1308882/1‑149 (MQ=255)
                                                                                                                                                                 cccTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                         <  2:762128/177‑1 (MQ=255)
                                                                                                                                                                    ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCtggctggc                                                                                                                                                                                                                                                                     >  3:63681/1‑114 (MQ=255)
                                                                                                                                                                    ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCt                                                                                                                                                                                                                                  <  6:727199/149‑1 (MQ=255)
                                                                                                                                                                    ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATNCTTATTCA‑GCACTGAACTGATCCCCCTCGCGCTTAATGCCGGGGTCATTTTTTTCAAACGCGTAGTTTTTcatgtcctcc                                                                                                >  3:230580/1‑269 (MQ=255)
                                                                                                                                                                     tCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGctggctgg                                                                                                                                                                                                                                                                      >  1:42608/1‑112 (MQ=255)
                                                                                                                                                                        aCCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATa                                                                                                                                                                                                                              <  7:247537/149‑1 (MQ=255)
                                                                                                                                                                        aCCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGAGCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATa                                                                                                                                                                                                                              <  5:1232863/149‑1 (MQ=255)
                                                                                                                                                                          cAGAATCAAATCCGGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGGTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCa                                                                                                 <  2:36513/271‑1 (MQ=255)
                                                                                                                                                                           aGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCt                                                                                                                                                                                                                           >  7:1452928/1‑149 (MQ=255)
                                                                                                                                                                            gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                          <  5:1680184/133‑1 (MQ=255)
                                                                                                                                                                             aaTCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTAt                                                                                                                                                                                                                         >  5:1359713/1‑149 (MQ=255)
                                                                                                                                                                               tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtg                                                                                                                                                                                                                                                         <  5:1466118/115‑1 (MQ=255)
                                                                                                                                                                               tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                                                                                                                                       >  7:120950/1‑149 (MQ=255)
                                                                                                                                                                                 aaaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCttt                                                                                                                                                                                                                     <  5:529291/149‑1 (MQ=255)
                                                                                                                                                                                    tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGcc                                                                                                                                                                                                                  <  5:55887/149‑1 (MQ=255)
                                                                                                                                                                                       aGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                                                                                                                                                              <  2:611969/150‑1 (MQ=255)
                                                                                                                                                                                          cttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATa                                                                                                                                                                                                                              >  7:190468/2‑131 (MQ=255)
                                                                                                                                                                                                gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                      >  7:1670239/1‑149 (MQ=255)
                                                                                                                                                                                                gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                      <  7:784568/149‑1 (MQ=255)
                                                                                                                                                                                                 tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa                                                                                                                                                                                                     >  6:1413679/1‑149 (MQ=255)
                                                                                                                                                                                                   tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                                                                                                                                                                                                                 >  1:151058/1‑119 (MQ=255)
                                                                                                                                                                                                   tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc                                                                                                                                                                                                   >  5:1094038/1‑149 (MQ=255)
                                                                                                                                                                                                   tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc                                                                                                                                                                                                   >  6:111912/1‑149 (MQ=255)
                                                                                                                                                                                                   tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTat                                                                    >  2:300293/1‑275 (MQ=255)
                                                                                                                                                                                                    ccATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCTGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTatc                                                                   <  2:183777/275‑1 (MQ=255)
                                                                                                                                                                                                       ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa                                                                                                                                                                                                     <  8:845847/143‑1 (MQ=255)
                                                                                                                                                                                                         cGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATc                                                                                                                                                                                             <  8:1659407/149‑1 (MQ=255)
                                                                                                                                                                                                             ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg                                                                                                                                                                                         >  5:812645/1‑149 (MQ=255)
                                                                                                                                                                                                                 gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                                                                                                                                     >  6:349431/1‑149 (MQ=255)
                                                                                                                                                                                                                 gagaAACAGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                         <  8:1088145/129‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAt                                                                                                                                                                                  <  7:163575/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAt                                                                                                                                                                                  <  6:611029/149‑1 (MQ=255)
                                                                                                                                                                                                                      aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAtac                                                                                                                                                                                <  5:400805/149‑3 (MQ=255)
                                                                                                                                                                                                                      aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                                <  7:422844/149‑1 (MQ=255)
                                                                                                                                                                                                                      aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                                <  6:1359713/149‑1 (MQ=255)
                                                                                                                                                                                                                      aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                                <  5:682778/149‑1 (MQ=255)
                                                                                                                                                                                                                             gtccatgtcttctTCGTGTGGTGTTTCCAGTC‑TCAAGGCCCGTCGTTGGTTGCTGGCAACGGTGGGTGGTGCCGCGCCGGCCTTTCGCAGCTCAGATCTTTCTTTGCCCGGCTTAACAACCGCCATGTCATGGGGCGTATGCTTAATCATGCCGTGG‑CGGAGCCCTCTTGCGCATAATGCCGGGATCAGTTTTTGCGAAAGCGGAGTTGTTCCGGTACCACCCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCGTCGGGCAAg                                          <  4:271713/268‑1 (MQ=255)
                                                                                                                                                                                                                             aTCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCa                                                                                                                                                                         >  6:44928/1‑149 (MQ=255)
                                                                                                                                                                                                                             aTCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCa                                                                                                                                                                         >  8:1365106/1‑149 (MQ=255)
                                                                                                                                                                                                                               ccAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑Gc                                                                                                                                                                      <  7:711500/149‑1 (MQ=255)
                                                                                                                                                                                                                                     cttcttCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATcc                                                                                                                                                                                                  <  2:476397/116‑1 (MQ=255)
                                                                                                                                                                                                                                      ttcttcGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg                                                                                                                                                                                          <  7:945146/123‑1 (MQ=255)
                                                                                                                                                                                                                                      ttcttcGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGAc                                                                                                                                                               <  5:1021913/149‑1 (MQ=255)
                                                                                                                                                                                                                                        cttcGTCTGGTGTTTCC‑GTCATCAAGACCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACgg                                                                                                                                                             >  5:432422/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tcGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGAt                                                                                                                                                           <  5:211151/149‑1 (MQ=255)
                                                                                                                                                                                                                                          tcGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGAt                                                                                                                                                           <  7:1539871/149‑1 (MQ=255)
                                                                                                                                                                                                                                             tCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATccc                                                                                                                                                        >  8:35155/1‑149 (MQ=255)
                                                                                                                                                                                                                                               tGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCctc                                                                                                                                                      >  6:1283820/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 ggggtgtcaggcatcaAGGCCCTGCGGTGGTGACGGGCAACTATGGTTTGTGCCGGGCCAGCCTGGCGCAGGTCCATGCAGTCTTGTCCTGGCTTATCAACATCCAGATCATGGTGCGAATGCTTGATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATCGCGGGTGCAGTGTTTGCAAAAGCGTGGTTTTTCCGGTACCGGTCCACGCGGGAAAGGCGAGTAACGGTATCATCGTATTTCCTGCCATCGGACAGGCTCCTCTGGCGTATTGAc                        <  4:363061/263‑1 (MQ=255)
                                                                                                                                                                                                                                                   gtTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGgc                                                                                                                                                  >  8:23139/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGAt                                                                                                                                                           >  5:1551691/1‑138 (MQ=255)
                                                                                                                                                                                                                                                     ttCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGAt                                                                                                                                                           <  6:1551691/138‑1 (MQ=255)
                                                                                                                                                                                                                                                        c‑GTCATCAAGGCCAATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑CCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACaaa               <  2:57132/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           cctccagtgCCCATTGTTGCTTGCTGTCACCGATTGTTTTGCCCGCGCCGGCCTTTCGCAGCTTATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATTGTGCGTATGCTTACTCA‑GCCCTGGACGGATCCCTCTTGCGCATAATTCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATg             <  4:402195/267‑1 (MQ=255)
                                                                                                                                                                                                                                                              tcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGcc                                                                                                                                        <  8:120950/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               caAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCg                                                                                                                                       <  7:177592/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                acaaaccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                                                                                                                                      >  7:851657/6‑149 (MQ=255)
                                                                                                                                                                                                                                                                aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                                                                                                                                      <  8:555099/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                                                                                                                                      <  8:738765/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGAt                                                                                                                                                           >  8:465359/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                  ggCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGa                                                                                                                                    <  6:812645/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAt                                                                                                                                   <  8:675016/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAt                                                                                                                                   >  8:1318037/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  <  6:278977/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                         >  6:540518/1‑78 (MQ=37)
                                                                                                                                                                                                                                                                     ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCa                                                                                                                                 >  8:1643060/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                                                                                                                                                            <  8:1482274/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCCCGCGGCAAAGGCGAGTAACGGTATCCTCGTTTTTCCTGCCATCGGGCAAGCTCCTCTCGCGTATTCCCGTTGACAAATGCACCTTTATGa  >  4:600574/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          gTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTCGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttt                                                                                                                            <  7:1107717/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ttGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGttttt                                                                                                                           <  7:639125/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ttGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGttttt                                                                                                                           <  8:945588/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt                                                                                                                          <  6:1722681/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                          
GCCAGCCCCTCAACATCATTAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCC‑GTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCA‑GCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040946‑4041481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑