Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,847,936 A→G 5.1% intergenic (+30/+61) yfhL → / ← mutM 4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,847,9360AG 5.1% 561.6 / 11.5 175intergenic (+30/+61)yfhL/mutM4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (87/79);  new base G (5/4);  total (92/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 182 total aligned reads displayed.

TACGATGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGC  >  NC_002947/5847787‑5848084
                                                                                                                                                     |                                                                                                                                                    
tACGATGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATACCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTga                                                                                                                                                      <  2:1457861/150‑1 (MQ=255)
   gATGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGagag                                                                                                                                                   >  2:2181706/1‑150 (MQ=255)
     tGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCg                                                                                                                                                 >  1:308983/1‑150 (MQ=255)
     tGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCg                                                                                                                                                 >  1:707973/1‑150 (MQ=255)
       aGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCggg                                                                                                                                               >  2:3149646/1‑150 (MQ=255)
       aGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCggg                                                                                                                                               <  2:3003933/150‑1 (MQ=255)
         ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTga                                                                                                                                                      <  1:2950514/141‑1 (MQ=255)
         ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTga                                                                                                                                                      >  2:2950514/1‑141 (MQ=255)
         ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCa                                                                                                                                             >  1:4456061/1‑150 (MQ=255)
         ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCa                                                                                                                                             >  2:1895753/1‑150 (MQ=255)
         ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCa                                                                                                                                             >  2:4551328/1‑150 (MQ=255)
           gCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATg                                                                                                                                           >  2:4670252/1‑150 (MQ=255)
               tGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTg                                                                                                                                       >  1:274285/1‑150 (MQ=255)
                gCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGc                                                                                                                                      <  1:1097577/150‑1 (MQ=255)
                 ccAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCt                                                                                                                                     >  1:489300/1‑150 (MQ=255)
                  cagcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTc                                                                                                                                    <  2:3101186/150‑1 (MQ=255)
                  cagcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTc                                                                                                                                    <  2:3336364/150‑1 (MQ=255)
                    gcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTccc                                                                                                                                  <  1:4547770/150‑1 (MQ=255)
                       ggTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCcaca                                                                                                                               <  2:94703/150‑1 (MQ=255)
                        gTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAg                                                                                                                              >  2:456147/1‑150 (MQ=255)
                         tCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGt                                                                                                                             >  2:856303/1‑150 (MQ=255)
                           tGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGatt                                                                                                                           <  2:4240203/150‑1 (MQ=255)
                            gCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                          <  2:3825959/150‑1 (MQ=255)
                                 gTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGttt                                                                                                                     >  2:3920230/1‑150 (MQ=255)
                                  tCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTc                                                                                                                    <  2:4376389/150‑1 (MQ=255)
                                     aCTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCt                                                                                                                 >  1:1735178/1‑150 (MQ=255)
                                          aTCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAAc                                                                                                            <  1:3739415/150‑1 (MQ=255)
                                           tCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACa                                                                                                           <  1:1828939/150‑1 (MQ=255)
                                           tCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACa                                                                                                           >  2:1716942/1‑150 (MQ=255)
                                             ccGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAgg                                                                                                         >  1:2519022/1‑150 (MQ=255)
                                                  ggATGAAGCGCATCAGGATACCGAAGTACAGCTGAGGGCCATTTATCGCCGGATCACTGGCAAGGGCTGTTCATTGCCCTTCCCTGCTCGGTCGCTGTGAGATCGGGGATGCGTGCTGACACAGTTTAGGTTTCTCTTCAACAGGTTCGc                                                                                                    <  2:1392781/150‑1 (MQ=255)
                                                   gATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCa                                                                                                   <  1:1505573/150‑1 (MQ=255)
                                                     tGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGc                                                                                                 <  2:789579/150‑1 (MQ=255)
                                                     tGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGc                                                                                                 <  1:3576864/150‑1 (MQ=255)
                                                         gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                             >  2:3473951/1‑150 (MQ=255)
                                                         gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                             >  2:2882988/1‑150 (MQ=255)
                                                              tCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                        >  2:83908/1‑150 (MQ=255)
                                                              tCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                        >  1:1919504/1‑150 (MQ=255)
                                                                 ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCt                                                                                                                                     >  2:2666937/1‑102 (MQ=255)
                                                                 ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATc                                                                                     >  1:3250121/1‑150 (MQ=255)
                                                                   aaaCCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCATGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGt                                                                                   <  2:2259525/150‑1 (MQ=255)
                                                                      ccGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                <  2:5107694/150‑1 (MQ=255)
                                                                       cGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGc                                                                               >  2:847604/1‑150 (MQ=255)
                                                                           gaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcg                                                                           >  1:3030890/1‑150 (MQ=255)
                                                                            aaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                             <  2:4186122/115‑1 (MQ=255)
                                                                            aaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc                                                                          >  2:203859/1‑150 (MQ=255)
                                                                            aaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc                                                                          >  2:5006899/1‑150 (MQ=255)
                                                                              cAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCgg                                                                        >  2:3744059/1‑150 (MQ=255)
                                                                              cAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCgg                                                                        <  1:4670252/150‑1 (MQ=255)
                                                                               aGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg                                                                       >  1:4255076/1‑150 (MQ=255)
                                                                               aGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGGTTCATCGTTGGCAGCGCggg                                                                       >  2:3087376/1‑150 (MQ=255)
                                                                                 cTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCa                                                                     <  1:1730803/150‑1 (MQ=255)
                                                                                   gATGGGCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCGACAGTTTAGGTTTCGCTTCATCAGGTTCGCAGCGAACGTTTCATCGTTGGCTGCGCGGGCAGt                                                                   >  2:1671460/1‑150 (MQ=255)
                                                                                   gATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGt                                                                   <  1:4156168/150‑1 (MQ=255)
                                                                                     tGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTac                                                                 <  2:4456061/150‑1 (MQ=255)
                                                                                       gCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTacac                                                               >  1:1409890/1‑150 (MQ=255)
                                                                                        ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACg                                                              >  1:5088707/1‑150 (MQ=255)
                                                                                         cAAGTATCGCCGGATCACTGGCACGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCGCCCACAGTTTAGGTTTCGCTTCCACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGc                                                             <  2:3538773/150‑1 (MQ=255)
                                                                                         cAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTGCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGTCAGTACACGc                                                             <  2:3538858/150‑1 (MQ=255)
                                                                                           aGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                           <  2:2481307/150‑1 (MQ=255)
                                                                                           aGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                           <  2:3889208/150‑1 (MQ=255)
                                                                                           aGTATCGCCGGATCACTGGCAAGGGCTGATCAGTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTc                                                           <  1:4140847/150‑1 (MQ=255)
                                                                                            gTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCg                                                          <  1:3170792/150‑1 (MQ=255)
                                                                                                    ggATCACTGGCAAGGGCTGACCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGc                                                  >  1:516006/1‑150 (MQ=255)
                                                                                                     gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGcc                                                 >  2:1519419/1‑150 (MQ=255)
                                                                                                     gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGcc                                                 >  1:4224885/1‑150 (MQ=255)
                                                                                                     gATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGcc                                                 >  2:1519347/1‑150 (MQ=255)
                                                                                                      aTCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGccc                                                <  2:3250121/150‑1 (MQ=255)
                                                                                                         aCTGGCAAGGGGTGATCATTGCCCTTCCCTGCACGGTCGCTGTGAGTGCGGGCATGGCTGCTCCGACAGTTTAGGTTTCGCTTCAAGAGGTTCGCAGCGAAGGTTTCTTGGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCGAg                                             >  2:4157738/1‑150 (MQ=255)
                                                                                                             gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTc                                         >  1:3421549/1‑150 (MQ=255)
                                                                                                             gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTc                                         <  1:2647054/150‑1 (MQ=255)
                                                                                                             gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTc                                         >  1:5067440/1‑150 (MQ=255)
                                                                                                              cAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCCCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCa                                        <  2:3382132/150‑1 (MQ=255)
                                                                                                               aaGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCCGTACACGCTCGCCCTCTGCCCCAGCTTCAc                                       >  1:3517083/1‑150 (MQ=255)
                                                                                                                   gCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcg                                   >  1:480698/1‑150 (MQ=255)
                                                                                                                   gCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcg                                   <  2:1151770/150‑1 (MQ=255)
                                                                                                                    cTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                  <  1:1058390/150‑1 (MQ=255)
                                                                                                                    cTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                  >  2:2913034/1‑150 (MQ=255)
                                                                                                                    cTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                  <  2:4493678/150‑1 (MQ=255)
                                                                                                                              cccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGt                        <  1:3653692/150‑1 (MQ=255)
                                                                                                                                 ttCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGcc                     >  1:3991832/1‑150 (MQ=255)
                                                                                                                                    ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                  <  1:4288687/150‑1 (MQ=255)
                                                                                                                                    ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                  <  1:2181706/150‑1 (MQ=255)
                                                                                                                                     cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCAt                 >  1:1373049/1‑150 (MQ=255)
                                                                                                                                      tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATa                <  1:4491112/150‑1 (MQ=255)
                                                                                                                                      tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATa                <  2:707973/150‑1 (MQ=255)
                                                                                                                                      tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATa                <  1:4954595/150‑1 (MQ=255)
                                                                                                                                         aCGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCt             >  1:420968/1‑150 (MQ=255)
                                                                                                                                            taaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                >  1:2530302/5‑80 (MQ=25)
                                                                                                                                            taaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                <  2:2530302/76‑1 (MQ=25)
                                                                                                                                            taaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                <  2:2530523/76‑1 (MQ=25)
                                                                                                                                              cccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCAc        <  2:274285/150‑1 (MQ=255)
                                                                                                                                               ccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                             >  1:2240180/1‑64 (MQ=255)
                                                                                                                                               ccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc                                                                                             <  2:2240180/64‑1 (MQ=255)
                                                                                                                                               acTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACg       >  2:1936820/2‑150 (MQ=255)
                                                                                                                                                cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg      >  2:2862580/1‑150 (MQ=255)
                                                                                                                                                  gtgAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCt    <  1:96914/150‑1 (MQ=255)
                                                                                                                                                  gtgAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCt    <  1:3804826/150‑1 (MQ=255)
                                                                                                                                                   tgAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTg   >  1:4652700/1‑150 (MQ=255)
                                                                                                                                                    gagaGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGc  <  1:648505/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACGATGAGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGC  >  NC_002947/5847787‑5848084

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