Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,611,517 A→G 5.0% S21G (AGC→GGC)  hyfJ → putative processing element hydrogenase 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,611,5170AG 5.0% 306.8 / 4.2 99S21G (AGC→GGC) hyfJputative processing element hydrogenase 4
Reads supporting (aligned to +/- strand):  ref base A (49/45);  new base G (3/2);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAG  >  NC_000913/2611369‑2611663
                                                                                                                                                    |                                                                                                                                                  
aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg                                                                                                                                                   >  2:695871/1‑150 (MQ=255)
aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg                                                                                                                                                   >  2:544678/1‑150 (MQ=255)
    gggTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcgac                                                                                                                                               <  2:1541056/150‑1 (MQ=255)
     ggTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcgacg                                                                                                                                              <  1:1170700/150‑1 (MQ=255)
       ttAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCTACTGGcgacaaa                                                                                                                                            <  2:1785753/150‑4 (MQ=255)
        gaaCACCTGGTGGCGCGATGCCGACGACCCACGTCTCAATCGTATCGTTCAGCAACTCTTTCGCGTACACCGGCGGTTACATGACTGAACAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCTACTGGcgacaaac                                                                                                                                           <  1:295655/149‑5 (MQ=255)
         aaCACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTAGCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGAt                                                                                                                                          <  1:746038/150‑1 (MQ=255)
           cacCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGc                                                                                                                                        >  2:521778/1‑150 (MQ=255)
           cacCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGc                                                                                                                                        <  1:416578/150‑1 (MQ=255)
             ccTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCg                                                                                                                                      <  1:724575/150‑1 (MQ=255)
             ccTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCg                                                                                                                                      >  1:995311/1‑150 (MQ=255)
               gggTTGCGCGCTGCGGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGa                                                                                                                                    <  2:1120827/149‑1 (MQ=255)
                 gTTGCGCGATGCCGACGATCCACGGCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCTGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAac                                                                                                                                  <  2:1746449/150‑1 (MQ=255)
                  ttGCGCGTTGCCGACGTTCCCCGTCTCAATAGTCTCGTTCAGCTACTCTTTCGCGTACTCCGGGGGTTACCTGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAaca                                                                                                                                 <  2:54570/150‑1 (MQ=255)
                      gcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGc                                                                                                                             >  2:1130085/1‑150 (MQ=255)
                      gcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGAACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGc                                                                                                                             <  1:930070/150‑1 (MQ=255)
                             cgacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAgg                                                                                                                      <  2:1115805/150‑1 (MQ=255)
                             cgacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAgg                                                                                                                      >  1:1755351/1‑150 (MQ=255)
                             cgacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAgg                                                                                                                      >  1:1304645/1‑150 (MQ=255)
                             cgacgaTCCACGTCTCAATAGTATAGTTCAGCAACTCTTTAGAGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGATGCGAAAAAAGTTAGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTATCAgg                                                                                                                      >  2:1432748/1‑150 (MQ=255)
                              gacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGt                                                                                                                     <  2:451177/150‑1 (MQ=255)
                              gacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGt                                                                                                                     <  1:1768645/150‑1 (MQ=255)
                                cgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTaa                                                                                                                   <  2:193425/150‑1 (MQ=255)
                                 gaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAAt                                                                                                                  >  2:1261818/1‑150 (MQ=255)
                                   tCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGACCACAGCTCTCAGGTAATGt                                                                                                                >  2:883288/1‑150 (MQ=255)
                                      aCGTCGCAATAGTATCGTTCAGGAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTAtt                                                                                                             <  1:669908/150‑1 (MQ=255)
                                         tctcAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACt                                                                                                          <  2:1752406/150‑1 (MQ=255)
                                            cAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgc                                                                                                       >  2:8917/1‑150 (MQ=255)
                                               tAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGTGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgctag                                                                                                    <  1:89653/150‑1 (MQ=255)
                                               tAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgctag                                                                                                    <  2:558602/150‑1 (MQ=255)
                                                    tCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATc                                                                                               >  1:1053230/1‑150 (MQ=255)
                                                      gTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTGGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCgg                                                                                             >  2:1087287/1‑150 (MQ=255)
                                                      gTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCgg                                                                                             >  1:1824275/1‑150 (MQ=255)
                                                      gTTCAGCAACTCTTTCGCGTACTCCGGGCGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGCCGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCgg                                                                                             >  1:969814/1‑150 (MQ=255)
                                                       ttCAGCAACTCTTTGGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGc                                                                                            <  1:643554/150‑1 (MQ=255)
                                                        tCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGcc                                                                                           <  1:1589012/150‑1 (MQ=255)
                                                          aGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCAt                                                                                         >  1:1287313/1‑150 (MQ=255)
                                                           gCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATc                                                                                        >  2:1462540/1‑150 (MQ=255)
                                                           gCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATc                                                                                        >  1:411530/1‑150 (MQ=255)
                                                           gCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATc                                                                                        >  1:78957/1‑150 (MQ=255)
                                                             aaCTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCAc                                                                                      >  2:188722/1‑150 (MQ=255)
                                                              aCTCTTTCGCGTACTCCGGTGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACg                                                                                     <  1:357874/150‑1 (MQ=255)
                                                               ctctTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCGCTCAGGTAATGTATTACTCGCTAGCTATCGGCCAacactt                                                                                    >  2:1124408/1‑148 (MQ=255)
                                                                 ctTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTg                                                                                  <  1:1488108/150‑1 (MQ=255)
                                                                 ctTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTg                                                                                  >  1:1617390/1‑150 (MQ=255)
                                                                 ctTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTg                                                                                  >  2:1006752/1‑150 (MQ=255)
                                                                   ttCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGc                                                                                >  2:1206575/1‑150 (MQ=255)
                                                                   ttCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGc                                                                                >  2:1131430/1‑150 (MQ=255)
                                                                    tCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCg                                                                               >  1:1720988/1‑150 (MQ=255)
                                                                      gcgTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGtg                                                                             >  2:1602426/1‑150 (MQ=255)
                                                                        ttagtCCGGTGTTTACATGACTGTAGAGTGAGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTATTGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGtgat                                                                           <  1:409326/146‑1 (MQ=255)
                                                                            tCCGGGGGTTACTTGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGtgattgat                                                                       <  1:29897/150‑1 (MQ=255)
                                                                            tCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGTCGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGtgattgat                                                                       <  2:628708/150‑1 (MQ=255)
                                                                            tCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGtgattgat                                                                       <  1:1130085/150‑1 (MQ=255)
                                                                             ccGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTgattgatt                                                                      >  1:1455252/1‑150 (MQ=255)
                                                                              cGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGattgattg                                                                     <  1:1261818/150‑1 (MQ=255)
                                                                              cGGGGGGTACATGACTGAAGAGTGCGGGGAAATTGTATTCTGGACGCTGCGAAAAACGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGattgattg                                                                     <  2:1125698/150‑1 (MQ=255)
                                                                               gggggTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGt                                                                    >  2:1667294/1‑150 (MQ=255)
                                                                                ggggTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTc                                                                   >  2:917108/1‑150 (MQ=255)
                                                                                    ttACATTACTGAAGAGTGCGTGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTCTCGGCCATCACGTTGGCGTGATTGATTGTCTGaa                                                               <  2:252334/150‑1 (MQ=255)
                                                                                     tACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAAt                                                              >  2:1633848/1‑150 (MQ=255)
                                                                                     tACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGAGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCATGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAAt                                                              >  1:1753606/1‑150 (MQ=255)
                                                                                      aCATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATg                                                             <  2:796685/150‑1 (MQ=255)
                                                                                          gACTGAAGAGTGCGGGGAAATTTTTTTCTGGAAGCTGCGAAAAAAGTTTGTCGACAGTAGCGACGAGATGCCTGAACACAGCTCTCAGGTAATGTATTACTCTATAGCTATCGGCCATCACGTTTGCGTGATTGATTGTCTGAATGTCGc                                                         >  2:1077234/1‑150 (MQ=255)
                                                                                           aCTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGcc                                                        <  2:1723057/150‑1 (MQ=255)
                                                                                            cTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGAAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCt                                                       <  2:763181/150‑1 (MQ=255)
                                                                                              gAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTc                                                     >  1:521857/1‑150 (MQ=255)
                                                                                                 gagTGCGGTGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGc                                                  <  2:1016093/150‑1 (MQ=255)
                                                                                                   gTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTg                                                <  1:1029739/150‑1 (MQ=255)
                                                                                                    tGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGc                                               <  1:683614/150‑1 (MQ=255)
                                                                                                        gggAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCAc                                           <  1:695871/150‑1 (MQ=255)
                                                                                                         ggAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGTAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACCTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACt                                          >  1:1660990/1‑150 (MQ=255)
                                                                                                         ggAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACt                                          >  2:1425797/1‑150 (MQ=255)
                                                                                                         ggAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACt                                          >  2:1264827/1‑150 (MQ=255)
                                                                                                          gAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTg                                         >  1:621734/1‑150 (MQ=255)
                                                                                                          gAAATTGTTTTCTGGACGCTGCGAAAAAAGGTTTTCGCCCGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTATCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTAGCCTTCCGCTGCCCACTg                                         <  2:313240/150‑1 (MQ=255)
                                                                                                                gTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGaa                                   >  2:1191843/1‑150 (MQ=255)
                                                                                                                  tttCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATa                                 >  2:431458/1‑150 (MQ=255)
                                                                                                                   ttCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATAc                                <  1:1510570/150‑1 (MQ=255)
                                                                                                                    tCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACg                               >  1:815168/1‑150 (MQ=255)
                                                                                                                      tGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGaa                             >  1:579971/1‑150 (MQ=255)
                                                                                                                       ggACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAg                            >  1:959157/1‑150 (MQ=255)
                                                                                                                        gACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGgaga                           <  2:1720988/150‑1 (MQ=255)
                                                                                                                        gACGCGGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGa                           <  2:1608429/150‑1 (MQ=255)
                                                                                                                          cGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGAtt                         >  2:669037/1‑150 (MQ=255)
                                                                                                                             tGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGc                      <  2:506812/150‑1 (MQ=255)
                                                                                                                              gCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGGTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCt                     >  2:355050/1‑150 (MQ=255)
                                                                                                                               cGAAAAATGTTTGTCGCGAGTAGCGAGGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTTAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCtt                    <  1:1432748/150‑1 (MQ=255)
                                                                                                                               cGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCtt                    >  1:1761840/1‑150 (MQ=255)
                                                                                                                                gAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTg                   >  2:734722/1‑150 (MQ=255)
                                                                                                                                    gcatctcgTCGCTACTGGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACt               >  2:1316483/8‑150 (MQ=255)
                                                                                                                                     aaGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTg              <  2:929906/150‑1 (MQ=255)
                                                                                                                                       ttttGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCACGTAATGTATTACTCGCTAGCTATCGGCCATCACTTTGTCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGt            <  1:141296/149‑1 (MQ=255)
                                                                                                                                       gTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCCCTGACGGAATACGAAGATTGGCTTGCACTGGt            >  1:808002/1‑150 (MQ=255)
                                                                                                                                        tttGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTc           <  1:1462540/150‑1 (MQ=255)
                                                                                                                                         tcgTCGCTACTGGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCg          >  1:1035957/3‑150 (MQ=255)
                                                                                                                                           gTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGaa        >  1:1129635/1‑150 (MQ=255)
                                                                                                                                                cAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAgagga   >  1:69534/1‑150 (MQ=255)
                                                                                                                                                 aGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAgaggag  <  1:5767/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGCTATCGGCCATCACGTTGGCGTGATTGATTGTCTGAATGTCGCCTTCCGCTGCCCACTGACGGAATACGAAGATTGGCTTGCACTGGTCGAAGAGGAG  >  NC_000913/2611369‑2611663

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