Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,982,389 G→T 5.2% L70M (CTG→ATG) ‡ cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,982,3890GT 5.2% 714.7 / 18.6 212L70M (CTG→ATG) ‡cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (109/92);  new base T (5/6);  total (114/98)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 218 total aligned reads displayed.

CGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTCTGGCCG  >  NC_002947/3982240‑3982537
                                                                                                                                                     |                                                                                                                                                    
     ggTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGa                                                                                                                                                 >  2:536/1‑150 (MQ=255)
       tCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGAcc                                                                                                                                               >  1:2078197/1‑150 (MQ=255)
        cATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACCTGTATGTCGAcca                                                                                                                                              <  1:1026402/150‑2 (MQ=255)
         aTCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGa                                                                                                                                             >  2:1594745/1‑150 (MQ=255)
         aTCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGa                                                                                                                                             >  1:2441017/1‑150 (MQ=255)
          tCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCCGGTCGACCTGTATGTCGAccagg                                                                                                                                            <  1:1863478/150‑4 (MQ=255)
            gCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGc                                                                                                                                          <  2:1706880/150‑1 (MQ=255)
              ccGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCAt                                                                                                                                        >  1:5050695/1‑150 (MQ=255)
               cGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATg                                                                                                                                       <  2:2186152/150‑1 (MQ=255)
               cGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATg                                                                                                                                       >  2:1785956/1‑150 (MQ=255)
                    ccAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCt                                                                                                                                  <  1:4491419/150‑1 (MQ=255)
                    ccAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCt                                                                                                                                  <  1:4328264/150‑1 (MQ=255)
                    ccAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCt                                                                                                                                  <  2:3157983/150‑1 (MQ=255)
                      aGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGc                                                                                                                                >  2:2922202/1‑150 (MQ=255)
                          tCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCCTGGTTCTGCACCt                                                                                                                            >  2:2755968/1‑150 (MQ=255)
                           cAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTg                                                                                                                           >  1:2260800/1‑150 (MQ=255)
                           cAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTg                                                                                                                           >  2:3214267/1‑150 (MQ=255)
                             gCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCa                                                                                                                         >  2:1258273/1‑150 (MQ=255)
                             gCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCa                                                                                                                         >  2:3840150/1‑150 (MQ=255)
                                    gCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGtt                                                                                                                  <  2:4966402/150‑1 (MQ=255)
                                    gCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGtt                                                                                                                  >  2:429746/1‑150 (MQ=255)
                                      cACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTgg                                                                                                                >  2:2436775/1‑150 (MQ=255)
                                        cGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGcc                                                                                                              >  2:1163119/1‑150 (MQ=255)
                                        cGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGcc                                                                                                              <  1:673147/150‑1 (MQ=255)
                                            cggcCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGgc                                                                                                          <  2:1908776/150‑1 (MQ=255)
                                                cAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAg                                                                                                      <  2:2783855/150‑1 (MQ=255)
                                                cAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAg                                                                                                      <  2:2066849/150‑1 (MQ=255)
                                                    tCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTg                                                                                                  <  1:1423081/150‑1 (MQ=255)
                                                    tCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTg                                                                                                  <  2:3646761/150‑1 (MQ=255)
                                                     cGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTgg                                                                                                 >  2:1753117/1‑150 (MQ=255)
                                                      gACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTggg                                                                                                <  2:4149442/150‑1 (MQ=255)
                                                      gACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTggg                                                                                                >  2:1654988/1‑150 (MQ=255)
                                                        ccAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGaa                                                                                              >  1:3228700/1‑150 (MQ=255)
                                                         cAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGaaa                                                                                             <  1:3359973/150‑1 (MQ=255)
                                                         cAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGaaa                                                                                             <  1:1887860/150‑1 (MQ=255)
                                                            ttgttCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTc                                                                                          >  2:4596169/1‑150 (MQ=255)
                                                              gttCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGACAGTCtt                                                                                        >  2:827270/1‑150 (MQ=255)
                                                               ttCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTc                                                                                       <  1:4895053/150‑1 (MQ=255)
                                                               ttCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTc                                                                                       <  1:727466/150‑1 (MQ=255)
                                                                tCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCg                                                                                      <  2:3583645/150‑1 (MQ=255)
                                                                tCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCg                                                                                      <  2:749261/150‑1 (MQ=255)
                                                                  gACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCggg                                                                                    <  2:1421910/150‑1 (MQ=255)
                                                                  gACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCggg                                                                                    >  1:1331985/1‑150 (MQ=255)
                                                                    cGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCa                                                                                  >  1:959235/1‑150 (MQ=255)
                                                                     gCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCaa                                                                                 >  2:280915/1‑150 (MQ=255)
                                                                     gCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCaa                                                                                 >  2:2134342/1‑150 (MQ=255)
                                                                      ccGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAAt                                                                                <  1:3214267/150‑1 (MQ=255)
                                                                       cGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAAtt                                                                               <  2:4846508/150‑1 (MQ=255)
                                                                           gcgcCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGctc                                                                           >  1:1566435/1‑150 (MQ=255)
                                                                             gcCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTg                                                                         >  1:2599846/1‑150 (MQ=255)
                                                                               cAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCg                                                                       >  2:1188238/1‑150 (MQ=255)
                                                                                 gTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGgt                                                                     <  1:3980954/150‑1 (MQ=255)
                                                                                    gCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTc                                                                  <  1:3934484/150‑1 (MQ=255)
                                                                                    gCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTc                                                                  >  2:4126035/1‑150 (MQ=255)
                                                                                    gCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGCCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTc                                                                  >  2:1959131/1‑150 (MQ=255)
                                                                                       aaCGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCCTCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGc                                                               >  2:1377136/1‑150 (MQ=255)
                                                                                       aaCGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGc                                                               <  1:1374161/150‑1 (MQ=255)
                                                                                         cGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGa                                                             <  1:2521721/150‑1 (MQ=255)
                                                                                         cGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGa                                                             <  1:2886334/150‑1 (MQ=255)
                                                                                         cGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGa                                                             <  1:422852/150‑1 (MQ=255)
                                                                                         cGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCAGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGa                                                             <  2:1322945/150‑1 (MQ=255)
                                                                                           ccACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGc                                                           >  1:2952976/1‑150 (MQ=255)
                                                                                              ccGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAg                                                        <  1:3137888/150‑1 (MQ=255)
                                                                                                  gCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCa                                                    >  1:2709396/1‑150 (MQ=255)
                                                                                                  gCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCa                                                    >  1:2746227/1‑150 (MQ=255)
                                                                                                      cGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGgtg                                                <  2:4850715/150‑1 (MQ=255)
                                                                                                      cGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGgtg                                                <  1:2699096/150‑1 (MQ=255)
                                                                                                        acacCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGtgtg                                              >  2:3835040/1‑150 (MQ=255)
                                                                                                        acacCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGtgtg                                              <  1:4551298/150‑1 (MQ=255)
                                                                                                            cAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATc                                          >  1:10901/1‑150 (MQ=255)
                                                                                                            cAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATc                                          >  1:3592938/1‑150 (MQ=255)
                                                                                                              ggATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATcgc                                        >  2:4709230/1‑150 (MQ=255)
                                                                                                                    ccTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGc                                  >  1:3828649/1‑150 (MQ=255)
                                                                                                                    ccTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGc                                  >  1:868144/1‑150 (MQ=255)
                                                                                                                     cTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCa                                 >  2:4908189/1‑150 (MQ=255)
                                                                                                                       tggcgtggcgCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCaga                               <  1:111472/150‑1 (MQ=255)
                                                                                                                        ggcgtggcgCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCagag                              >  2:3027297/1‑150 (MQ=255)
                                                                                                                        ggcgtggcgCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCagag                              <  2:2367327/150‑1 (MQ=255)
                                                                                                                            tggcgCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATc                          >  1:61600/1‑150 (MQ=255)
                                                                                                                              gcgcAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCAc                        <  2:3326883/150‑1 (MQ=255)
                                                                                                                              gcgcAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCAc                        <  1:4676433/150‑1 (MQ=255)
                                                                                                                               cgcAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCAcc                       <  1:1759535/150‑1 (MQ=255)
                                                                                                                                gcagcaCCTGGTCGACATACAGGTCGACCGAGGCCTGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCa                      >  2:196959/1‑150 (MQ=255)
                                                                                                                                      ctcggTCGACCTGTATGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTcgg                >  2:4192269/4‑150 (MQ=255)
                                                                                                                                       tcggTCGACCTGTATGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGtt                                                                                                                  <  2:3887534/49‑1 (MQ=21)
                                                                                                                                       tcggTCGACCTGTATGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGtt                                                                                                                  >  1:3887534/3‑51 (MQ=21)
                                                                                                                                       tcggTCGACCTGTATGTCGACCAAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTc                                                                                       >  1:4784310/3‑78 (MQ=37)
                                                                                                                                       tcggTCGACCTGTATGTCGACCAAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTc                                                                                       <  2:4784310/76‑1 (MQ=37)
                                                                                                                                       cTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGaaa                                                                                             >  1:2335349/1‑72 (MQ=255)
                                                                                                                                       cTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGaaa                                                                                             <  2:2335349/72‑1 (MQ=255)
                                                                                                                                        tGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTcggcg              >  2:2730737/1‑150 (MQ=255)
                                                                                                                                          gTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCggg            >  2:1478349/1‑150 (MQ=255)
                                                                                                                                           tCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGgt           <  2:4039029/150‑1 (MQ=255)
                                                                                                                                           tCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGgt           >  2:4974957/1‑150 (MQ=255)
                                                                                                                                              aCATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTc        >  2:4417222/1‑150 (MQ=255)
                                                                                                                                              aCATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTc        >  2:1147128/1‑150 (MQ=255)
                                                                                                                                              aCATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTc        >  2:741156/1‑150 (MQ=255)
                                                                                                                                                  aCAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCAcc                       <  1:1739019/131‑1 (MQ=255)
                                                                                                                                                    aGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTcggcgg             <  1:3046473/139‑1 (MQ=255)
                                                                                                                                                    aGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTCTGGCCg  >  2:592117/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCGAGGCATGGTTCTGCACCTGCATCGGGTTGGCCAGGCGCAGGCTGGGAAAGTCTTCGGGCAATTGCTCTGCGGTGTCGGCGAGCAAGGCCAGGTGTGAATCGCCGCTGCAGAGAATCACCAGCTCGGCGGGTGTCTGGCCG  >  NC_002947/3982240‑3982537

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