Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,187 T→C 5.2% L737P (CTG→CCG)  mngB → alpha‑mannosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,1870TC 5.2% 520.2 / 11.8 154L737P (CTG→CCG) mngBalpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base T (38/108);  new base C (4/4);  total (42/112)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACG  >  NC_000913/770038‑770336
                                                                                                                                                     |                                                                                                                                                     
tGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCt                                                                                                                                                       >  2:1625767/1‑150 (MQ=255)
tGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCt                                                                                                                                                       <  1:706390/150‑1 (MQ=255)
tGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCt                                                                                                                                                       <  1:706387/150‑1 (MQ=255)
  aaTTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGc                                                                                                                                                     <  2:239992/150‑1 (MQ=255)
   aTTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCt                                                                                                                                                    <  2:282059/150‑1 (MQ=255)
   aTTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCt                                                                                                                                                    <  1:2179357/150‑1 (MQ=255)
    tttGAAGTCATCGGTGAAGAGAAGAAAAACTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGTTCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCtt                                                                                                                                                   <  2:419842/150‑1 (MQ=255)
      tGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCCGCTTTc                                                                                                                                                 >  1:296483/1‑150 (MQ=255)
           tCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTc                                                                                                                                            <  2:1580568/150‑1 (MQ=255)
           tCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTc                                                                                                                                            <  1:816753/150‑1 (MQ=255)
           tCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTc                                                                                                                                            <  1:1256640/150‑1 (MQ=255)
                 gTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTaa                                                                                                                                      >  2:869587/1‑150 (MQ=255)
                   gaagagaagaAAACCTTTGCCATTACGGTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTATAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGt                                                                                                                                    >  2:1613375/1‑150 (MQ=255)
                      gagatgaAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTTCTGGGCAATGTAGCTCTGCTTTTAAGGGCTGGGCGGCCTTCGGTAATTAAAATGCCAGTCCCGGACTCACAACTAGGTGGTCTGCTTTCTTGTCGCCTAAGTTTa                                                                                                                                 <  2:1884115/150‑1 (MQ=255)
                        gaagaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTAtt                                                                                                                               <  1:475002/150‑1 (MQ=255)
                        gaagaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTAtt                                                                                                                               <  1:45735/150‑1 (MQ=255)
                        gaagaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTAtt                                                                                                                               <  1:1900896/150‑1 (MQ=255)
                          agaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGa                                                                                                                             <  1:1741577/150‑1 (MQ=255)
                          agaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGa                                                                                                                             <  1:1775221/150‑1 (MQ=255)
                          agaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCATCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGa                                                                                                                             <  1:90982/150‑1 (MQ=255)
                           gaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAg                                                                                                                            >  2:1750369/1‑150 (MQ=255)
                           gaaAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAg                                                                                                                            <  1:224070/150‑1 (MQ=255)
                            aaaaCCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGt                                                                                                                           >  1:454331/1‑150 (MQ=255)
                             aaaCCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGtt                                                                                                                          >  1:1151/1‑150 (MQ=255)
                                   ttGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCg                                                                                                                    <  1:168128/150‑1 (MQ=255)
                                     gCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGt                                                                                                                  >  1:633968/1‑150 (MQ=255)
                                     gCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGt                                                                                                                  >  1:686184/1‑150 (MQ=255)
                                             gTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGACTTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAc                                                                                                          >  1:243834/1‑150 (MQ=255)
                                              ttGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAcc                                                                                                         <  2:670031/150‑1 (MQ=255)
                                              ttGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAcc                                                                                                         <  2:1974425/150‑1 (MQ=255)
                                                 ctgcGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAccgcc                                                                                                      <  1:1585994/150‑1 (MQ=255)
                                                           tggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAg                                                                                            <  1:36167/150‑1 (MQ=255)
                                                           tggGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAg                                                                                            <  2:1404034/150‑1 (MQ=255)
                                                             ggCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGt                                                                                                                  <  2:1171633/126‑1 (MQ=255)
                                                             ggCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGt                                                                                                                  >  1:1171633/1‑126 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            <  1:1533081/150‑1 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            >  2:404837/1‑150 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            >  2:353367/1‑150 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            >  1:697487/1‑150 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            >  1:1308887/1‑150 (MQ=255)
                                                                           agaagaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGc                                                                            <  1:2047647/150‑1 (MQ=255)
                                                                              agaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATg                                                                         <  2:598115/150‑1 (MQ=255)
                                                                               gaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATgg                                                                        <  1:94852/150‑1 (MQ=255)
                                                                               gaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATgg                                                                        <  2:1926340/150‑1 (MQ=255)
                                                                               gaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATgg                                                                        <  1:1979777/150‑1 (MQ=255)
                                                                               gaTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGGCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCGAAGTTTATTGAGTTATACCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCAAgg                                                                        <  1:248251/150‑1 (MQ=255)
                                                                                aTCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGc                                                                       <  1:562495/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTCCGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  2:1795199/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACTCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  2:61597/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  2:1917275/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  2:1917217/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  1:1990738/150‑1 (MQ=255)
                                                                                 tCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCt                                                                      <  2:1288031/150‑1 (MQ=255)
                                                                                  cTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTg                                                                     >  2:2368415/1‑150 (MQ=255)
                                                                                   tGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGa                                                                    <  1:396993/150‑1 (MQ=255)
                                                                                   tGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGa                                                                    <  1:1329315/150‑1 (MQ=255)
                                                                                   tGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGa                                                                    <  1:944857/150‑1 (MQ=255)
                                                                                   tGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGa                                                                    <  2:81593/150‑1 (MQ=255)
                                                                                     cTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACt                                                                  <  1:2086761/150‑1 (MQ=255)
                                                                                        ttAAGGCCTGGGCGGCCTTCGGGAATTAAAATGTCAGTCCCGGACTCACAACTACGTGGTCCGCTTTCTTGTCGCCTAAGTTTATTGTGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCa                                                               <  1:541994/150‑1 (MQ=255)
                                                                                        ttAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCCGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCa                                                               <  1:2024353/150‑1 (MQ=255)
                                                                                        ttAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCCGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCa                                                               <  2:296483/150‑1 (MQ=255)
                                                                                           aGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTTTTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                            <  2:1572958/150‑1 (MQ=255)
                                                                                           aGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                            <  2:2316499/150‑1 (MQ=255)
                                                                                           aGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                            <  1:1115841/150‑1 (MQ=255)
                                                                                           aGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                            <  1:2282242/150‑1 (MQ=255)
                                                                                           aGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTa                                                            <  1:235927/150‑1 (MQ=255)
                                                                                             gCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACa                                                          <  2:1372871/150‑1 (MQ=255)
                                                                                             gCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACa                                                          <  2:2276160/150‑1 (MQ=255)
                                                                                                  ggcggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGtt                                                     <  1:1147771/150‑1 (MQ=255)
                                                                                                  ggcggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGtt                                                     <  1:1809970/150‑1 (MQ=255)
                                                                                                   gcggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTa                                                    <  2:1080103/150‑1 (MQ=255)
                                                                                                    cggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTac                                                   <  2:2062782/150‑1 (MQ=255)
                                                                                                     ggcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTaca                                                  <  2:1403929/150‑1 (MQ=255)
                                                                                                       ccTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTacaac                                                <  2:1530649/150‑1 (MQ=255)
                                                                                                         ttCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTAcaacaa                                              >  1:1736631/1‑150 (MQ=255)
                                                                                                         ttCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTAcaacaa                                              >  1:260872/1‑150 (MQ=255)
                                                                                                            gggAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAAt                                           >  2:1219693/1‑150 (MQ=255)
                                                                                                              gAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATcc                                         <  2:2100167/150‑1 (MQ=255)
                                                                                                                 ttAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCAt                                      >  2:833322/1‑150 (MQ=255)
                                                                                                                  tAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATg                                     <  2:999419/150‑1 (MQ=255)
                                                                                                                  tAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATg                                     <  2:1326968/150‑1 (MQ=255)
                                                                                                                              cccGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGaa                         <  1:76619/150‑1 (MQ=255)
                                                                                                                               ccGGACTCCTAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTCTGGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACGCCAGTACAGTGTTACAACACAATCCCATGGGATGTGATGAAg                        <  2:1910560/150‑1 (MQ=255)
                                                                                                                               ccGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAg                        <  2:864613/150‑1 (MQ=255)
                                                                                                                               ccGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAg                        <  1:1300802/150‑1 (MQ=255)
                                                                                                                               ccGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAg                        <  1:1506353/150‑1 (MQ=255)
                                                                                                                                cGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGc                       >  1:2332119/1‑150 (MQ=255)
                                                                                                                                     tCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGACGCTcaac                  >  1:2155353/1‑150 (MQ=255)
                                                                                                                                     tCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTcaac                  <  1:415393/150‑1 (MQ=255)
                                                                                                                                       acaACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTcaacaa                >  2:2162552/1‑150 (MQ=255)
                                                                                                                                       acaACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTcaacaa                >  1:1076748/1‑150 (MQ=255)
                                                                                                                                       acaACTACGTGGTCTGCTTTCTTGTCGCCTAAGATTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGACTTACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTcaacaa                >  2:53320/1‑150 (MQ=255)
                                                                                                                                        caACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACaaa               <  2:615712/150‑1 (MQ=255)
                                                                                                                                         aaCTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAg              <  2:725812/150‑1 (MQ=255)
                                                                                                                                              cGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGa         <  2:228343/150‑1 (MQ=255)
                                                                                                                                                   tCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCaa    >  1:1534555/1‑150 (MQ=255)
                                                                                                                                                     tGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACg  >  1:1581314/1‑150 (MQ=255)
                                                                                                                                                     tGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACg  >  1:2240491/1‑150 (MQ=255)
                                                                                                                                                     tGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACg  >  2:2278199/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAATCCCATGGGATGTGATGAAGCTCAACAAAGCCGGATTCAACG  >  NC_000913/770038‑770336

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