Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,005 G→A 5.6% R81H (CGC→CAC) ‡ rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,0050GA 5.6% 677.7 / 18.3 215R81H (CGC→CAC) ‡rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (128/75);  new base A (7/5);  total (135/80)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 221 total aligned reads displayed.

GTGGTCGTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGA  >  NC_002947/541859‑542154
                                                                                                                                                  |                                                                                                                                                     
gTGGTCGTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCgcgg                                                                                                                                                    >  2:615881/1‑150 (MQ=255)
      gTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAAt                                                                                                                                              >  1:3574106/1‑150 (MQ=255)
          ggTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGc                                                                                                                                          <  1:4937672/150‑1 (MQ=255)
          ggTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGc                                                                                                                                          >  1:4320298/1‑150 (MQ=255)
            tGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCga                                                                                                                                        <  2:694063/150‑1 (MQ=255)
            tGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCga                                                                                                                                        <  2:637418/150‑1 (MQ=255)
                 ttAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGt                                                                                                                                   >  2:2523413/1‑150 (MQ=255)
                   aagaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgc                                                                                                                                 >  2:550749/1‑150 (MQ=255)
                     gaagCCGGAGACGATCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTgcgg                                                                                                                               <  2:2352132/150‑1 (MQ=255)
                     gaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcgg                                                                                                                               >  1:2094814/1‑150 (MQ=255)
                        gCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGt                                                                                                                            >  1:2809999/1‑150 (MQ=255)
                         ccGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCACGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGtt                                                                                                                           >  2:3933832/1‑150 (MQ=255)
                           gtaaaCCATCAACTACGGTAGGTTCAAGACTGAGCGTGACTGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGACTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAGGTGCGGCGTTgt                                                                                                                         >  2:3380578/3‑149 (MQ=255)
                           ggAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGa                                                                                                                         >  2:2899830/1‑150 (MQ=255)
                            gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                        >  1:2452419/1‑150 (MQ=255)
                            gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                        <  1:394584/150‑1 (MQ=255)
                                 cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  2:4594766/1‑150 (MQ=255)
                                 cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  2:4549246/1‑150 (MQ=255)
                                 cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   <  1:1870162/150‑1 (MQ=255)
                                 cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  1:3727140/1‑150 (MQ=255)
                                 cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  2:4281109/1‑150 (MQ=255)
                                  aTCAACTACCGTACGTTAAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGcc                                                                                                                  >  2:4878692/1‑150 (MQ=255)
                                   tCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGccc                                                                                                                 >  2:1706104/1‑150 (MQ=255)
                                                  tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                  <  2:3063773/150‑1 (MQ=255)
                                                  tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                  <  1:3073476/150‑1 (MQ=255)
                                                  tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                  >  1:4302273/1‑150 (MQ=255)
                                                   cAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtcg                                                                                                 >  1:2420117/1‑150 (MQ=255)
                                                    aaGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTcgtcgt                                                                                                <  1:436264/150‑1 (MQ=255)
                                                    aaGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTcgtcgt                                                                                                >  1:2771896/1‑150 (MQ=255)
                                                      gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              <  1:4139503/150‑1 (MQ=255)
                                                      gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              <  1:3224855/150‑1 (MQ=255)
                                                      gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              <  1:3105031/150‑1 (MQ=255)
                                                      gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              >  1:1089079/1‑150 (MQ=255)
                                                        cTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAgc                                                                                            >  2:3603291/1‑150 (MQ=255)
                                                         tGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgcg                                                                                           >  2:2149021/1‑150 (MQ=255)
                                                         tGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgcg                                                                                           >  1:3890132/1‑150 (MQ=255)
                                                           aGCGTGACGGCCTGATCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTAAGAGTGACTGAGCGGCAAGTACAAGCGCCTCAACCACCGCGGCGAAATCTGCTAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGATCGTCGTGAGCGCa                                                                                         >  2:319419/1‑150 (MQ=255)
                                                                       tGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCg                                                                             >  1:4848246/1‑150 (MQ=255)
                                                                        gTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCCCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGa                                                                            >  1:2854351/1‑150 (MQ=255)
                                                                         ttCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAg                                                                           >  1:43284/1‑150 (MQ=255)
                                                                           cTGCGCCACGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCt                                                                         <  2:436807/150‑1 (MQ=255)
                                                                           cTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCt                                                                         <  2:1916581/150‑1 (MQ=255)
                                                                              cgcgcgtaTCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGc                                                                      <  2:707491/143‑1 (MQ=255)
                                                                               gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGcc                                                                     >  1:4767111/1‑150 (MQ=255)
                                                                               gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGcc                                                                     >  2:3910735/1‑150 (MQ=255)
                                                                                ccAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCt                                                                    >  2:3181185/1‑150 (MQ=255)
                                                                                ccAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCt                                                                    >  2:1258346/1‑150 (MQ=255)
                                                                                 cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                   >  1:224057/1‑150 (MQ=255)
                                                                                  aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                  >  2:1889494/1‑150 (MQ=255)
                                                                                  aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                  <  2:1955696/150‑1 (MQ=255)
                                                                                   aGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGc                                                                 >  1:2395644/1‑150 (MQ=255)
                                                                                        tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                            >  1:435347/1‑150 (MQ=255)
                                                                                        tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                            <  1:4214573/150‑1 (MQ=255)
                                                                                             cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTAGCCca                                                       >  1:2944669/1‑150 (MQ=255)
                                                                                              ccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCcac                                                      >  2:3035637/1‑150 (MQ=255)
                                                                                                aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAt                                                    >  2:2486776/1‑150 (MQ=255)
                                                                                                aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAt                                                    <  2:4188601/150‑1 (MQ=255)
                                                                                                 gTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATc                                                   <  2:4622653/150‑1 (MQ=255)
                                                                                                    aaGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTgg                                                >  1:4749015/1‑150 (MQ=255)
                                                                                                    aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  1:1449361/1‑150 (MQ=255)
                                                                                                    aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                <  2:1331379/150‑1 (MQ=255)
                                                                                                    aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  1:2190457/1‑150 (MQ=255)
                                                                                                    aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGACCACATCTgg                                                >  2:4715061/1‑150 (MQ=255)
                                                                                                     aGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGt                                               >  2:2747295/1‑150 (MQ=255)
                                                                                                       gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                             <  1:4878692/150‑1 (MQ=255)
                                                                                                         cTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCt                                           <  2:3504269/150‑1 (MQ=255)
                                                                                                         cTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCt                                           >  1:740884/1‑150 (MQ=255)
                                                                                                          tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                          >  2:2134638/1‑150 (MQ=255)
                                                                                                          tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                          >  2:3803463/1‑150 (MQ=255)
                                                                                                          tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                          >  1:2341981/1‑150 (MQ=255)
                                                                                                            cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                        <  2:4114242/150‑1 (MQ=255)
                                                                                                             gaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAg                                       >  2:1625166/4‑150 (MQ=255)
                                                                                                              aGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGt                                      >  2:4988533/1‑150 (MQ=255)
                                                                                                                  ccTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCt                                  <  1:1570223/150‑1 (MQ=255)
                                                                                                                    tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCACATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGc                                >  1:1518651/1‑150 (MQ=255)
                                                                                                                     gtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGcc                               >  1:2470417/1‑150 (MQ=255)
                                                                                                                        cGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTc                            <  1:3035637/150‑1 (MQ=255)
                                                                                                                        cGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTc                            <  1:3449996/150‑1 (MQ=255)
                                                                                                                          gCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATcgc                          <  1:5223/150‑3 (MQ=255)
                                                                                                                          gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCAt                                                                                        >  1:2879627/1‑88 (MQ=255)
                                                                                                                          gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                          <  2:798312/150‑1 (MQ=255)
                                                                                                                              gTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTAt                      >  2:3299153/1‑150 (MQ=255)
                                                                                                                               tACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATc                     >  2:2219174/1‑150 (MQ=255)
                                                                                                                                aCAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCg                    <  1:673350/150‑1 (MQ=255)
                                                                                                                                 cAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCgg                   >  2:1687737/1‑150 (MQ=255)
                                                                                                                                 cAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCgg                   <  1:555483/150‑1 (MQ=255)
                                                                                                                                  aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  <  1:4432747/150‑1 (MQ=255)
                                                                                                                                  aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  >  2:4099119/1‑150 (MQ=255)
                                                                                                                                  aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  >  1:2283611/1‑150 (MQ=255)
                                                                                                                                     cgcCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCACCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTg               <  2:2464807/150‑1 (MQ=255)
                                                                                                                                      gcCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc              >  1:1639655/1‑150 (MQ=255)
                                                                                                                                         tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa           <  1:4549246/150‑1 (MQ=255)
                                                                                                                                         tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa           >  2:1096000/1‑150 (MQ=255)
                                                                                                                                             gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa       >  2:3313904/1‑150 (MQ=255)
                                                                                                                                             gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa       >  1:4650259/1‑150 (MQ=255)
                                                                                                                                              cACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAc      <  2:1019140/150‑1 (MQ=255)
                                                                                                                                               aCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACa     >  2:4938237/1‑150 (MQ=255)
                                                                                                                                                  gcggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa  >  2:394160/1‑150 (MQ=255)
                                                                                                                                                  gcggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa  <  2:414998/150‑1 (MQ=255)
                                                                                                                                                  gcggcgTAATCTGCGAGAAGTGCAGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa  >  2:1756539/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGTCGTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGA  >  NC_002947/541859‑542154

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