Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,916,385 G→C 6.1% intergenic (+87/‑28) PP_1715 → / → PP_1716 hypothetical protein/Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,916,3850GC 6.1% 211.5 / 3.7 66intergenic (+87/‑28)PP_1715/PP_1716hypothetical protein/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base C (2/2);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

CGCGTGGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAG  >  NC_002947/1916239‑1916532
                                                                                                                                                  |                                                                                                                                                   
cgcgTGGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCt                                                                                                                                                   <  8:190923/149‑1 (MQ=255)
    tGGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTaa                                                                                                                                               <  6:185928/149‑1 (MQ=255)
    tGGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTaa                                                                                                                                               <  4:242250/149‑1 (MQ=255)
       tCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCTGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACAGCATGTAATcc                                                                                                                                            <  3:109978/149‑1 (MQ=255)
         gAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCtt                                                                                                                                          >  8:81970/1‑149 (MQ=255)
         gAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCtt                                                                                                                                          >  3:229202/1‑149 (MQ=255)
          aGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCttt                                                                                                                                         <  1:212422/149‑1 (MQ=255)
          aGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACAGCATGTAATccgta                                                                                                                                         <  7:4040/149‑4 (MQ=255)
            cgcCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCttttt                                                                                                                                       <  7:81970/149‑1 (MQ=255)
                tGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAg                                                                                                                                   >  8:32790/1‑149 (MQ=255)
                  ccGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGTTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAggg                                                                                                                                 <  8:47847/149‑1 (MQ=255)
                    ggCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGt                                                                                                                               <  8:35997/149‑1 (MQ=255)
                     gCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGtt                                                                                                                              <  4:102476/149‑1 (MQ=255)
                      cTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGttt                                                                                                                             <  4:71368/149‑1 (MQ=255)
                      cTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGGTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGttt                                                                                                                             <  7:168334/149‑1 (MQ=255)
                         ccccGGAAGCGCTGCGTCTTGATCTAAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTTCACATAACTCGGAGTTTGTAATATAATTTTCATAATTGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTg                                                                                                                            <  3:161713/147‑1 (MQ=255)
                          cccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTaa                                                                                                                         >  7:37492/1‑149 (MQ=255)
                             ggAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATgg                                                                                                                      <  1:250351/149‑1 (MQ=255)
                                  cgcTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTAtt                                                                                                                 <  8:37492/149‑1 (MQ=255)
                                   gcTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTAtt                                                                                                                 >  4:226302/1‑148 (MQ=255)
                                     tGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTc                                                                                                                <  3:226302/147‑1 (MQ=255)
                                       cGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGt                                                                                                            <  8:263875/149‑1 (MQ=255)
                                          cTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTaaa                                                                                                                                    <  4:233841/122‑1 (MQ=255)
                                          cTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTaaa                                                                                                                                    >  3:233841/1‑122 (MQ=255)
                                           tGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGtcat                                                                                                         >  5:174744/1‑148 (MQ=255)
                                             gAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGtcatc                                                                                                        <  3:52150/147‑1 (MQ=255)
                                              aGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGc                                                                                                     >  7:271863/1‑149 (MQ=255)
                                                   aacaaGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGa                                                                                                <  6:174744/149‑1 (MQ=255)
                                                      aaGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTa                                                                                                                   <  2:91258/127‑1 (MQ=255)
                                                      aaGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACgg                                                                                             >  7:30643/1‑149 (MQ=255)
                                                      aaGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATAGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTa                                                                                                                   >  1:91258/1‑127 (MQ=255)
                                                               aaCCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGa                                                                                    <  8:271863/149‑1 (MQ=255)
                                                                 tctgCCAGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAAAATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGt                                                                                  <  6:291865/146‑1 (MQ=255)
                                                                         aaTCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCg                                                                            <  4:287129/147‑1 (MQ=255)
                                                                          aTCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTa                                                                          <  4:229202/148‑1 (MQ=255)
                                                                             tcACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATCACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTc                                                                      >  8:52771/1‑149 (MQ=255)
                                                                              cACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCaa                                                                                         <  3:211836/129‑1 (MQ=255)
                                                                              cACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCaa                                                                                         >  4:211836/1‑129 (MQ=255)
                                                                              cACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTcc                                                                     >  6:278331/1‑149 (MQ=255)
                                                                               aCCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGtt                                                                       <  2:173696/146‑1 (MQ=255)
                                                                                 cTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTg                                                                  <  1:287364/149‑1 (MQ=255)
                                                                                    aaCTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTgg                                                                 <  3:219416/147‑1 (MQ=255)
                                                                                         cacaGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCaa                                                           >  5:242730/1‑148 (MQ=255)
                                                                                            aGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATACTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATAc                                                        >  1:57651/1‑148 (MQ=255)
                                                                                              aCCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAg                                                     >  4:201318/1‑149 (MQ=255)
                                                                                                cTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTAAGGTCATCAGCCGGGACGGTCAACCCGAGTAAGCGGTAGTTCCGTGGGAGCAATACCAGGc                                                   >  4:174400/1‑149 (MQ=255)
                                                                                                cTCGGAGTTGGTAATATAATTGGAATAAGGGTTATGCCTAAGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGCGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGACGTAGTTCCCTGGGAGCAATACCAGGc                                                   >  6:47715/1‑149 (MQ=255)
                                                                                                   ggAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATCACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGt                                                 <  7:52771/148‑1 (MQ=255)
                                                                                                         ggTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGtcat                                                                                                         >  4:201721/1‑86 (MQ=255)
                                                                                                         ggTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGtcat                                                                                                         <  3:201721/86‑1 (MQ=255)
                                                                                                         ggTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCa                                           <  8:30643/148‑1 (MQ=255)
                                                                                                          gTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAgc                                         >  5:163554/1‑149 (MQ=255)
                                                                                                                taATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTgc                                    <  5:278331/148‑1 (MQ=255)
                                                                                                                  aTTGGCATATGGGGTATTGCTACTGATTACATGCTGTAATCCTTTTTAAATGGGTTTGTAATGGGTTTTCTGGGCATCAGCCTGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGTGAGCAATACCAGGCGTTGCTCAGCGCTgcagg                                 <  5:47715/149‑1 (MQ=255)
                                                                                                                   ttGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTgcaggg                                >  6:312107/1‑149 (MQ=255)
                                                                                                                        aTAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTgcagggcagg                            >  3:52156/1‑148 (MQ=255)
                                                                                                                           aGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCAc                         <  3:201318/148‑1 (MQ=255)
                                                                                                                               ttATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAACAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGcc                    >  1:101983/1‑149 (MQ=255)
                                                                                                                                     aaaaggATTACAGCATGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCg              >  3:192291/5‑149 (MQ=255)
                                                                                                                                      aaaggATTACAGCATGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCgg             >  5:301740/4‑149 (MQ=255)
                                                                                                                                           aTTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTc         >  2:273702/1‑148 (MQ=255)
                                                                                                                                            ttACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCGGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCa       <  5:312107/149‑1 (MQ=255)
                                                                                                                                              aCATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAg      <  2:57651/148‑1 (MQ=255)
                                                                                                                                               cATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGt     >  6:146081/1‑148 (MQ=255)
                                                                                                                                               cATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGt     >  2:193605/1‑148 (MQ=255)
                                                                                                                                                 tGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAg  <  3:174400/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCGTGGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAG  >  NC_002947/1916239‑1916532

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