Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,629,456 A→C 11.9% *322Y (TAA→TAC mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,629,4560AC11.9% 182.8 / 17.1 67*322Y (TAA→TACmucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (5/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

CCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTC  >  NC_002947/1629309‑1629602
                                                                                                                                                   |                                                                                                                                                  
ccAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAAt                                                                                                                                                   >  3:157003/1‑149 (MQ=255)
ccAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAAt                                                                                                                                                   >  8:284362/1‑149 (MQ=255)
   gCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGgc                                                                                                                                                <  7:284362/149‑1 (MQ=255)
       ggCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGCACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGgcgcgc                                                                                                                                            >  4:205715/1‑149 (MQ=255)
         cccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCtt                                                                                                                                          >  3:238615/1‑149 (MQ=255)
         cccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcct                                                                                                                                          >  3:357149/1‑147 (MQ=255)
          cctcccTCGGCAGGCTCGCGTCGCCTGAATACACCTAAGGGAAAAGTGATGGTGACCGTGGGTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCATTACTGgcgcgcctg                                                                                                                                         <  8:351249/145‑4 (MQ=255)
          ccGCCCTCGGCAGTCTCGCGTCGCCTGACTCCACCTCCGGGCAAAGTGAGGGTGACCGTGGGTGGCGAACTTCCGCTGGGTACGGCAGCACGCGTCGCCCTGGCGATGCGCCCCCAGGAGGCCCAGGCGCGCCATTACTGgcgcgcctg                                                                                                                                         <  8:346294/149‑4 (MQ=255)
          ccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCttt                                                                                                                                         >  7:283938/1‑149 (MQ=255)
           cGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                                        <  4:357149/149‑5 (MQ=255)
           cGACCTCCGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCttgg                                                                                                                                        <  3:213310/149‑3 (MQ=255)
              ccTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCg                                                                                                                                     <  1:287039/149‑1 (MQ=255)
                       tctcGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGATGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGc                                                                                                                             >  4:39839/1‑148 (MQ=255)
                         tcGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGcc                                                                                                                            >  6:312640/1‑147 (MQ=255)
                            cgTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGc                                                                                                                       <  3:74204/149‑1 (MQ=255)
                                 ccTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                                  <  2:117939/149‑1 (MQ=255)
                                  cTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTgg                                                                                                                 <  3:123723/149‑1 (MQ=255)
                                  cTGAATACACCCAAGGGCAAAGTGATGATGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTgg                                                                                                                 <  3:39839/149‑1 (MQ=255)
                                    gAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTggg                                                                                                                >  3:123317/1‑148 (MQ=255)
                                    gAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGc                                                                                                               >  5:102630/1‑149 (MQ=255)
                                    gAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGAAGAACGCGTCCCCCTGTCGATGCCCCCCCAGGATGCCCAGGCGCGGCAGTAATGGCGCGCGTTTGCGTGCCGCGCCAAGGCAGGTggg                                                                                                                >  1:106335/1‑148 (MQ=255)
                                     aaTACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCa                                                                                                              >  1:128528/1‑149 (MQ=255)
                                       tACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATa                                                                                                            >  1:318321/1‑149 (MQ=255)
                                          acCAAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGATGGGTAAGGCAGAACGCGTCGCCCTGTCGATGCGCCCCAAGGATGCCCAGGAGCGCCATTAATGGCGCGATTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGc                                                                                                         >  6:381559/1‑149 (MQ=255)
                                            ccAAGGGCAAATTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                       >  1:127251/1‑149 (MQ=255)
                                             cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                       >  5:241999/1‑148 (MQ=255)
                                             cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                       >  3:131479/1‑148 (MQ=255)
                                               aGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaa                                                                                                    <  3:205715/149‑1 (MQ=255)
                                               aGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaa                                                                                                    <  6:241999/149‑1 (MQ=255)
                                                gggCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaa                                                                                                   <  8:283938/149‑1 (MQ=255)
                                                 ggCTAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGTGTACGGCAGAACGCGTCGACCTGTCGATGCTGGCCCAGGATGCGCTGGCGCGCCAGTAATGGCGCTCTTTTGGGTGGAGAGCCATGGGAGGTGGGGATAGGGGtgaaa                                                                                                   <  2:127251/148‑1 (MQ=255)
                                                   cAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatga                                                                                                <  5:312640/149‑1 (MQ=255)
                                                   cAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatg                                                                                                 >  2:243709/1‑148 (MQ=255)
                                                    aaaGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatgaa                                                                                               <  4:369129/149‑1 (MQ=255)
                                                      aGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTgaaatgaaat                                                                                             <  1:243709/149‑1 (MQ=255)
                                                       gTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAAtt                                                                                            <  4:157003/149‑1 (MQ=255)
                                                                            gAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATc                                                                       >  4:289527/1‑149 (MQ=255)
                                                                            gAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATc                                                                       >  8:200412/1‑149 (MQ=255)
                                                                                ttCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaa                                                                                                   <  8:9542/117‑1 (MQ=255)
                                                                                ttCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaa                                                                                                   >  7:9542/1‑117 (MQ=255)
                                                                                            cGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTc                                                        <  2:128528/148‑1 (MQ=255)
                                                                                              gCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCtat                                                     <  2:318321/149‑1 (MQ=255)
                                                                                                   aCGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTg                                                                                   <  3:239499/114‑1 (MQ=255)
                                                                                                   aCGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTg                                                                                   >  4:239499/1‑114 (MQ=255)
                                                                                                       gTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCAGCGGCAAATTTTCTATAGGTCAGGc                                            >  3:32755/1‑149 (MQ=255)
                                                                                                              tGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTg                                        >  7:95528/1‑146 (MQ=255)
                                                                                                                  gATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGc                                                                >  8:50289/1‑118 (MQ=255)
                                                                                                                  gATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGc                                                                <  7:50289/118‑1 (MQ=255)
                                                                                                                  gATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGcc                                 <  4:123317/149‑1 (MQ=255)
                                                                                                                    tGCGCCGCCAGGATTCCCAGGCGCGGCAGTAATGTCGCGATTTTGAGTGCCGAGCCAAGTCAGGTGTGCATACGACTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTg                               <  5:381559/149‑1 (MQ=255)
                                                                                                                    tGCGCCCCCAGGATGCCCATGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCt                                <  3:195111/148‑1 (MQ=255)
                                                                                                                       gcCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAttt                                                                             <  7:191436/100‑1 (MQ=255)
                                                                                                                       gcCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAttt                                                                             >  8:191436/1‑100 (MQ=255)
                                                                                                                           ccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCtttt                        <  3:289527/149‑1 (MQ=255)
                                                                                                                           ccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGCCAGGCTCTTCAAAGCCTGACCtata                        <  4:238615/149‑4 (MQ=255)
                                                                                                                                aTGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCg                   <  8:95528/149‑1 (MQ=255)
                                                                                                                                  gcaaaagcgcgcCATTACTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGtt                 >  6:198083/7‑149 (MQ=255)
                                                                                                                                    ccAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGc                                                                                                         >  2:297061/1‑59 (MQ=255)
                                                                                                                                    ccAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGc                                                                                                         <  1:297061/59‑1 (MQ=255)
                                                                                                                                    aaaagcgcgcCATTACTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgt               >  2:358122/5‑149 (MQ=255)
                                                                                                                                     cAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgt               >  5:305501/1‑148 (MQ=255)
                                                                                                                                      aagcgcgcCATTACTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTTGGCATAAGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTATATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTATGCGTtgtg              >  1:269153/3‑148 (MQ=255)
                                                                                                                                      aagcgcgcCATTACTGGCGCGCTTTTGCGTGCAGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg              >  1:234265/3‑148 (MQ=255)
                                                                                                                                      aGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg              >  8:132490/1‑148 (MQ=255)
                                                                                                                                       ggCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg              >  2:34425/1‑147 (MQ=255)
                                                                                                                                            gcCAGTAATGGCGCGCTTTTGAGTGCCGAGCCAAGGCACGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTAAGGCTTTGAAGAGCCTGCCCTTTTCTGCGTTGTGCAGGGGt       >  3:209295/1‑149 (MQ=255)
                                                                                                                                                 tAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCTTACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTc  <  6:102630/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTC  >  NC_002947/1629309‑1629602

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