Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,540,955 T→C 5.9% L23P (CTG→CCG)  PP_3128 → exopolysaccharide biosynthesis/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,540,9550TC 5.9% 201.5 / 2.3 70L23P (CTG→CCG) PP_3128exopolysaccharide biosynthesis/transport protein
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base C (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCTGGTGCACGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACGGGCCGATC  >  NC_002947/3540810‑3541100
                                                                                                                                                  |                                                                                                                                                 
gCTGGTGCACGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGt                                                                                                                                                 >  7:303999/1‑148 (MQ=255)
  tGGTGCACGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGtctc                                                                                                                                              <  8:303999/149‑1 (MQ=255)
      gCACGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGc                                                                                                                                          >  4:28510/1‑149 (MQ=255)
         cGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGcccc                                                                                                                                       >  4:47368/1‑149 (MQ=255)
          gtgtTCTACATGCCACTGGATGTGCTGTTGAAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGccccc                                                                                                                                      >  3:11850/1‑149 (MQ=255)
            gtTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCgg                                                                                                                                    >  2:246167/1‑149 (MQ=255)
              tCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTa                                                                                                                                  >  3:289260/1‑149 (MQ=255)
               cTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTa                                                                                                                                  <  3:47368/148‑1 (MQ=255)
                         cTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCAAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGcc                                                                                                                       >  3:421572/1‑149 (MQ=255)
                          tGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCAGTTGGCAACAATATTCATCAGGTTACACCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGcc                                                                                                                       >  7:391970/1‑148 (MQ=255)
                             aTGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAg                                                                                                                   <  7:320259/149‑1 (MQ=255)
                               gtgCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTg                                                                                                                 >  4:210484/1‑149 (MQ=255)
                                tgCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGc                                                                                                                <  7:360869/149‑1 (MQ=255)
                                   tgttgtAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             >  6:8765/1‑149 (MQ=255)
                                       gtAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGt                                                                                                           >  5:127963/1‑147 (MQ=255)
                                       gtAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTc                                                                                                          <  4:183488/148‑1 (MQ=255)
                                        tAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGt                                                                                                           >  1:20555/1‑146 (MQ=255)
                                        tAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGt                                                                                                           <  2:20555/146‑1 (MQ=255)
                                            aaGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTc                                                                                                     >  1:222290/1‑148 (MQ=255)
                                                  tGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACg                                                                                              <  1:251640/149‑1 (MQ=255)
                                                    tctcGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGt                                                                                             >  7:27326/1‑148 (MQ=255)
                                                      tcGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGaa                                                                                          >  1:395757/1‑149 (MQ=255)
                                                      tcGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGaa                                                                                          >  3:362423/1‑149 (MQ=255)
                                                        ggTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGaa                                                                                          <  4:11850/147‑1 (MQ=255)
                                                        ggTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAaca                                                                                        >  8:24276/1‑149 (MQ=255)
                                                         gTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAacac                                                                                       >  5:124917/1‑149 (MQ=255)
                                                              tataGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             <  8:256042/122‑1 (MQ=255)
                                                              tataGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             >  7:256042/1‑122 (MQ=255)
                                                              tataGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAg                                                                                   <  2:41397/148‑1 (MQ=255)
                                                               ataGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAgg                                                                                  <  1:246167/148‑1 (MQ=255)
                                                                taGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAggt                                                                                 >  1:71419/1‑149 (MQ=255)
                                                                   ggAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTa                                                                             >  6:407784/1‑149 (MQ=255)
                                                                   ggAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTa                                                                             >  3:162873/1‑149 (MQ=255)
                                                                         cAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             <  6:416696/111‑1 (MQ=255)
                                                                         cAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             <  6:237888/111‑1 (MQ=255)
                                                                         cAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             >  5:416696/1‑111 (MQ=255)
                                                                         cAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCgg                                                                                                             >  5:237888/1‑111 (MQ=255)
                                                                          aGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGc                                                                      >  6:181703/1‑149 (MQ=255)
                                                                              aTGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTg                                                                   <  7:24276/148‑1 (MQ=255)
                                                                              aTGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTg                                                                   <  5:341438/148‑1 (MQ=255)
                                                                                 gACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGa                                                               <  3:28510/149‑1 (MQ=255)
                                                                                 gACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCATGCAGa                                                               >  4:393173/1‑149 (MQ=255)
                                                                                  aCAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACCGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCaa                                                                                    <  2:429165/127‑1 (MQ=255)
                                                                                  aCAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACCGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCaa                                                                                    >  1:429165/1‑127 (MQ=255)
                                                                                  aCAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACCGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCaa                                                                                    <  2:435696/127‑1 (MQ=255)
                                                                                  aCAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACCGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCaa                                                                                    >  1:435696/1‑127 (MQ=255)
                                                                                      gATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGaccacc                                                          <  6:124917/149‑1 (MQ=255)
                                                                                       aTCAAGCCGGCAGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGaccacc                                                          <  8:391970/148‑1 (MQ=255)
                                                                                        tCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTg                                                        >  5:67427/1‑149 (MQ=255)
                                                                                              cGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCg                                                   >  5:357371/1‑148 (MQ=255)
                                                                                                          aacaaTATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGt                                       >  8:398252/1‑149 (MQ=255)
                                                                                                           acaaTATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTg                                      <  2:395757/149‑1 (MQ=255)
                                                                                                              atatTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCg                                   <  3:210484/149‑1 (MQ=255)
                                                                                                              atatTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCg                                   >  5:174020/1‑149 (MQ=255)
                                                                                                                 ttCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTac                                <  2:222290/149‑1 (MQ=255)
                                                                                                                 ttCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTa                                 <  4:162873/148‑1 (MQ=255)
                                                                                                                    atcaGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTAc                                                                                               >  1:270183/1‑83 (MQ=255)
                                                                                                                    atcaGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTAc                                                                                               <  2:270183/83‑1 (MQ=255)
                                                                                                                        ggTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGa                         <  7:398252/149‑1 (MQ=255)
                                                                                                                          tgaCCGCAGGTACTCCCCACGCACAGTCTCCCGCCCCCGGTACTGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAAc                       <  8:68387/147‑1 (MQ=255)
                                                                                                                           tACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACt                      >  6:199958/1‑149 (MQ=255)
                                                                                                                              cccAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTcc                    >  1:205715/1‑148 (MQ=255)
                                                                                                                              cccAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTcc                    >  1:435256/1‑148 (MQ=255)
                                                                                                                                cAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCaa                  >  2:258432/1‑148 (MQ=255)
                                                                                                                                      cTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAAc           <  6:127963/149‑1 (MQ=255)
                                                                                                                                       tCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACg          <  4:421572/149‑1 (MQ=255)
                                                                                                                                       tCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACg          >  3:50742/1‑149 (MQ=255)
                                                                                                                                         cccAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACggg        >  4:107407/1‑149 (MQ=255)
                                                                                                                                         cccAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACgg         <  6:357371/148‑1 (MQ=255)
                                                                                                                                               ctcTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACGGGCCGATc  <  4:289260/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGGTGCACGTGTTCTACATGCCACTGGATGTGCTGTTGTAAAAAGTCATGTCTCGGTTTGTATAGGGAAAGCAGTGATGGACAGGATCAAGCCGGCTGTTGGCAACAATATTCATCAGGTTACCCCAGGTACTCCCCAGGCACTGTCTCCCGCCCCCGGTACGGCGCTGGCCGAAGTGCCGGGTCAGTTCGATTACGTGAACACCAAGGTGGTACCGCTGCGTGCAGACCACCTGGAGCGTAACCGCATCGTGGCGTACAACAAGAACTCCAACATGAACGGGCCGATC  >  NC_002947/3540810‑3541100

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: