Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,748,944 A→T 9.4% intergenic (‑277/‑55) 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,9440AT 9.4% 174.6 / 8.5 64intergenic (‑277/‑55)PP_1548/PP_5744hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base T (2/4);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTCGCCGCTGG  >  NC_002947/1748801‑1749092
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                          >  1:161262/1‑149 (MQ=255)
                       tCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATg                                                                                                                          >  3:271337/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                                                                                                        <  2:198190/145‑1 (MQ=255)
                                             ctgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGa                                                                                                        <  4:271337/145‑1 (MQ=255)
                                              tgGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCttt                                                                                                   >  2:313334/1‑149 (MQ=255)
                                                gTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCg                                                                                                 >  8:140428/1‑149 (MQ=255)
                                                gTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCg                                                                                                 <  5:109525/149‑1 (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:66202/149‑1 (MQ=255)
                                                                                                    gAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAg                                             <  1:297202/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                                    >  6:282359/1‑149 (MQ=255)
                                                                                                             tttATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGc                                    <  7:132780/149‑1 (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)
                                                                                                                               gTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGAAGTTAGTATGGAAGAGCTAACGAGACTTGAGCTGATGGGCTTGATGGTTTCAGCgat                   >  8:93804/1‑149 (MQ=255)
                                                                                                                                   ttttAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCgatga                 <  6:280518/148‑1 (MQ=255)
                                                                                                                                    tttAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCgatgat                <  7:233602/148‑1 (MQ=255)
                                                                                                                                     ttAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATc               <  5:282359/148‑1 (MQ=255)
                                                                                                                                       aaaTGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCATTGTGTAGGAGCTTTCGGATGTCGGCGCTCGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCaa             <  7:93804/148‑1 (MQ=255)
                                                                                                                                       aaaTGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCaa             <  5:81033/148‑1 (MQ=255)
                                                                                                                                            cAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAgg                              >  2:35355/1‑126 (MQ=255)
                                                                                                                                            cAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAgg                              <  1:35355/126‑1 (MQ=255)
                                                                                                                                            cAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTcgc        <  8:26150/148‑1 (MQ=255)
                                                                                                                                            cAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTcgc        <  8:31019/148‑1 (MQ=255)
                                                                                                                                                cACGCTTGTGTAATTTTCTCTTACTGTATGGAGGTCCAGTGTGTAGGAGCTTTCGGATGTCTGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTCGCCGCTg   <  8:116113/149‑1 (MQ=255)
                                                                                                                                                 aCGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTCGCCGCTgg  <  6:247392/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAGAACTCAATCGCCGGATTAACTCGGCGCCTGGTGGCTCATCTGCTGGTTTTTAAAGAGCGGTGGCTGCCGGAGCTGCCAGCTCGATGCAAATGCAAGCGAACTTGCATTTATAAAAGCATGCTTGTTTTTTAAATGCAAGCACGCTTGTGTAATTTTCTCTTACTGTATGGATGTCCAGTGTGTAGGAGCTTTCGGATGTCGGCGCAGGAAAGAGCAGTTCGTATGGAAGAGCTAACGAGCCTTGAGCGGCTGGGCTTGAGGGTTTCAGCGATGATCAACTCGCCGCTGG  >  NC_002947/1748801‑1749092

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: