Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,640,394 G→T 8.8% intergenic (‑32/+26) PP_4105 ← / ← gabP‑IV hypothetical protein/gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,640,3940GT 8.8% 175.9 / 4.1 61intergenic (‑32/+26)PP_4105/gabP‑IVhypothetical protein/gamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base T (2/3);  total (29/32)
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.88e-01

ACCGTCGGCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAGCC  >  NC_002947/4640246‑4640542
                                                                                                                                                    |                                                                                                                                                    
acCGTCGGCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGgcgcg                                                                                                                                                      >  5:109141/1‑149 (MQ=255)
       gCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGcc                                                                                                                                               <  3:66641/149‑1 (MQ=255)
            ccAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCa                                                                                                                                           <  6:241548/148‑1 (MQ=255)
             cAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCAATAGCAAGTACACACAGGGGCCGAGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTACTGTGCTTGTTGGAATGGGGGCGCGCAACGCACCCaa                                                                                                                                          >  1:84224/1‑147 (MQ=255)
             cAGTGGCTTTGCAGGTGCAGGCCACTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGTTGCGCGCCCCCAtt                                                                                                                                         <  3:102784/149‑1 (MQ=255)
              aGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCAtt                                                                                                                                         >  8:345275/1‑148 (MQ=255)
               gTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTaa                                                                                                                                       <  8:15243/149‑1 (MQ=255)
                tGGGTTTGCAGGTTCTTGGCGGTGGGGTTGGAACAGCCGCATGCATTGCAAGTACCCTCAGGGGGCGCGATGGCTTGAATGGTAATTGTTTTATCAGTCATACGCTTCCTGTTCTTGTTTGAATGGGGGCGTTGCGCGCCCCCAttc                                                                                                                                        <  2:116412/147‑2 (MQ=255)
                tGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTaaa                                                                                                                                      <  6:109141/149‑1 (MQ=255)
                tGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGAATGAAGGGAAAGTGTGAGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGTTGCGCGCCCCCAttcca                                                                                                                                      >  1:116412/1‑146 (MQ=255)
                 ggCTTTGCAGGTGCAGGCCGCTGGGCTGCCAAAAGCCGCAGCAATATCAACTACCCTCAGGGCACGCGAGCGCGAGAAGGGAAACTCAGTGATCAGTCCTACTCTTCCTGTGCTTGATTGAATAGTCTCGCGCAAAGCCCCCATTaaat                                                                                                                                     >  8:17104/1‑148 (MQ=255)
                   cTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCa                                                                                                                                   <  8:150545/149‑1 (MQ=255)
                       gcaggtgcaggCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGAGATCAGACATACGATACCTGTGCATGTTGGAATGGGGGCGCACAACGCCCCCAATAAACCACaa                                                                                                                                >  8:237690/1‑148 (MQ=255)
                           gtgcaggCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATc                                                                                                                           <  2:236666/149‑1 (MQ=255)
                             gcaggCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGc                                                                                                                                                >  4:203487/1‑126 (MQ=255)
                             gcaggCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGc                                                                                                                                                <  3:203487/126‑1 (MQ=255)
                              caggCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAgg                                                                                                                        >  6:388484/1‑149 (MQ=255)
                                    gCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCt                                                                                                                  <  7:345275/149‑1 (MQ=255)
                                     cTGGGGTGGGAACAGCCGCAGCCATAGGAAGTACCCTCAGTGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCGTCCTGTGCTTGTTGGAATGGGCGCGAGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGGCTg                                                                                                                 <  3:212910/149‑1 (MQ=255)
                                      tGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTgg                                                                                                                <  5:121745/149‑1 (MQ=255)
                                            gggAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAAc                                                                                                           >  2:330420/1‑148 (MQ=255)
                                             ggAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGTTGCGCGCCCCCAttcc                                                                                                                                       <  7:310898/119‑3 (MQ=255)
                                             ggAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGTTGCGCGCCCCCAttcc                                                                                                                                       >  8:310898/1‑117 (MQ=255)
                                                aCAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTgg                                                                                                                >  6:380940/1‑139 (MQ=255)
                                                aCAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTgg                                                                                                                <  5:380940/139‑1 (MQ=255)
                                                  aGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTc                                                                                                    <  8:295602/149‑1 (MQ=255)
                                                    ccGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCg                                                                                                   >  2:253526/1‑148 (MQ=255)
                                                    ccGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGc                                                                                                  <  4:211171/149‑1 (MQ=255)
                                                    ccGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGc                                                                                                  <  1:227819/149‑1 (MQ=255)
                                                      gCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCAc                                                                                                >  4:120131/1‑149 (MQ=255)
                                                        aGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCAcgc                                                                                              >  4:96804/1‑149 (MQ=255)
                                                                cAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGt                                                                                      >  5:226359/1‑149 (MQ=255)
                                                                   gTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATGAAAGCAGACACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAg                                                                                   >  5:109119/1‑149 (MQ=255)
                                                                    tACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGc                                                                                  <  7:17104/149‑1 (MQ=255)
                                                                        cTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGaaa                                                                              <  2:108984/149‑1 (MQ=255)
                                                                           aGGGGCCGCGTGGGCTTGTAGTGAAAGTGTGTGTTGAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACTTCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAAc                                                                             <  6:199614/147‑1 (MQ=255)
                                                                           aGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAAc                                                                             <  5:266693/147‑1 (MQ=255)
                                                                           aGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAAc                                                                             <  5:388484/147‑1 (MQ=255)
                                                                           aGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCa                                                                           <  1:330420/149‑1 (MQ=255)
                                                                               gCCGCGAGGGCTTGAAGGCAAAGTGTGTGATCAGTCATACGCTTCCTCTGCTTGTTGGAATGGGGGCGCGCACCGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCTTCTTCGCACGCCGCTGGCTCAGCGAAACCAGCa                                                                        >  6:273809/1‑148 (MQ=255)
                                                                                  gcgAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAAt                                                                    >  4:337804/1‑149 (MQ=255)
                                                                                         cTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAAcgcg                                                             <  1:253526/149‑1 (MQ=255)
                                                                                         cTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAAcgcg                                                             <  3:337804/149‑1 (MQ=255)
                                                                                              aGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCAc                                                        >  4:83733/1‑149 (MQ=255)
                                                                                               gggAAAGTGTGTGATCAGCCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCAcc                                                       >  2:201739/1‑149 (MQ=255)
                                                                                                   aaGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCgg                                                   >  2:2906/1‑149 (MQ=255)
                                                                                                                   attcGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTc                                   <  7:346653/146‑1 (MQ=255)
                                                                                                                   aTACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTc                                   >  5:274479/1‑149 (MQ=255)
                                                                                                                    tACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCg                                  >  6:303647/1‑149 (MQ=255)
                                                                                                                      cGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTa                                <  1:2906/149‑1 (MQ=255)
                                                                                                                      agCTTCCTGTGGTTGTTGGAATGTTTGAGCGCAACTCCCACATTAAAGCAGAAACATCAGGCCGTCTGTGCAAGCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTa                                <  2:84224/148‑1 (MQ=255)
                                                                                                                       gCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAg                               >  4:12408/1‑149 (MQ=255)
                                                                                                                           ccTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCg                                  >  5:64808/1‑142 (MQ=255)
                                                                                                                           ccTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGCTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCg                                  <  6:64808/142‑1 (MQ=255)
                                                                                                                                  ttgttgGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATc                    >  7:268255/1‑149 (MQ=255)
                                                                                                                                     ttgGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCt                 <  6:274479/149‑1 (MQ=255)
                                                                                                                                       ggAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTg                <  3:96804/148‑1 (MQ=255)
                                                                                                                                        gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  3:195558/1‑149 (MQ=255)
                                                                                                                                        gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  8:14484/1‑149 (MQ=255)
                                                                                                                                        gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              <  3:83733/149‑1 (MQ=255)
                                                                                                                                        gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  2:280223/1‑149 (MQ=255)
                                                                                                                                              gggCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCAt        >  8:326657/1‑149 (MQ=255)
                                                                                                                                                gcgcgcAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAg           >  7:198588/1‑144 (MQ=255)
                                                                                                                                                gcgcgcAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAg           <  8:198588/144‑1 (MQ=255)
                                                                                                                                                  gcgcAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAg    <  1:201739/149‑1 (MQ=255)
                                                                                                                                                    gcAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAGcc  >  2:71163/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGTCGGCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAGCC  >  NC_002947/4640246‑4640542

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: