Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,640,397 A→G 7.7% intergenic (‑35/+23) 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,3970AG 7.7% 182.2 / 4.9 65intergenic (‑35/+23)PP_4105/gabP‑IVhypothetical protein/gamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (2/3);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAGCC  >  NC_002947/4640253‑4640542
                                                                                                                                                |                                                                                                                                                 
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                                                                                                  <  1:227819/149‑1 (MQ=255)
                                             ccGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGc                                                                                                  <  4:211171/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:83733/149‑1 (MQ=255)
                                                                                                                                 gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  8:14484/1‑149 (MQ=255)
                                                                                                                                 gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  2:280223/1‑149 (MQ=255)
                                                                                                                                 gAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGaa              >  3:195558/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)
                                                                                                                                                aCGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCAt        <  1:239501/140‑1 (MQ=255)
                                                                                                                                                aCGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCAt        >  2:239501/1‑140 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCGTCCCAGTGGCTTTGCAGGTGCAGGCCGCTGGGGTGGGAACAGCCGCAGCCATAGCAAGTACCCTCAGGGGCCGCGAGGGCTTGAAGGGAAAGTGTGTGATCAGTCATACGCTTCCTGTGCTTGTTGGAATGGGGGCGCGCAACGCCCCCATTAAAGCAGAAACATCAGGCCGTCTGGGCAACCATCTTCGCACGCCGCTGGGTCAGCGAAACCAGCAAAATGAACGCGGTCACCCCGGCAGTCATCAGCGTTTCGTAGCGATAGGCATCGCTGAACAGCATGTAGCC  >  NC_002947/4640253‑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: