Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,031,138 C→G 11.6% I134M (ATC→ATG mgtA → ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,031,1380CG11.6% 231.0 / 17.5 76I134M (ATC→ATGmgtAATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (42/25);  new base G (4/5);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCTGGTGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAG  >  NC_002947/3030994‑3031286
                                                                                                                                                |                                                                                                                                                    
ccTGGTGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGAc                                                                                                                                                 <  2:325488/150‑1 (MQ=255)
    gTGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATg                                                                                                                                             <  2:785600/150‑1 (MQ=255)
    gTGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATg                                                                                                                                             <  2:644565/150‑1 (MQ=255)
     tGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGACCTTgg                                                                                                                                            <  1:347771/150‑1 (MQ=255)
        aaCTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGACCTTGGTCa                                                                                                                                         >  1:905112/1‑150 (MQ=255)
          cTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGc                                                                                                                                       >  1:418824/1‑150 (MQ=255)
              tCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGAccttg                                                                                                                                             >  1:743134/1‑137 (MQ=255)
              tCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGACCTTg                                                                                                                                             <  2:743134/140‑1 (MQ=255)
              tCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTca                                                                                                                                   >  1:414681/1‑150 (MQ=255)
              tCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTca                                                                                                                                   >  2:786799/1‑150 (MQ=255)
              tCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTca                                                                                                                                   >  2:121648/1‑150 (MQ=255)
               cAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGACCTTGGTCAggtcggc                                                                                                                                  <  2:905112/150‑7 (MQ=38)
                aaGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTcagc                                                                                                                                 <  2:905251/150‑1 (MQ=255)
                aaGGCCCTGAACAACCCATTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTcagc                                                                                                                                 >  1:510397/1‑150 (MQ=255)
                 aGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCTAGCTGGCCAAAGTCATCAGCATGATGACCATGGTCAGCCTcagca                                                                                                                                <  1:1247688/150‑1 (MQ=255)
                   gCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCagcagc                                                                                                                              <  1:759085/150‑1 (MQ=255)
                           caacCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCg                                                                                                                      <  2:414681/150‑1 (MQ=255)
                            aacCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGa                                                                                                                     >  2:1305277/1‑150 (MQ=255)
                             acCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGAt                                                                                                                    <  1:1305277/150‑1 (MQ=255)
                              cccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGAtt                                                                                                                   >  1:1137418/1‑150 (MQ=255)
                               ccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATTGTCGGCCTCAGCAGCCTCCTGCGATTc                                                                                                                  >  2:233423/1‑150 (MQ=255)
                               ccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTc                                                                                                                  >  1:863006/1‑150 (MQ=255)
                               ccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTc                                                                                                                  >  1:1122274/1‑150 (MQ=255)
                                cGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCt                                                                                                                 >  1:898773/1‑150 (MQ=255)
                                 gTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATg                                                                                                                                             <  1:309466/121‑1 (MQ=255)
                                 gTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATg                                                                                                                                             >  2:309466/1‑121 (MQ=255)
                                 gTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTg                                                                                                                >  1:626284/1‑150 (MQ=255)
                                     aTCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                            >  2:1127281/1‑150 (MQ=255)
                                     aTCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                            <  1:1226348/150‑1 (MQ=255)
                                     aTCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                            >  1:1406958/1‑150 (MQ=255)
                                       cTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGa                                                                                                          >  2:1256567/1‑150 (MQ=255)
                                             gaagctgATCCCGGCCAGGGTCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCg                                                                                                    >  2:804251/4‑150 (MQ=255)
                                                 ctgACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGtcatcatcatca                                                                                                                                                     >  1:349997/1‑97 (MQ=255)
                                                 ctgACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGgcatcatcatca                                                                                                                                                     <  2:349997/97‑1 (MQ=255)
                                                    aCCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAAc                                                                                             >  2:517050/1‑150 (MQ=255)
                                                      ccTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACaa                                                                                           >  1:1188148/1‑150 (MQ=255)
                                                      ccTGGCCGGGATCAGCTTCGTCACCGACTACGGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACaa                                                                                           >  1:1102013/1‑150 (MQ=255)
                                                       cTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAg                                                                                          <  1:1256567/150‑1 (MQ=255)
                                                       cTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAg                                                                                          <  2:1137418/150‑1 (MQ=255)
                                                          gCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCt                                                                                       <  2:418824/150‑1 (MQ=255)
                                                              ggATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCg                                                                                   >  2:1377440/1‑150 (MQ=255)
                                                               gATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGa                                                                                  >  1:849794/1‑150 (MQ=255)
                                                                     cTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCt                                                                            >  2:389425/1‑150 (MQ=255)
                                                                      ttCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCa                                                                                                             <  1:237730/116‑1 (MQ=255)
                                                                      ttCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCa                                                                                                             >  2:237730/1‑116 (MQ=255)
                                                                      ttCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTc                                                                           >  2:1158851/1‑150 (MQ=255)
                                                                          tCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAATCAGACGACGCCGACCTGACCAAGGGCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGATCAAAg                                                                       >  1:390045/1‑150 (MQ=255)
                                                                             ccGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGc                                                                                     >  1:197747/1‑133 (MQ=255)
                                                                             ccGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGc                                                                                     <  2:197747/133‑1 (MQ=255)
                                                                             ccGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATAACCGTTGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCaa                                                                    <  1:804251/150‑1 (MQ=255)
                                                                              cGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAAt                                                                   >  1:1303880/1‑150 (MQ=255)
                                                                               gACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATg                                                                  <  2:510397/150‑1 (MQ=255)
                                                                                    ctGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTgcg                                                             >  1:199576/1‑150 (MQ=255)
                                                                                     tGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTgcgc                                                            >  2:465301/1‑150 (MQ=255)
                                                                                       gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAGGCCCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGAc                                                                                 <  1:449042/127‑1 (MQ=255)
                                                                                       gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGAc                                                                                 >  2:449042/1‑127 (MQ=255)
                                                                                       gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCAc                                                          >  2:1042652/1‑150 (MQ=255)
                                                                                       gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCAc                                                          >  1:92019/1‑150 (MQ=255)
                                                                                         tGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTa                                                        >  1:698445/1‑150 (MQ=255)
                                                                                                     cAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGgtgt                                            <  1:389425/150‑1 (MQ=255)
                                                                                                                gacgacGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAg                                 >  2:684924/1‑150 (MQ=255)
                                                                                                                      gCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                            <  2:307219/69‑1 (MQ=255)
                                                                                                                      gCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                            >  1:307219/1‑69 (MQ=255)
                                                                                                                      gCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATc                           >  2:812314/1‑150 (MQ=255)
                                                                                                                        cGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCgg                         <  1:517050/150‑1 (MQ=255)
                                                                                                                        cGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCgg                         <  2:897567/150‑1 (MQ=255)
                                                                                                                           gcTGACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCggg                        <  2:1243760/147‑1 (MQ=39)
                                                                                                                           gcTGACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCggg                        >  1:1243760/2‑148 (MQ=39)
                                                                                                                            cTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGAc                                                                                 <  2:404054/90‑1 (MQ=255)
                                                                                                                            cTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGAc                                                                                 >  1:404054/1‑90 (MQ=255)
                                                                                                                             tGACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAg                    >  2:1070906/1‑150 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTa                                                        >  1:1188844/1‑110 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTa                                                        <  2:1188844/110‑1 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCaa                     <  1:1222722/145‑1 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCaa                     >  2:1222722/1‑145 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGcccc                >  1:786489/1‑150 (MQ=255)
                                                                                                                                 cAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGcccc                >  1:840588/1‑150 (MQ=255)
                                                                                                                                   aGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCcgc              <  2:849794/150‑1 (MQ=255)
                                                                                                                                      tcatcatcatcaTGAACATCGTTAGTGTCACCAGCCTTCTGCGCTTCGATCAGGAGGACCGCGCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGcc           >  2:15324/1‑150 (MQ=255)
                                                                                                                                         tcatcatcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTgc        <  2:863006/150‑1 (MQ=255)
                                                                                                                                           atcatcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGTCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTgcgc      >  1:178580/1‑150 (MQ=255)
                                                                                                                                              atcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCACGCCCCGCGCCTGCGCGaa   <  2:1128665/150‑1 (MQ=255)
                                                                                                                                              atcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGaa   <  2:840588/150‑1 (MQ=255)
                                                                                                                                               tcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAg  <  2:1122274/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCTGGTGCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAG  >  NC_002947/3030994‑3031286

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