Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,525 T→C 5.5% T427A (ACC→GCC)  mgtE ← magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,077,5250TC 5.5% 218.4 / 2.4 73T427A (ACC→GCC) mgtEmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGACACTT  >  NC_002947/5077380‑5077674
                                                                                                                                                 |                                                                                                                                                     
gCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGTGTGCCAAGGTCATTGTGGTCAGCATGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATc                                                                                                                                                   <  1:31059/150‑1 (MQ=255)
gCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATc                                                                                                                                                   >  2:164751/1‑150 (MQ=255)
  cAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGc                                                                                                                                                 >  2:991950/1‑150 (MQ=255)
     cccAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGc                                                                                                                                              >  2:801412/1‑150 (MQ=255)
     cccAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGc                                                                                                                                              >  1:1569474/1‑150 (MQ=255)
           gAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGtcat                                                                                                                                        <  1:745220/150‑1 (MQ=255)
            aaGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGtcatc                                                                                                                                       <  1:739787/150‑1 (MQ=255)
                 gaagaaACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATcaccac                                                                                                                                  >  2:1189015/1‑150 (MQ=255)
                    gaaACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCAcccc                                                                                                                               >  2:194181/1‑150 (MQ=255)
                    gaaACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCAcccc                                                                                                                               <  1:991950/150‑1 (MQ=255)
                            cGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGAcc                                                                                                                       >  2:1304612/1‑150 (MQ=255)
                                gTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCt                                                                                                                   <  2:1427933/150‑1 (MQ=255)
                                gTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCt                                                                                                                   <  1:1189015/150‑1 (MQ=255)
                                        cagcagTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACa                                                                                                           <  2:1316189/150‑1 (MQ=255)
                                          gcagTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACTACCCCCAGCGACCAGCTGCCATACAGc                                                                                                         >  1:661464/1‑150 (MQ=255)
                                               gATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGaa                                                                                                    >  1:103144/1‑150 (MQ=255)
                                                    tcaCACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCCGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGGCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGccc                                                                                               >  1:317450/1‑149 (MQ=255)
                                                    tcaCACTCGCCCCCATCGCAGGGGCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCa                                                                                               <  2:51866/150‑1 (MQ=255)
                                                     cacaCTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCAc                                                                                              <  2:1197370/150‑1 (MQ=255)
                                                        acTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGc                                                                                           <  2:697752/150‑1 (MQ=255)
                                                          tCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAACTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCAc                                                                                         <  1:1050263/150‑1 (MQ=255)
                                                            gCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCa                                                                                       <  2:1394007/150‑1 (MQ=255)
                                                              ccccATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGtcatc                                                                                                                                       <  2:554742/100‑1 (MQ=255)
                                                              ccccATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGtcatc                                                                                                                                       >  1:554742/1‑100 (MQ=255)
                                                              ccccATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAAt                                                                                     >  2:69810/1‑150 (MQ=255)
                                                                ccATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGa                                                                                   >  1:1002579/1‑150 (MQ=255)
                                                                 cATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGAc                                                                                  >  2:1641373/1‑150 (MQ=255)
                                                                 cATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGAc                                                                                  >  1:960684/1‑150 (MQ=255)
                                                                             gcgcCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGccc                                                                           >  1:639519/1‑145 (MQ=255)
                                                                             gcgcCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGccc                                                                           <  2:639519/145‑1 (MQ=255)
                                                                               gcCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCACAAcccc                                                                    >  2:344475/1‑149 (MQ=255)
                                                                                ccccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTTCGCCCCAATCCAg                                                                   <  2:475228/150‑1 (MQ=255)
                                                                                 cccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAgg                                                                  >  1:524409/1‑150 (MQ=255)
                                                                                 cccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGCCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCAcc                                                                                        >  1:1326062/1‑128 (MQ=255)
                                                                                 cccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGCCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCAcc                                                                                        >  1:1325938/1‑128 (MQ=255)
                                                                                 cccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGCCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCAcc                                                                                        <  2:1326062/128‑1 (MQ=255)
                                                                                 cccAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGCCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCAcc                                                                                        <  2:1325938/128‑1 (MQ=255)
                                                                                     aGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGt                                                              <  1:164751/150‑1 (MQ=255)
                                                                                        ccgCCAAGATCATCGGGACCAGCACGTCCCTCAGCGCCGCCAGTAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGa                                                           <  2:676647/148‑1 (MQ=255)
                                                                                          gCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGATCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATc                                                         >  1:568080/1‑150 (MQ=255)
                                                                                             aaGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGc                                                      <  1:801412/150‑1 (MQ=255)
                                                                                             aaGGGCATGGGGATCAGCACGCCCATCAGCGCCGCCAATAGCTAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGc                                                      <  2:804855/150‑1 (MQ=255)
                                                                                               ggTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGa                                                    <  2:103144/150‑1 (MQ=255)
                                                                                                   aTCGGTATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAAt                                                                                     <  2:1585413/113‑1 (MQ=255)
                                                                                                   aTCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAAt                                                                                     >  1:1585413/1‑113 (MQ=255)
                                                                                                   aTCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGa                                                                                   >  2:477585/1‑115 (MQ=255)
                                                                                                   aTCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGa                                                                                   <  1:477585/115‑1 (MQ=255)
                                                                                                     cGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGc                                              <  1:1304612/150‑1 (MQ=255)
                                                                                                        gATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCaa                                           >  2:1533425/1‑150 (MQ=255)
                                                                                                            aGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGc                                                      <  2:623231/135‑1 (MQ=255)
                                                                                                            aGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGc                                                      >  1:623231/1‑135 (MQ=255)
                                                                                                                      tCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCa                             >  1:377611/1‑150 (MQ=255)
                                                                                                                           gccgccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTg                                    >  1:1084216/1‑138 (MQ=255)
                                                                                                                           gccgccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTg                                    <  2:1084216/138‑1 (MQ=255)
                                                                                                                            ccgccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAAc                                                                      <  2:1325051/103‑1 (MQ=255)
                                                                                                                            ccgccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAAc                                                                      >  1:1325051/1‑103 (MQ=255)
                                                                                                                               ccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATTAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTgg                    <  1:1641373/150‑1 (MQ=255)
                                                                                                                               ccAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTgg                    >  2:57506/1‑150 (MQ=255)
                                                                                                                                cAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACt                                                                                             <  1:1131508/76‑1 (MQ=255)
                                                                                                                                cAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACt                                                                                             >  2:1131508/1‑76 (MQ=255)
                                                                                                                                 aaTAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGc                       <  2:38851/145‑1 (MQ=255)
                                                                                                                                 aaTAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGc                       >  1:38851/1‑145 (MQ=255)
                                                                                                                                      cAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGcc             >  1:1432786/1‑150 (MQ=255)
                                                                                                                                      cAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGcc             >  1:1050277/1‑150 (MQ=255)
                                                                                                                                       aaGATCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCt            >  2:1033997/1‑150 (MQ=255)
                                                                                                                                          ttCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTgg          >  2:1062022/1‑149 (MQ=255)
                                                                                                                                          ttCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTgg          <  1:1062022/149‑1 (MQ=255)
                                                                                                                                            cAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGa       >  2:1400898/1‑150 (MQ=255)
                                                                                                                                             aaGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCa                                                        <  2:1368045/100‑1 (MQ=255)
                                                                                                                                             aaGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCa                                                        >  1:1368045/1‑100 (MQ=255)
                                                                                                                                              aGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCCACTCCTTGCCCATCAAGCGGCTGGTATTGCCTGGCGacc     >  2:654905/1‑149 (MQ=255)
                                                                                                                                               ggTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGacac    <  2:1569474/150‑1 (MQ=255)
                                                                                                                                               ggTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGacac    <  2:1050277/150‑1 (MQ=255)
                                                                                                                                                 tCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGACACtt  >  1:1034869/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCATCACACTCGCCCCCATCGCAGGGTCGCGCCCCAAGCGGGCCAAGGTCATCGGGATCAGCACGCCCATCAGCGCCGCCAATAACAAGTTCAAGGTCATCGCTGCGGTCATCACCACCCCCAGCGACCAGCTGCCATACAGCCAGAAAGCCACTGCACCAATGACTCCGCCCCAAACCAGGCCGTTGATCAGCGACACCGCCAACTCCTTGCGCATCAAGCGGCTGGTATTGCCTGGCGACACTT  >  NC_002947/5077380‑5077674

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