Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,030,765 A→G 6.1% Q10R (CAG→CGG)  mgtA → ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,030,7650AG 6.1% 288.8 / 6.4 98Q10R (CAG→CGG) mgtAATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (53/39);  new base G (3/3);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

TTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGAACAGG  >  NC_002947/3030620‑3030914
                                                                                                                                                 |                                                                                                                                                     
ttGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCgg                                                                                                                                                   <  1:1781910/150‑1 (MQ=255)
 tGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGc                                                                                                                                                  <  2:248991/150‑1 (MQ=255)
   tttgAGTAGAAAGTTCGCTGCCGACACAACTTACTACTCTCGGTAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCCTGCCGCTGCGGCGGCac                                                                                                                                                <  2:445227/147‑2 (MQ=255)
   cgcgAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTGAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCTGCGGCGGCac                                                                                                                                                <  2:1007853/150‑2 (MQ=255)
   cgcgAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAg                                                                                                                                                <  1:1609407/150‑1 (MQ=255)
      gAGTCGACAGCTAGCTGCAGACACACCTTCCTACTCTCGGCATCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAg                                                                                                                                             <  2:1123734/150‑1 (MQ=255)
      gAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCNGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCCTGCCGCTGCGGCGGCACGAg                                                                                                                                             >  1:445227/1‑150 (MQ=255)
      gAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCGCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGCTCCTGCCGCTGCGGCGGCacgaca                                                                                                                                            <  2:621188/150‑6 (MQ=255)
         tAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAg                                                                                                                                                      <  1:1287515/138‑1 (MQ=255)
         tAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAg                                                                                                                                                      >  2:1287515/1‑138 (MQ=255)
         tAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCACCGACCCGAAATGCCTTCGTGCACCGCCTTAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGcac                                                                                                                                          >  1:2090357/1‑150 (MQ=255)
          aGAAAGTTCGCGGCCGACACAACCTCCGAATCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCGTTATTCAGCGACCCGACATGCCTTCTTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGcact                                                                                                                                         <  2:506102/150‑2 (MQ=255)
           gAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAg                                                                                                                                        >  2:221146/1‑150 (MQ=255)
               gTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAAc                                                                                                                                    <  1:221146/150‑1 (MQ=255)
                   gCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                                >  2:855245/1‑150 (MQ=255)
                     tGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACt                                                                                                                              <  2:2348526/150‑1 (MQ=255)
                          acacAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTATGTGCCGTGGATCGTCAGCTGCAGGAGCACAGCAACATGACTGTCaa                                                                                                                         >  1:1965120/1‑150 (MQ=255)
                                cTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCg                                                                                                                   >  1:155492/1‑150 (MQ=255)
                                 ttCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGt                                                                                                                  >  2:1293927/1‑150 (MQ=255)
                                   ccGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTaa                                                                                                                <  1:2040195/150‑1 (MQ=255)
                                     gACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAaca                                                                                                              >  2:848457/1‑150 (MQ=255)
                                        tctcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCt                                                                                                           >  2:1046038/1‑150 (MQ=255)
                                          tcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCa                                                                                                         >  1:1815290/1‑150 (MQ=255)
                                              cAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCgg                                                                                                     <  1:1293927/150‑1 (MQ=255)
                                               aGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGc                                                                                                    >  2:362797/1‑150 (MQ=255)
                                                  ggACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCaa                                                                                                 >  1:1596386/1‑150 (MQ=255)
                                                   gacagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                                <  1:1240343/150‑1 (MQ=255)
                                                    acagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCAGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACAGTAACACCTCATCGGCCAAgg                                                                                               >  2:336631/1‑150 (MQ=255)
                                                      agacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGcc                                                                                             >  1:1788153/1‑150 (MQ=255)
                                                        acaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCATTCATTCAAAGCCAGAGTGAAACGTCTGTGTTTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGa                                                                                           >  1:1500146/1‑150 (MQ=255)
                                                          aaCTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAc                                                                                         >  1:1031189/1‑150 (MQ=255)
                                                            cTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGCCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAcgc                                                                                       <  2:542353/150‑1 (MQ=255)
                                                             tGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTaa                                                                                                                <  2:1237219/124‑1 (MQ=255)
                                                             tGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTaa                                                                                                                >  1:1237219/1‑124 (MQ=255)
                                                             tGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAcgcg                                                                                      >  2:97989/1‑150 (MQ=255)
                                                              ctcagAAAGAACCAGCATTATTCAGTGACCCGACATTCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGTAGCGGCAGGAGCACAGCAACATGACTTGCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCga                                                                                     <  2:82307/146‑1 (MQ=255)
                                                               tCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCgag                                                                                    >  2:1511619/1‑150 (MQ=255)
                                                                 ggAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCgagag                                                                                  >  1:82307/1‑149 (MQ=255)
                                                                 ggAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCga                                                                                     <  2:82681/147‑1 (MQ=255)
                                                                 ggAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCga                                                                                     >  1:82681/1‑147 (MQ=255)
                                                                  gAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                                <  1:608164/103‑1 (MQ=255)
                                                                  gAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                                >  2:608164/1‑103 (MQ=255)
                                                                  gAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCGAAGGCCGAACGCGAGAcc                                                                                 >  1:503895/1‑150 (MQ=255)
                                                                  gAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGAcc                                                                                 >  2:1283887/1‑150 (MQ=255)
                                                                   aaaGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGAccc                                                                                >  2:152142/1‑150 (MQ=255)
                                                                    aaGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCg                                                                               <  1:1824388/150‑1 (MQ=255)
                                                                    aaGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGCCCAAGGCCGAACGCGAGccccc                                                                               >  2:1025947/1‑149 (MQ=255)
                                                                       aaCCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCt                                                                            >  1:1796191/1‑150 (MQ=255)
                                                                       aaCCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCt                                                                            <  1:935621/150‑1 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACAcc                                                                                                            >  1:705766/1‑115 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACAcc                                                                                                            <  2:705766/115‑1 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                                >  1:175852/1‑127 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                                <  2:175852/127‑1 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTc                                                                         >  2:1185158/1‑150 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTc                                                                         >  2:199918/1‑150 (MQ=255)
                                                                          cAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCAGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTc                                                                         >  1:449729/1‑150 (MQ=255)
                                                                           aGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTcc                                                                        >  1:1483737/1‑150 (MQ=255)
                                                                           aGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTcc                                                                        >  1:1185369/1‑150 (MQ=255)
                                                                           aGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTcc                                                                        <  1:362797/150‑1 (MQ=255)
                                                                              attattCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCAcc                                                                     >  1:725758/1‑150 (MQ=255)
                                                                                     aGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGcgccg                                                              >  1:951674/1‑150 (MQ=255)
                                                                                     aGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGcgccg                                                              >  2:672409/1‑150 (MQ=255)
                                                                                      gCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCg                                                                               <  2:1080625/132‑1 (MQ=255)
                                                                                      gCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCg                                                                               >  1:1080625/1‑132 (MQ=255)
                                                                                             gACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCgaga                                                                                   <  2:305332/121‑1 (MQ=255)
                                                                                             gACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCgaga                                                                                   >  1:305332/1‑121 (MQ=255)
                                                                                                    cTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGcaca                                                                                                                                         <  2:2217605/60‑1 (MQ=255)
                                                                                                    cTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGcaca                                                                                                                                         >  1:2217605/1‑60 (MQ=255)
                                                                                                    cTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCa                                               >  1:700428/1‑150 (MQ=255)
                                                                                                    cTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCa                                               >  2:726234/1‑150 (MQ=255)
                                                                                                           cACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCt                                        >  1:490757/1‑150 (MQ=255)
                                                                                                              cGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGc                                     <  1:97989/150‑1 (MQ=255)
                                                                                                               gCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCg                                    <  1:1283887/150‑1 (MQ=255)
                                                                                                                ccTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCgg                                   <  1:672409/150‑1 (MQ=255)
                                                                                                                    aGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCAc                               <  2:1788153/150‑1 (MQ=255)
                                                                                                                     gTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCAcc                              <  2:2090357/150‑1 (MQ=255)
                                                                                                                        aaaCGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTc                           >  1:1530241/1‑150 (MQ=255)
                                                                                                                         aaCGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCg                          <  1:578025/150‑1 (MQ=255)
                                                                                                                          aCGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGc                         <  2:1031189/150‑1 (MQ=255)
                                                                                                                              cTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAAc                     <  1:1046038/150‑1 (MQ=255)
                                                                                                                                 tgtgtCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTc                  >  1:1354289/1‑150 (MQ=255)
                                                                                                                                    gtCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGAt               <  2:1796191/150‑1 (MQ=255)
                                                                                                                                    ctcctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGGAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGAt               >  2:1784197/5‑150 (MQ=255)
                                                                                                                                     tcctGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGaa                                                      >  1:1191605/4‑110 (MQ=255)
                                                                                                                                     tcctGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCAGCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGaa                                                      <  2:1191605/107‑1 (MQ=255)
                                                                                                                                     tCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATTACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATg              >  1:680378/1‑150 (MQ=255)
                                                                                                                                      cGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGc             <  2:1530241/150‑1 (MQ=255)
                                                                                                                                       gTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAAcc                    >  2:2085531/1‑142 (MQ=255)
                                                                                                                                       gTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAAcc                    <  1:2085531/142‑1 (MQ=255)
                                                                                                                                       gTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGcc            >  1:183127/1‑150 (MQ=255)
                                                                                                                                       ctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGcc            >  2:1359783/2‑150 (MQ=255)
                                                                                                                                        tGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCa           <  1:848457/150‑1 (MQ=255)
                                                                                                                                          ccgccgCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGc         <  1:199918/150‑1 (MQ=255)
                                                                                                                                           cgccgcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCg        >  2:1867261/1‑150 (MQ=255)
                                                                                                                                           cgccgcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACAATCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCg        >  1:1363827/1‑150 (MQ=255)
                                                                                                                                               gcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAgg                                                                                               <  1:1656008/59‑1 (MQ=255)
                                                                                                                                               gcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAgg                                                                                               >  2:1656008/1‑59 (MQ=255)
                                                                                                                                               gcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCACTGTCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGAACa    >  1:22505/1‑150 (MQ=255)
                                                                                                                                                 aGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAAGCGAGACCCGCCTGTCCCCCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCCTTGCCAGCGAACAgg  <  2:825750/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGAACAGG  >  NC_002947/3030620‑3030914

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