Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,561,589 T→A 5.1% Q294L (CAG→CTG)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,561,5890TA 5.1% 730.6 / 6.7 235Q294L (CAG→CTG) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (95/120);  new base A (15/2);  total (113/122)
Fisher's exact test for biased strand distribution p-value = 5.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

ATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCC  >  NC_002947/5561441‑5561738
                                                                                                                                                    |                                                                                                                                                     
aTCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGAGATCGGCGTAGCGACCGCCCTg                                                                                                                                                      >  2:4391598/1‑150 (MQ=255)
 tCGCCACTGTCCGGCATCCAGATGCCGCAAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTATTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTgg                                                                                                                                                     <  1:397747/150‑1 (MQ=255)
 tCGCCAATGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCTGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTAGTAGCGACCGCCCTgg                                                                                                                                                     <  2:681793/150‑1 (MQ=255)
     cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGATATCGTCGTAGCGACCCCCCAGGGcga                                                                                                                                                 >  1:147712/1‑150 (MQ=255)
     cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCCTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGATATCGTCGTGGCGACCGCCCTGGGcga                                                                                                                                                 >  1:501723/1‑150 (MQ=255)
      aCTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat                                                                                                                                                <  2:3136478/150‑1 (MQ=255)
        tGTCCGTCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGATCGCCCAGGGcggtcg                                                                                                                                              <  1:353775/150‑1 (MQ=255)
         tggcGGCATCTGGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatcga                                                                                                                                             >  1:4446959/4‑150 (MQ=255)
         ggccGGCCTCCAGATGCCGCCAGTTGGCAATACGACCGCCGCTCGCCCCAGTGTGACCCGGGTCTTCAACTCCGTGGAGAAACCCGGGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGAGAGCGTCGTAGCGACCGCCCTGGGcgatcga                                                                                                                                             <  1:974210/148‑1 (MQ=255)
          tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatcgat                                                                                                                                            >  2:1216984/1‑150 (MQ=255)
            cGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGCTATCGTCGTAGCGACCGCCCAGGGCGATCGATTg                                                                                                                                          >  1:1264549/1‑150 (MQ=255)
            cGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCCAAGTCGGCGCCGCTAACCGCGTAGGGACCGCCCTGGGCGATCGATTg                                                                                                                                          >  1:3659085/1‑150 (MQ=255)
            cGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTg                                                                                                                                          <  1:3031947/150‑1 (MQ=255)
            cGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTg                                                                                                                                          <  2:1626008/150‑1 (MQ=255)
             ggCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTgg                                                                                                                                         >  1:1864043/1‑150 (MQ=255)
                    aGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCTGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGc                                                                                                                                  <  1:1768198/150‑1 (MQ=255)
                    aGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGGGCGCCCGGCCCAAGTCGGCGCCGCTATCGTCGTGGCGACCGCCCTGGGCGATCGATTGGCCTTCGc                                                                                                                                  >  1:830329/1‑150 (MQ=255)
                       tGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTg                                                                                                                               <  2:802549/150‑1 (MQ=255)
                           gccAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCCGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCAc                                                                                                                           <  1:4005959/150‑1 (MQ=255)
                              aGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                        <  2:1828069/150‑1 (MQ=255)
                              aGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGaa                                                                                                                        <  2:1861577/150‑1 (MQ=255)
                                   gCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACgg                                                                                                                   <  1:537790/150‑1 (MQ=255)
                                   gCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACgg                                                                                                                   <  2:3884053/150‑1 (MQ=255)
                                    cAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGc                                                                                                                  <  2:4019351/150‑1 (MQ=255)
                                    cAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGc                                                                                                                  <  2:1765063/150‑1 (MQ=255)
                                       tACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGGCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGCATTGCCTTCGCCTGGCACGAACACGGCGaa                                                                                                               >  1:445670/1‑150 (MQ=255)
                                          gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAACTCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAacac                                                                                                            <  2:3685925/150‑1 (MQ=255)
                                          gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAacac                                                                                                            >  2:4195494/1‑150 (MQ=255)
                                          gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAacac                                                                                                            >  1:2772350/1‑150 (MQ=255)
                                          gACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAacac                                                                                                            <  2:489948/150‑1 (MQ=255)
                                           aCCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAAcacc                                                                                                           >  1:4504103/1‑150 (MQ=255)
                                             caccGCTCGCCCCAGTGTGACCAGGGTCTACAAAACCGAGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAAcaccac                                                                                                         <  1:2199905/148‑1 (MQ=255)
                                             cTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAAcaccac                                                                                                         <  2:3138655/150‑1 (MQ=255)
                                              tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACg                                                                                                        >  2:2770403/1‑150 (MQ=255)
                                              tCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCAGGGCGATCGATTGGCCTTCGCCAGGCACGAACACCGCGAAcaccacc                                                                                                        >  1:2187650/1‑149 (MQ=255)
                                                     gCCCCAGTGTGACCAGGGTCTTCAGATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                 <  1:4241333/150‑1 (MQ=255)
                                                     cccccATTTTCACCAGTGTATTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCgtgt                                                                                                 <  1:2331766/149‑1 (MQ=255)
                                                      ccccAGTGTGACCCGGGTCTTCAAAGCCGTGGAGAACCCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGtgtg                                                                                                <  1:3212818/150‑1 (MQ=255)
                                                      ccccAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGGGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGtgtg                                                                                                <  2:2391858/150‑1 (MQ=255)
                                                         cAGTGTGACCAGGTTCTTAAAATCCGTGGTTAAACCTGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTgtggta                                                                                             <  1:1346315/150‑1 (MQ=255)
                                                         cAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCTGGGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTgtggta                                                                                             <  1:2418696/150‑1 (MQ=255)
                                                         cAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTgtggta                                                                                             >  2:507213/1‑150 (MQ=255)
                                                            tgtgACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTgtggtagtg                                                                                          <  2:2013564/150‑1 (MQ=255)
                                                              tgACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTgtggtagtggt                                                                                        <  2:2812961/150‑1 (MQ=255)
                                                                aCCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTaa                                                                                      <  1:1105128/150‑1 (MQ=255)
                                                                   aGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCa                                                                                   >  2:358614/1‑150 (MQ=255)
                                                                   aGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCCCCCGGCCACAGTCGGCGCCGATATCGTCGTAGCGCCCCCCCTGGGCGAGCGATTGGCCTTCGACTGGCACGCACACCGCGAACACCACCCCCGTCTCGTAGTGGGAACCa                                                                                   >  2:404824/1‑150 (MQ=255)
                                                                     ggTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTTGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACg                                                                                 >  1:2478428/1‑150 (MQ=255)
                                                                      gTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGc                                                                                <  2:754845/150‑1 (MQ=255)
                                                                         ttCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGAAATCGTCGTAGCGACCGCCCAGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGAGTGGTAGTGGTAACCACGCAGc                                                                             >  1:2141113/1‑150 (MQ=255)
                                                                            aaaTCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCg                                                                          >  2:3493600/1‑150 (MQ=255)
                                                                                ccctGGAGAAACCAGTGGCCGGGCGCGCCCCGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGa                                                                      <  2:2499461/147‑1 (MQ=255)
                                                                                  gTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAgg                                                                    <  2:2599223/150‑1 (MQ=255)
                                                                                   tGGAGACACCCGGGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGt                                                                   <  1:1161421/150‑1 (MQ=255)
                                                                                     gagaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCg                                                                 >  1:3067164/1‑150 (MQ=255)
                                                                                     gagaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCGTCGCCTGGCACGAACACGGCGGACACCACGCCCGTGTGGTAGGGGTAACCACGCAGCTCGCCGAGGGCg                                                                 >  1:1793480/1‑150 (MQ=255)
                                                                                       gaAACCTGGGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                               <  1:2851419/150‑1 (MQ=255)
                                                                                       gaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                               <  2:2689021/150‑1 (MQ=255)
                                                                                       gaAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                               <  1:2355170/150‑1 (MQ=255)
                                                                                       gaAACCAGTGGCCGGGCGCGCCCGGCCAAAGACGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGaa                                                               >  1:3699347/1‑150 (MQ=255)
                                                                                        aaaCCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCGCGCCGAGGTCGAAg                                                              >  2:3438701/1‑150 (MQ=255)
                                                                                         aaCCCGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGt                                                             <  1:2918346/150‑1 (MQ=255)
                                                                                          aCCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTa                                                            <  1:3516800/150‑1 (MQ=255)
                                                                                           ccAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGAGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTAc                                                           >  1:2496661/1‑150 (MQ=255)
                                                                                            cAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACa                                                          >  1:773580/1‑150 (MQ=255)
                                                                                              gTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGc                                                        <  1:3861142/150‑1 (MQ=255)
                                                                                                 gCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGc                                                     >  1:1065223/1‑150 (MQ=255)
                                                                                                 gCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTAGAAGTACAGCGGc                                                     <  1:4047940/150‑1 (MQ=255)
                                                                                                  ccGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGcc                                                    >  1:4343540/1‑149 (MQ=255)
                                                                                                         cgcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCgg                                             <  1:4287087/150‑1 (MQ=255)
                                                                                                          gcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCggg                                            <  2:3923247/150‑1 (MQ=255)
                                                                                                          gcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCggg                                            >  2:4407059/1‑150 (MQ=255)
                                                                                                          gcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCggg                                            <  2:1729163/150‑1 (MQ=255)
                                                                                                          gcCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCggg                                            <  2:1493439/150‑1 (MQ=255)
                                                                                                            ccGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCGTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTa                                          <  1:1627114/150‑1 (MQ=255)
                                                                                                             cGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAg                                         <  1:4202209/150‑1 (MQ=255)
                                                                                                               gCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAgcg                                       >  1:1172970/1‑150 (MQ=255)
                                                                                                                 cAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGcgcg                                     >  1:1010102/1‑150 (MQ=255)
                                                                                                                  aaaGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGa                                    >  2:2911042/1‑150 (MQ=255)
                                                                                                                          cgcCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGc                            <  1:3941746/150‑1 (MQ=255)
                                                                                                                               atatCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCa                       <  2:2453466/150‑1 (MQ=255)
                                                                                                                               atatCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCa                       >  1:3396931/1‑150 (MQ=255)
                                                                                                                               atatCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCa                       <  1:2030406/150‑1 (MQ=255)
                                                                                                                                 atCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGc                     <  2:1301708/150‑1 (MQ=255)
                                                                                                                                 atCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGc                     <  2:3762457/150‑1 (MQ=255)
                                                                                                                                   cgtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGa                   <  1:2885116/150‑1 (MQ=255)
                                                                                                                                   cgtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGa                   <  2:4504103/150‑1 (MQ=255)
                                                                                                                                    gtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGAt                  >  2:1179937/1‑150 (MQ=255)
                                                                                                                                    gtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGAt                  >  1:56471/1‑150 (MQ=255)
                                                                                                                                       gtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGc               >  2:2040852/1‑150 (MQ=255)
                                                                                                                                        aatcGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGcc              >  1:2003020/4‑150 (MQ=255)
                                                                                                                                            gACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCGCTGATCGcctggc          <  1:3744164/150‑5 (MQ=255)
                                                                                                                                            gACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagca          <  2:3305850/150‑1 (MQ=255)
                                                                                                                                             aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagcag         >  2:3866671/1‑150 (MQ=255)
                                                                                                                                             aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagcag         >  2:582888/1‑150 (MQ=255)
                                                                                                                                                gCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTc      <  2:1547526/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCa     <  2:219810/150‑1 (MQ=255)
                                                                                                                                                  ccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCAt    >  1:4366308/1‑150 (MQ=255)
                                                                                                                                                  ccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCAt    >  1:2950694/1‑150 (MQ=255)
                                                                                                                                                   cTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATc   <  2:605326/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCC  >  NC_002947/5561441‑5561738

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