Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,869 C→T 5.8% P392P (CCC→CCT ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,8690CT 5.8% 363.6 / 7.4 122P392P (CCC→CCTybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (62/51);  new base T (5/2);  total (68/54)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 127 total aligned reads displayed.

CTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTG  >  NC_002947/4479720‑4480017
                                                                                                                                                     |                                                                                                                                                     
    tgGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCa                                                                                                                                                   >  1:246740/1‑150 (MQ=255)
      gAGNACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCAcc                                                                                                                                                 >  2:1505397/1‑150 (MQ=255)
        gCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCtt                                                                                                                                               >  1:2001562/1‑150 (MQ=255)
            aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGTGCAGCACCTGGcg                                                                                                                                                   >  1:403663/1‑141 (MQ=255)
            aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCGGGTGCAGCACCTGGCa                                                                                                                                                   <  2:403663/142‑1 (MQ=255)
                  gCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGGGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     >  2:1285852/1‑150 (MQ=255)
                  gCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCTGCAAGGTCATCCTCGTGCGAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     >  2:1248159/1‑150 (MQ=255)
                  gCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     <  2:897317/150‑1 (MQ=255)
                  gCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     >  1:684521/1‑150 (MQ=255)
                  gCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     >  2:470402/1‑150 (MQ=255)
                         gcCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACctct                                                                                                                             >  1:650314/1‑150 (MQ=255)
                            tGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCg                                                                                                                                                      >  2:1736259/1‑123 (MQ=255)
                            tGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTaga                                                                                                                          <  1:1867108/150‑1 (MQ=255)
                             gCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTagag                                                                                                                         <  2:253430/150‑1 (MQ=255)
                                  cGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGATACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAggagga                                                                                                                    >  1:60523/1‑150 (MQ=255)
                                    gAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAAt                                                                                                                  <  2:603542/150‑1 (MQ=255)
                                     aaCTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGccc                                                                                                                                         <  2:342516/127‑1 (MQ=255)
                                     aaCTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGccc                                                                                                                                         >  1:342516/1‑127 (MQ=255)
                                     aaCTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATg                                                                                                                 >  2:67513/1‑150 (MQ=255)
                                        tGATCGTTACTCACCACTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATc                                                                                                              >  2:580957/1‑150 (MQ=255)
                                          aTCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTg                                                                                                                                              >  2:1759255/1‑117 (MQ=255)
                                          aTCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTg                                                                                                                                              <  1:1759255/117‑1 (MQ=255)
                                                aCTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGTCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACgg                                                                                                      >  2:1100383/1‑150 (MQ=255)
                                                  tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtc                                                                                                    <  2:106784/150‑1 (MQ=255)
                                                  tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtc                                                                                                    >  1:1367194/1‑150 (MQ=255)
                                                    aCCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGc                                                                                                                                                    >  1:1820483/1‑101 (MQ=255)
                                                    aCCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGc                                                                                                                                                    <  2:1820483/101‑1 (MQ=255)
                                                         ttAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAg                                                                                             <  2:246740/150‑1 (MQ=255)
                                                              ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGt                                                                                                     >  1:1603047/1‑137 (MQ=255)
                                                              ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGt                                                                                                     <  2:1603047/137‑1 (MQ=255)
                                                              ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACt                                                                                        >  1:759638/1‑150 (MQ=255)
                                                               cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGt                                                                                                     <  2:680748/136‑1 (MQ=255)
                                                               cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                       >  2:712562/1‑150 (MQ=255)
                                                               cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                       >  2:1768665/1‑150 (MQ=255)
                                                               cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACAATAGCACACGTCTGAACTCAAGTGATGTATGAGCATGAACGGTCTCGCCAGTTACtt                                                                                       >  1:2377739/1‑150 (MQ=255)
                                                                tAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACct                                                                                                                               >  2:818353/1‑109 (MQ=255)
                                                                         cggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGcc                                                                             <  1:470402/150‑1 (MQ=255)
                                                                            ccATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTa                                                                          <  2:2001562/150‑1 (MQ=255)
                                                                               tGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGTCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCg                                                                       <  1:1100383/150‑1 (MQ=255)
                                                                               tGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGTCCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCg                                                                       <  2:2292986/150‑1 (MQ=255)
                                                                                  gCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCt                                                                                                                                        <  1:1034104/83‑1 (MQ=255)
                                                                                  gCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCt                                                                                                                                        >  2:1034104/1‑83 (MQ=255)
                                                                                    tACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGATGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAt                                                                  >  2:1186488/1‑150 (MQ=255)
                                                                                         aaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                             <  1:1505397/150‑1 (MQ=255)
                                                                                            gTTCGCTCAAAAGTATTATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                                          <  2:60523/150‑1 (MQ=255)
                                                                                             ttCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCt                                                         >  2:1334178/1‑150 (MQ=255)
                                                                                                gATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCAt                                                      >  1:469638/1‑150 (MQ=255)
                                                                                                 aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTCGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATc                                                     >  2:1186122/1‑150 (MQ=255)
                                                                                                 aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACTGTCTCGCCAGTTACTTCTGTATTGCCGTAGCGGCAATTCTTCTGCTCATc                                                     <  2:684521/150‑1 (MQ=255)
                                                                                                 aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATc                                                     >  1:660382/1‑150 (MQ=255)
                                                                                                    aaaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGcc                                                                             <  2:1892252/123‑1 (MQ=255)
                                                                                                    aaaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGcc                                                                             >  1:1892252/1‑123 (MQ=255)
                                                                                                    aaaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                  >  1:1737307/1‑150 (MQ=255)
                                                                                                      aaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGa                                                                                  >  2:941267/1‑116 (MQ=255)
                                                                                                      aaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGa                                                                                  <  1:941267/116‑1 (MQ=255)
                                                                                                      aaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCt                                                <  1:712562/150‑1 (MQ=255)
                                                                                                         cAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCAGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTa                                                                                          <  2:2223307/105‑1 (MQ=255)
                                                                                                         cAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCAGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTa                                                                                          >  1:2223307/1‑105 (MQ=255)
                                                                                                          aGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTgg                                            <  2:650314/150‑1 (MQ=255)
                                                                                                          aGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTgg                                            <  2:469638/150‑1 (MQ=255)
                                                                                                            gATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGc                                          >  2:1457330/1‑150 (MQ=255)
                                                                                                                  tGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     <  1:1922335/54‑1 (MQ=255)
                                                                                                                  tGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGac                                                                                                                                     >  2:1922335/1‑54 (MQ=255)
                                                                                                                  tGCCGCAAGGTCATACTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGt                                    >  2:2033217/1‑150 (MQ=255)
                                                                                                                    ccGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCa                                  <  2:1367194/150‑1 (MQ=255)
                                                                                                                    ccGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCATGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCa                                  >  1:238889/1‑150 (MQ=255)
                                                                                                                      gCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGc                                <  1:2136413/150‑1 (MQ=255)
                                                                                                                        aaGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCgt                              >  1:878130/1‑150 (MQ=255)
                                                                                                                           gTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGcc                                                                             <  2:1226365/100‑1 (MQ=255)
                                                                                                                           gTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGcc                                                                             >  1:1226365/1‑100 (MQ=255)
                                                                                                                            tCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTc                          <  2:1570886/150‑1 (MQ=255)
                                                                                                                              aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                        >  1:1801139/1‑150 (MQ=255)
                                                                                                                              aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                        >  2:1386427/1‑150 (MQ=255)
                                                                                                                                caaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGCGGAGGCATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTCTGGGCCATCGTCAGCGTGTTCCGGa                      >  2:785574/5‑150 (MQ=255)
                                                                                                                                caaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGa                      >  2:941183/5‑150 (MQ=255)
                                                                                                                                 cTCGTGCCAGGTGTTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAg                     >  1:1236467/1‑150 (MQ=255)
                                                                                                                                 cTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCAt                                                      <  1:536217/117‑1 (MQ=255)
                                                                                                                                 cTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTc                                   >  1:410087/1‑136 (MQ=255)
                                                                                                                                 cTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTc                                   <  2:410087/136‑1 (MQ=255)
                                                                                                                                  tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                  <  1:1181576/120‑1 (MQ=255)
                                                                                                                                  tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTgg                                            <  1:870381/126‑1 (MQ=255)
                                                                                                                                  tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc                    <  1:1334178/150‑1 (MQ=255)
                                                                                                                                  tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc                    <  1:1248159/150‑1 (MQ=255)
                                                                                                                                    gTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc                    <  1:387139/148‑1 (MQ=255)
                                                                                                                                    gTGCCAGGTGCTGCACCCGACACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGa                  >  1:1554891/1‑150 (MQ=255)
                                                                                                                                     tGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGAc                 >  1:2120260/1‑150 (MQ=255)
                                                                                                                                       ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAAGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAg                     >  2:720509/1‑144 (MQ=255)
                                                                                                                                       ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAAGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAg                     <  1:720509/144‑1 (MQ=255)
                                                                                                                                        cAGGTGCTGCACCCGGCACCTTGTGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCTTGTTCCGGAGCGACAAg              >  2:1184846/1‑150 (MQ=255)
                                                                                                                                         aGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGa             <  1:1186122/150‑1 (MQ=255)
                                                                                                                                          ggTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGAc            >  1:947825/1‑150 (MQ=255)
                                                                                                                                          ggTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGAc            >  1:1457471/1‑150 (MQ=255)
                                                                                                                                          ggTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGAc            >  1:1692156/1‑150 (MQ=255)
                                                                                                                                          ggTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACAGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGAc            >  2:1237390/1‑150 (MQ=255)
                                                                                                                                           gTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGAcc           >  2:329702/1‑150 (MQ=255)
                                                                                                                                                 cACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGcc     >  1:2452386/1‑150 (MQ=255)
                                                                                                                                                  aCCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                             <  2:1363114/93‑1 (MQ=255)
                                                                                                                                                    ccGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCt                                                <  2:2287087/104‑1 (MQ=255)
                                                                                                                                                    ccGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTg  >  1:1678845/1‑150 (MQ=255)
                                                                                                                                                    ccGGCACCTTGGGCCCTGACCACCTATCGCTGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTg  >  1:663754/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTG  >  NC_002947/4479720‑4480017

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