Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,819,831 Δ1 bp 5.7% intergenic (+26/‑449) rpoS → / → PP_1624 RNA polymerase sigma S factor/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,819,8310G. 5.7% 761.5 / 39.1 228intergenic (+26/‑449)rpoS/PP_1624RNA polymerase sigma S factor/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (104/111);  new base . (9/4);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 230 total aligned reads displayed.

GATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1819683‑1819978
                                                                                                                                                    |                                                                                                                                                      
gaTGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCgg                                                                                                                                                       >  1:4564990/1‑150 (MQ=255)
    ttGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGAt                                                                                                                                                   >  1:2811052/1‑150 (MQ=255)
       ggCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGAtc‑t‑c                                                                                                                                              >  1:1670438/1‑150 (MQ=255)
        gCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGGGGAAGGGCTCCAGCGCTTGCGCGAGAGCCTTGAAAAGAATGGTCTGTCCCGCGAGTCCGTGTGCCAGTAGCTCGCGCGGTAACACAACACGCCCCGGTGAGC‑T‑Cg                                                                                                                                             >  2:3954832/1‑150 (MQ=255)
         cTGGAGCTCGGCCTTACCCGTGAACGGGTGCGGCAAAGCCAGGTGGACGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑Cgg                                                                                                                                            <  2:208803/150‑1 (MQ=255)
         cTGGAGAGCGGCCTTACCCGTGAACGGGTGCGGCAAATCCTGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑Cgg                                                                                                                                            <  2:3785907/150‑1 (MQ=255)
          tGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑Cggg                                                                                                                                           >  1:4356448/1‑150 (MQ=255)
            gagaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCACGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGCGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGc                                                                                                                                         <  1:2171391/150‑1 (MQ=255)
              gaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGACATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑t                                                                                                                                      <  2:2931438/150‑1 (MQ=255)
               aTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑tt                                                                                                                                     >  2:2621755/1‑150 (MQ=255)
                  ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑ttttt                                                                                                                                  >  1:1522633/1‑150 (MQ=255)
                   gCCTTACCCGTGAACGGGTGCGGCAAATCCTGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑ttttt                                                                                                                                  <  1:4401432/149‑1 (MQ=255)
                   gCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑ttttt                                                                                                                                  >  2:4401432/1‑149 (MQ=255)
                   gCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑ttgttg                                                                                                                                 >  2:4821426/1‑150 (MQ=255)
                     cTTACCGGTGAACGGGTGCGGCATATCCTGGTGGAAGTGCTGAAGCGCTTGCGCTAGATCCTTGAAAAGATTGGTCTGTCCTTCTTTTCATTGTTCCAGTAGCTAGCTCTTTAAAAAAACACGCCCC‑GAGATCAC‑CGGGGCG‑ttttttgt                                                                                                                               <  2:3720494/150‑4 (MQ=255)
                       tACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTtgtgtg                                                                                                                             >  1:3259830/1‑150 (MQ=255)
                          ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTgcgc                                                                                                                          >  1:1110886/1‑150 (MQ=255)
                               caCGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGc                                                                                                                     <  2:2542182/149‑1 (MQ=255)
                                  gggTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATg                                                                                                                  <  2:3617488/150‑1 (MQ=255)
                                   ggTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATgg                                                                                                                 >  1:3703615/1‑150 (MQ=255)
                                      gCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGgcgc                                                                                                              >  2:2907419/1‑150 (MQ=255)
                                           aaaTCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                         <  2:3307684/150‑1 (MQ=255)
                                             aTCCAGGTGGAAGGGCTCAAGCGCTGGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTgcgc                                                                                                       <  2:4225843/150‑1 (MQ=255)
                                              tCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGcgcg                                                                                                      >  2:3207003/1‑150 (MQ=255)
                                              tCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGcgcg                                                                                                      >  1:1791563/1‑150 (MQ=255)
                                               ccAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGt                                                                                                     >  1:588412/1‑150 (MQ=255)
                                                cAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTa                                                                                                    <  2:5022731/150‑1 (MQ=255)
                                                 aGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                   >  2:2495332/1‑150 (MQ=255)
                                                    tGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcg                                                                                                >  1:633535/1‑150 (MQ=255)
                                                       aaGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGAGCGTATGCGCAA                                                                                             >  1:82172‑M1/1‑150 (MQ=255)
                                                        aGGGCTCACGCGCTTGCGCGAGATCCTTTAAAATAATGGTCTGTCATGCGTGTCTTTGTTCCAGTAGCTACCGCTTTAACAAAACACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAC                                                                                            <  1:4622595‑M1/150‑1 (MQ=255)
                                                        aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAC                                                                                            <  2:841826‑M1/150‑1 (MQ=255)
                                                           gCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCg                                                                                         >  2:2323064/1‑150 (MQ=255)
                                                           gCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCg                                                                                         >  1:5228943/1‑150 (MQ=255)
                                                               aaGCGCTTGCGCGAGATCCTTGAAAAGACTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGGATGCGCAACCCGCTTg                                                                                     >  1:693324/1‑150 (MQ=255)
                                                                 gcgcTTGCGCGAGATCCTTGAAAACAATGGTCGGTCCAGCGAGTCATTGTTCCAGTCGCTAGCGCGGTAACAAAACACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGt                                                                                   <  2:1401941/150‑1 (MQ=38)
                                                                    cTTGCGCGATACCCTTGAAAAGACTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                                <  1:677725/150‑1 (MQ=255)
                                                                     ttGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgggggtg                                                                               >  2:2777630/1‑145 (MQ=255)
                                                                        cgcgAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                            <  2:2121700/150‑1 (MQ=255)
                                                                         gcgAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                           >  2:4108071/1‑150 (MQ=255)
                                                                          ctagaTCCTTTAAAATAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGc                                                                          <  1:2001374/148‑1 (MQ=255)
                                                                               tCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCTGTAGCTAGCGCTTTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                     <  2:4273365/150‑1 (MQ=255)
                                                                                ccTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTg                                                                    <  2:4157626/150‑1 (MQ=255)
                                                                                ccTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTg                                                                    <  1:3542936/150‑1 (MQ=255)
                                                                                ccTTCAAAATAATGGTCTGTCCAGCGAGGCATTGTTCCAGTAGCTACCGCTGTAACAAAAAACGCCCCGGTGATC‑A‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTAGCCAAGTg                                                                    <  1:3950621/150‑1 (MQ=25)
                                                                                 cTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                   >  2:2184263/1‑150 (MQ=255)
                                                                                  ttGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTggg                                                                  >  1:1267473/1‑150 (MQ=255)
                                                                                   tGAAAAGAAGGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                 <  1:5350896/150‑1 (MQ=255)
                                                                                    gAAAAGAATGGTCTGTCCAGCTTGTAATTGTTCCAGTAGCTTGCGATGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCacttggccaccacaaccaagcgggtt                                                                <  1:3119535‑M1/150‑27 (MQ=255)
                                                                                       aaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACt                                                             >  1:2134517/1‑150 (MQ=255)
                                                                                          aaTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                          >  1:3484959/1‑150 (MQ=255)
                                                                                          aaTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                          <  1:3413805/150‑1 (MQ=255)
                                                                                           aTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTTGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                         <  1:679813/150‑1 (MQ=255)
                                                                                           aTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                         <  2:2658924/150‑1 (MQ=255)
                                                                                              gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGCTC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                      <  2:4611949/150‑1 (MQ=255)
                                                                                              gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                      <  1:4535024/150‑1 (MQ=255)
                                                                                                 tGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaa                                                   >  1:2418202/1‑150 (MQ=255)
                                                                                                 tGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaa                                                   >  1:117420/1‑150 (MQ=255)
                                                                                                  gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  1:2360107/150‑1 (MQ=255)
                                                                                                  gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  2:82172/150‑1 (MQ=255)
                                                                                                     cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                               <  2:1670438/150‑1 (MQ=255)
                                                                                                       gCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                             <  1:4021526/150‑1 (MQ=255)
                                                                                                           gTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCACCCCGCGTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTac                                         >  2:3112538/1‑150 (MQ=255)
                                                                                                           gTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTac                                         <  2:2418202/150‑1 (MQ=255)
                                                                                                            tCATTGTTCCAGTAGCTAGCGCTGTAAAAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                        <  1:2535452/150‑1 (MQ=255)
                                                                                                              aTTGTTCCAGTAGCTAGCGCTGTAACAAAACACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTac                                         <  2:1379313/147‑1 (MQ=255)
                                                                                                              aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacaca                                      >  1:5075350/1‑150 (MQ=255)
                                                                                                              aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTac                                         >  1:1379313/1‑147 (MQ=255)
                                                                                                                  ttCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCg                                  >  1:3890525/1‑150 (MQ=255)
                                                                                                                    ccAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                                <  2:668652/150‑1 (MQ=255)
                                                                                                                      aGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCac                              >  2:2579996/1‑149 (MQ=255)
                                                                                                                      aGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa                              <  2:1175056/150‑1 (MQ=255)
                                                                                                                       gTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCAAGGCGCGGTGCGCGTATGCGCAACCCCCGTGGTGGTGGTGGCCAAGTGGGTGACGTGGTGGTGGAAGTCCTCTACACACCCGGCAAg                             >  1:258166/1‑150 (MQ=18)
                                                                                                                        tagctagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                            <  2:3938428/150‑1 (MQ=255)
                                                                                                                          gctcgcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                          <  1:3954832/150‑1 (MQ=255)
                                                                                                                            tcgcgcTGTAACAAAACACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                        <  1:130225/148‑1 (MQ=255)
                                                                                                                            tagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGTGGTGAAAGTCCTCTACACACCCGGCAAggggca                        >  1:5014256/1‑148 (MQ=255)
                                                                                                                            tagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                        <  1:1245096/150‑1 (MQ=255)
                                                                                                                                  tGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa                  <  1:1197518/150‑1 (MQ=255)
                                                                                                                                   gcacacaaaaaCGCCCC‑GAGATCAC‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                 >  1:3027118/7‑150 (MQ=11)
                                                                                                                                    tAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                        <  2:2019378/110‑1 (MQ=255)
                                                                                                                                     acacaaaaaCGCCCC‑GAGATCAC‑CGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                >  2:4763832/5‑150 (MQ=11)
                                                                                                                                     acacaaaaaCGCCCC‑GAGATCAC‑CGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCGTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGACGTGTTAGc                >  1:3928348/5‑150 (MQ=11)
                                                                                                                                     aaCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGcc               <  1:4762365/150‑1 (MQ=255)
                                                                                                                                      aCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                    >  2:3190103/1‑144 (MQ=255)
                                                                                                                                      aCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                    <  1:3190103/144‑1 (MQ=255)
                                                                                                                                       acaaaaaCGCCCC‑GAGATCAC‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  >  1:4880048/3‑113 (MQ=11)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCggc            >  2:3807124/1‑149 (MQ=255)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            >  1:4040724/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            >  2:1515087/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            <  2:1749259/150‑1 (MQ=255)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            <  1:4242288/150‑1 (MQ=255)
                                                                                                                                        aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            <  2:3161916/150‑1 (MQ=255)
                                                                                                                                          aaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTGGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCACGGGGACGTGGTAGCCGGAgg          >  2:4637120/1‑150 (MQ=255)
                                                                                                                                           aaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc         >  1:2864425/1‑150 (MQ=255)
                                                                                                                                           aaaCGCCCC‑GAGATC‑ACCGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc         >  1:697778/1‑150 (MQ=11)
                                                                                                                                            aaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCa        >  2:5173626/1‑150 (MQ=255)
                                                                                                                                             aCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa       >  1:881840/1‑150 (MQ=255)
                                                                                                                                               gCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                <  1:2090695/139‑1 (MQ=255)
                                                                                                                                                  ccGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtg  >  2:4096481/1‑150 (MQ=40)
                                                                                                                                                    |                                                                                                                                                      
GATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1819683‑1819978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑