Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,748,925 T→A 7.4% intergenic (‑258/‑74) PP_1548 ← / → PP_5744 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,748,9250TA 7.4% 186.9 / 3.2 63intergenic (‑258/‑74)PP_1548/PP_5744hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GCTACGTTGCGTAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGG  >  NC_002947/1748780‑1749065
                                                                                                                                                  |                                                                                                                                             
gCTACGTTGCGTAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGTCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGt                                                                                                                                           >  5:272733/1‑149 (MQ=255)
gCTACGTTGCGTAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGt                                                                                                                                           >  7:159888/1‑149 (MQ=255)
  tACGTTGCGTAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGttt                                                                                                                                         >  6:109525/1‑149 (MQ=255)
           tAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATg                                                                                                                                 >  1:286404/1‑148 (MQ=255)
           tAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATg                                                                                                                                 <  2:286404/148‑1 (MQ=255)
            aaCGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGc                                                                                                                                >  2:297202/1‑148 (MQ=255)
             aCGTCCCCCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTGTTTAAAGCGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGc                                                                                                                                >  5:230536/1‑147 (MQ=255)
             aCGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGGGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCaa                                                                                                                              <  5:307114/149‑1 (MQ=255)
             aCGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCaa                                                                                                                              <  2:275554/149‑1 (MQ=255)
               gTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTAATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAACGAACTAGCATTTATAAAAGCATGAATGTTTTTTAAAAGCAAg                                                                                                                             >  4:299594/1‑148 (MQ=37)
                     caGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGTGCTTGCATTTAAAAAACAAGCATGCTTTTATAAATGCAAGTTCGCtt                                                                                                                      <  5:221018/149‑1 (MQ=25)
                       gAAGTCAATCGCCGGATTAACTCGTCGCCTGGTGGCTCATCTGCTGTTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCTAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTtgt                                                                                                                    <  3:299594/149‑1 (MQ=255)
                         aCTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGATCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGTGCTTGCATTTAAAAAACAAGCATGCTTTTATAAATGCAAGTTCGCTtgcat                                                                                                                  <  6:167323/149‑4 (MQ=21)
                           tCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTaa                                                                                                                >  8:280292/1‑149 (MQ=255)
                                  ccGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTtctc                                                                                                         <  8:159888/149‑1 (MQ=255)
                                  ccGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTtctc                                                                                                         >  6:297406/1‑149 (MQ=255)
                                     gATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTa                                                                                                      <  1:233335/149‑1 (MQ=255)
                                       ttAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACt                                                                                                    >  7:116113/1‑149 (MQ=255)
                                          aCTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGTGCTTGCATTTAAAAAACAAGCATGCTTTTATAAATGCAAGTTCGCTtgca                                                                                                                   <  6:133968/131‑3 (MQ=18)
                                          aCTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGTGCTTGCATTTAAAAAACAAGCATGCTTTTATAAATGCAAGTTCGCTtgca                                                                                                                   >  5:133968/1‑129 (MQ=18)
                                          aCTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCTCGCTTGTGTAATTTTCTCTTACTGTa                                                                                                 >  6:142336/1‑149 (MQ=255)
                                            tCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATg                                                                                               >  3:271337/1‑149 (MQ=255)
                                            tCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATg                                                                                               >  1:161262/1‑149 (MQ=255)
                                              ggCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTTCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATCGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTAt                                                                                                <  1:240294/147‑1 (MQ=255)
                                                          cTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGTCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAgtgtgt                                                                                 <  6:272733/149‑1 (MQ=255)
                                                             aTCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTa                                                                                <  7:280292/147‑1 (MQ=255)
                                                                 gctgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGTAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCt                                                                          >  8:132780/1‑149 (MQ=255)
                                                                  ctgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGa                                                                             <  4:271337/145‑1 (MQ=255)
                                                                  ctgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGa                                                                             <  2:198190/145‑1 (MQ=255)
                                                                   tgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCttt                                                                        >  2:313334/1‑149 (MQ=255)
                                                                     gTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCg                                                                      <  5:109525/149‑1 (MQ=255)
                                                                     gTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCg                                                                      >  8:140428/1‑149 (MQ=255)
                                                                        tttAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGa                                                                    <  7:112074/148‑1 (MQ=255)
                                                                                    tGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAgg                                                        >  2:66202/1‑148 (MQ=255)
                                                                                         gCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAagag                                                  <  6:230536/149‑1 (MQ=255)
                                                                                                   ccAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTa                                         <  3:209397/148‑1 (MQ=255)
                                                                                                   ccAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTAt                                        >  4:288188/1‑149 (MQ=255)
                                                                                                     aGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATg                                       >  5:247392/1‑148 (MQ=255)
                                                                                                       cTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGaa                                    <  7:140428/149‑1 (MQ=255)
                                                                                                          gATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAagag                                 >  7:222541/1‑149 (MQ=255)
                                                                                                          gATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAaga                                  <  3:288188/148‑1 (MQ=255)
                                                                                                                aatgcaaGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAAc                            >  5:280518/1‑148 (MQ=255)
                                                                                                                  tgcaaGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAg                         >  6:81033/1‑149 (MQ=255)
                                                                                                                        cGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCTCGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGa                   <  5:142336/149‑1 (MQ=255)
                                                                                                                         gAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAg                  <  1:297202/149‑1 (MQ=255)
                                                                                                                         gAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAg                  <  1:66202/149‑1 (MQ=255)
                                                                                                                          aaCTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGc                 <  2:258308/149‑1 (MQ=255)
                                                                                                                             ttGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTAt                                        >  3:246567/1‑123 (MQ=255)
                                                                                                                             ttGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTAt                                        <  4:246567/123‑1 (MQ=255)
                                                                                                                             ttGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGc              >  6:242800/1‑149 (MQ=255)
                                                                                                                               gCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTg            >  6:236964/1‑149 (MQ=255)
                                                                                                                                cATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTgg           <  5:297406/149‑1 (MQ=255)
                                                                                                                                 aTTTATAATAGCTTGCTTGTTTTTTAAATGCAAGCACGGTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGATAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTgg           <  5:28557/148‑1 (MQ=255)
                                                                                                                                 aTTTATAAAAGCATGCTTGATTTTTAAATGCAAGCACGCTTGTGTAATTTTCACTTACTGTAAGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGCCGCAGGAAAGAGCAGTTAGTATGGAAGAGCTAACGAGCCTTGAGCCGCTgg           >  6:28557/1‑148 (MQ=255)
                                                                                                                                  tttATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGc         <  7:132780/149‑1 (MQ=255)
                                                                                                                                  tttATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGc         >  6:282359/1‑149 (MQ=255)
                                                                                                                                   ttATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCt        >  3:62091/1‑149 (MQ=255)
                                                                                                                                    tataAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCtt       >  5:97459/1‑149 (MQ=255)
                                                                                                                                      taAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGa     <  2:161262/149‑1 (MQ=255)
                                                                                                                                        aaaGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAgg   <  1:313334/149‑1 (MQ=255)
                                                                                                                                         aaGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAggg  >  1:158770/1‑149 (MQ=255)
                                                                                                                                         aaGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTACGAGCCTTGAGCGGCTGGGCTTGAgg   <  4:9933/147‑1 (MQ=255)
                                                                                                                                         aaGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTACGAGCCTTGAGCGGCTGGGCTTGAgg   >  3:9933/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCTACGTTGCGTAACGTCCACCAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGG  >  NC_002947/1748780‑1749065

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: