Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,986,592 A→T 6.0% intergenic (+81/‑84) flgA → / → flgM flagella basal body P‑ring formation protein/flagellin biosynthesis negative regulator FlgM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,986,5920AT 6.0% 253.8 / 2.9 81intergenic (+81/‑84)flgA/flgMflagella basal body P‑ring formation protein/flagellin biosynthesis negative regulator FlgM
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base T (3/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCTCAATTCCAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGG  >  NC_002947/4986445‑4986734
                                                                                                                                                   |                                                                                                                                              
ccTCAATTCCAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAg                                                                                                                                               >  8:185071/1‑149 (MQ=255)
ccTCAATTCCAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAg                                                                                                                                               >  2:11665/1‑149 (MQ=255)
       tCCAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCaa                                                                                                                                        >  5:73497/1‑149 (MQ=255)
         cAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAAt                                                                                                                                       <  5:125829/148‑1 (MQ=255)
              gcgTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCAttt                                                                                                                                 <  6:73497/149‑1 (MQ=255)
              gcgTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCAttt                                                                                                                                 <  4:311587/149‑1 (MQ=255)
              gcgTGGGCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGCACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCAttt                                                                                                                                 <  2:93561/149‑1 (MQ=255)
               cgTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTtg                                                                                                                                >  3:137915/1‑149 (MQ=255)
                gTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTtgt                                                                                                                               <  2:346687/149‑1 (MQ=255)
                  ggTCCAGGCCAGGGTGCCCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGc                                                                                                                             <  8:292919/149‑1 (MQ=255)
                  ggTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGc                                                                                                                             <  7:86887/149‑1 (MQ=255)
                       aGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAATACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTa                                                                                                                          >  2:196421/1‑147 (MQ=255)
                          ccAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACTGGCGGCAATGCATTTGTGCCTAAAGtt                                                                                                                     <  1:11665/149‑1 (MQ=255)
                            aGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTta                                                                                                                   <  3:94477/149‑1 (MQ=255)
                              ggTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTtata                                                                                                                 >  5:279128/1‑149 (MQ=255)
                                 gACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCg                                                                                                              <  5:11915/149‑1 (MQ=255)
                                  aCCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCgg                                                                                                             <  4:189490/149‑1 (MQ=255)
                                   ccGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCggg                                                                                                            >  1:201816/1‑149 (MQ=255)
                                              cAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGaaaa                                                                                                 <  2:201816/149‑1 (MQ=255)
                                              cAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGaaaa                                                                                                 >  7:266293/1‑149 (MQ=255)
                                                      ggTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGc                                                                                                                             <  5:95325/113‑1 (MQ=255)
                                                      ggTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGc                                                                                                                             >  6:95325/1‑113 (MQ=255)
                                                      ggTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCaa                                                                                         <  2:7345/149‑1 (MQ=255)
                                                         cAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCg                                                                                      >  4:43345/1‑149 (MQ=255)
                                                               gTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCaaa                                                                                >  7:28855/1‑149 (MQ=255)
                                                                tAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCaaa                                                                                <  7:185071/148‑1 (MQ=255)
                                                                     tAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATAcc                                                                            >  1:269541/1‑147 (MQ=255)
                                                                     tAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATAcc                                                                            >  4:17175/1‑147 (MQ=255)
                                                                       gCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTTAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCa                                                                          <  1:129738/147‑1 (MQ=255)
                                                                        ccaccaGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCagag                                                                       <  3:43345/149‑1 (MQ=255)
                                                                         caccaGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCAGTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAgg                                                                      >  6:58798/1‑149 (MQ=255)
                                                                          accaGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCATTGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGt                                                                     >  4:118946/1‑149 (MQ=255)
                                                                          accaGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGCCTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTAGGCCGAAAACAAGCCAACCGTCCAAATACCCAGAGGt                                                                     >  1:158891/1‑149 (MQ=255)
                                                                            caGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGttt                                                                   >  6:1628/1‑149 (MQ=255)
                                                                             aGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTc                                                                  <  2:29891/149‑1 (MQ=255)
                                                                             aGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTc                                                                  <  6:187141/149‑1 (MQ=255)
                                                                                ttGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGa                                                               <  8:266293/149‑1 (MQ=255)
                                                                                  gcgcTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGAt                                                              <  2:269541/148‑1 (MQ=255)
                                                                                   cgcTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCa                                                            >  1:296207/1‑149 (MQ=255)
                                                                                     cTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATg                                                          <  7:104171/149‑1 (MQ=255)
                                                                                      tGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCaa                                                                                                                                        <  4:342840/70‑1 (MQ=255)
                                                                                      tGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCaa                                                                                                                                        >  3:342840/1‑70 (MQ=255)
                                                                                           gTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATc                                                    <  4:333345/149‑1 (MQ=255)
                                                                                                  aCCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACtt                                               >  4:232809/1‑147 (MQ=255)
                                                                                                   ccGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAg                                            <  3:232809/149‑1 (MQ=255)
                                                                                                     gCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCaa                                                                                         <  5:130185/102‑1 (MQ=255)
                                                                                                     gCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCaa                                                                                         >  6:130185/1‑102 (MQ=255)
                                                                                                      cTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCg                                         >  6:197107/1‑149 (MQ=255)
                                                                                                      cTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCg                                         >  6:74500/1‑149 (MQ=255)
                                                                                                          tCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGttt                                      <  6:279128/148‑1 (MQ=255)
                                                                                                            cTAAACTGTGACACAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAAAGAATTTCTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCCTCCAAATACCCAGAGGATTCTGATCATGGTCATAGACTTCAGTCgaatga                                    >  6:161750/1‑143 (MQ=255)
                                                                                                             tAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGaa                                   >  6:217378/1‑148 (MQ=255)
                                                                                                             tAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGACGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGaa                                   >  3:91799/1‑148 (MQ=255)
                                                                                                                aCTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTTCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATaa                                <  7:3794/148‑1 (MQ=255)
                                                                                                                aCTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATAcc                                                                            >  3:355350/1‑104 (MQ=255)
                                                                                                                aCTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATAcc                                                                            <  4:355350/104‑1 (MQ=255)
                                                                                                                aCTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAAc                               >  6:68346/1‑149 (MQ=255)
                                                                                                                aCTGTGTCAGAACACGGGTTCGCGGGCTTGCCGAAAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAAc                               >  6:73753/1‑149 (MQ=255)
                                                                                                                    tgtCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAActct                            <  5:197107/148‑1 (MQ=255)
                                                                                                                    tgtCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAActct                            <  4:137915/148‑1 (MQ=255)
                                                                                                                      tCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCg                                                                                                              >  3:324917/1‑64 (MQ=255)
                                                                                                                      tCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCg                                                                                                              <  4:324917/64‑1 (MQ=255)
                                                                                                                      tCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAg                                            >  5:377240/1‑130 (MQ=255)
                                                                                                                      tCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAg                                            <  6:377240/130‑1 (MQ=255)
                                                                                                                      tCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCtccg                         >  1:125788/1‑149 (MQ=255)
                                                                                                                        aGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCtccgtc                       >  2:215781/1‑149 (MQ=255)
                                                                                                                        aGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCtccgt                        <  2:296207/148‑1 (MQ=255)
                                                                                                                         gAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCtccgtcc                      >  6:281723/1‑149 (MQ=255)
                                                                                                                         gAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCtccgtcc                      >  8:181464/1‑149 (MQ=255)
                                                                                                                             acGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTATGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCa                  >  6:33492/1‑149 (MQ=255)
                                                                                                                                ggTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACgg               <  5:1628/149‑1 (MQ=255)
                                                                                                                                    cgcgGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGggcg            <  8:28855/148‑1 (MQ=255)
                                                                                                                                    aaatgcattGCCGCCTGTCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGggcgg           >  2:64898/8‑149 (MQ=255)
                                                                                                                                     gcgGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGCCCGTGacggacgga          <  5:68346/149‑2 (MQ=255)
                                                                                                                                      atgcattGCCGCCTGTCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACggg              >  1:10712/6‑144 (MQ=255)
                                                                                                                                      atgcattGCCGCCTGTCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACggg              <  2:10712/139‑1 (MQ=255)
                                                                                                                                       gggCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGCCCGTGacggacgga          <  5:73753/147‑2 (MQ=255)
                                                                                                                                         gCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAATCAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTgc      <  5:161750/149‑1 (MQ=255)
                                                                                                                                         gCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTgc      >  5:98843/1‑149 (MQ=255)
                                                                                                                                         cattGCCGCCTGTCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTgc      >  6:217560/3‑149 (MQ=255)
                                                                                                                                              ccGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCgg  <  5:217378/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCTCAATTCCAAACGCGTGGTCAAGGCCAGGGTGACCGGCCCGGGCCAGGTCGAGGTCAGTATGTAGGCTAGCCACCAGATTGCGCTGGCAGTGAGGAACCGCTTTTCCTAAACTGTGTCAGAACACGGGTTCGCGGGCTTGCCGACAGGCGGCAATGCATTTGTGCCTAAAGTTTATATCGGGTTGGCCGAAAACAAGGCAAGCGTCCAAATACCCAGAGGTTTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGG  >  NC_002947/4986445‑4986734

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: