Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,214,658	0	A	G	22.6%	162.4 / 14.9	84	E538E (GAA→GAG)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (0/19);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 2.37e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCC  >  NC_000913/3214511‑3214806                                                                                                                                                    |                                                                                                                                                      gCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTNGNCCTGCCCCTGGATTCTNCGACCACCGaaa                                                                                                                                                      >  1:626637/1‑149 (MQ=255)    ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGcc                                                                                                                                                   >  1:1402355/1‑149 (MQ=255)         aCAAGATCCGCAAAGTGCTGAAGATCGCCACAGCGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATCCGCATCCGGGGGCTTTCATCGAGGATACCACCCCCGAGCCGCCGCTGGATTCTGCGACCCCCGAACGCCggcgg                                                                                                                                              >  2:780518/1‑148 (MQ=255)          cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtg                                                                                                                                             >  1:1406582/1‑149 (MQ=255)              aTCCGCAAAGTGCTGAAGATCGCCAAATAGCCAATCTCCATGGCAACGACTATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCTCCGAAAGCCTGCGTGCGGc                                                                                                                                         <  1:1127234/149‑1 (MQ=255)               tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGGTCGGTTATGATGAAGATTCGCATCTAGGGGATTTCATCGAGGATACCACCCTCGGGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTCCGGCa                                                                                                                                        >  2:1024436/1‑149 (MQ=255)                  gCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                                                                                                                                     <  1:545087/149‑1 (MQ=255)                   cAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCGGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGc                                                                                                                                    >  1:890376/1‑149 (MQ=255)                          cTGAAGATCGCCAAAGAGCCAATCTCCGTGGAACCGCCGATCGGGGATGATGAAGTTTCGCATCTGGGGGATTTCATCGAGGTTACCACCCTTGTTTTTTTTAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGt                                                                                                                             <  2:135645/149‑1 (MQ=255)                          cTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCCGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCACCCCACGCCGt                                                                                                                             >  2:147151/1‑149 (MQ=255)                            gAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGGgc                                                                                                                           >  1:700754/1‑149 (MQ=255)                              aGACCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATCTCAGCGAGGATACCCCCCTTGTTTTTTTTAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctg                                                                                                                         <  2:1549274/149‑1 (MQ=255)                                  cGCCAAAGAGCCGACGTCCATGGACGCGCCGATCGGTGATGATGTAGGTTCGTATCTGGGGGACTTCATCGAGGTTACCGCCCTTGTTTTTTTTCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                                                                     <  2:254738/149‑1 (MQ=255)                                  cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                                                                     <  1:975281/149‑1 (MQ=255)                                    ccAAAGAGCCATTCTCCAGGGAATCGCCGATCGGTTATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCGGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                                   <  1:932660/149‑1 (MQ=255)                                    ccAAAGAGCCAATCTTCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                                   <  1:1390063/149‑1 (MQ=255)                                               aTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACcgcgcg                                                                                                       <  1:627124/149‑1 (MQ=255)                                                 gtcCATGGACACGCCGATCTTTGTTTTGTGAGATTCGCATCTGCGGGATTGCATAGAGGATACCACCCTCGTTTTTTTTAATGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTg                                                                                                     <  2:546333/148‑1 (MQ=255)                                                    cATGGAAACGCCGATCGGTGTTGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTTTTGTTGTTGCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAg                                                                                                  <  2:1479539/149‑1 (MQ=255)                                                    cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAg                                                                                                  <  1:694325/149‑1 (MQ=255)                                                     aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGc                                                                                                 >  1:1477651/1‑149 (MQ=255)                                                      tGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCa                                                                                                >  1:530944/1‑149 (MQ=255)                                                         aaaCGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAGCGCACGACGTGCTGGCTGGCCAG‑ACCGCGCGTGAAGCaaaa                                                                                             >  2:1235008/1‑149 (MQ=255)                                                          acgcccGATGGGGGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCCCCCTCTTGTTTTCTAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAg                                                                                            <  2:172021/145‑1 (MQ=255)                                                             gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTc                                                                                         >  1:1262972/1‑149 (MQ=255)                                                             gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTc                                                                                         >  1:1233293/1‑149 (MQ=255)                                                                gctCGGTGGTGTGGCAGGTTGGCCTCTGGGGGATCTCCTCGAGGCTACCACCCTTGTGTTTCCTCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGc                                                                                      <  2:40919/147‑1 (MQ=255)                                                                 ttCGGTTATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCAACATAATCACTTTTGTTTCTGCGACCCCCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGCGAAGCAAAAGTTCTGCg                                                                                     <  1:362518/148‑1 (MQ=255)                                                                    ggTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCTGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTAt                                                                                  <  2:379514/149‑1 (MQ=255)                                                                     gTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATg                                                                                 <  1:443584/149‑1 (MQ=255)                                                                     gTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATg                                                                                 >  1:673922/1‑149 (MQ=255)                                                                      tgatgttgaAGGTGCGCATGTGGGGGTTTTCATCGGGGATACCACCCTCGAGCTGCTGCTTGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGc                                                                                <  2:626637/149‑1 (MQ=255)                                                                           atgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGCGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAACGTTCTGCGCATGCGttcc                                                                           >  2:1395976/1‑147 (MQ=255)                                                                             gaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTGTGCGTTTCgg                                                                         >  1:576743/1‑149 (MQ=255)                                                                             gaAGATTCGCATCTGGGGGATTTCATCCAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTCCGTCTGCGTTTcag                                                                         >  1:1303778/1‑147 (MQ=255)                                                                                gtttggCATCTGGGGGATTTCATCGAGGATACCCCCCTTGTGCTGTCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTAt                                                                      <  2:372742/144‑1 (MQ=255)                                                                                gATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGCCTGGCCTGTATCGTAC‑TCATGCAAAAGTTCCGCGTATGCGTTTCGGTAt                                                                      <  1:1145801/149‑1 (MQ=255)                                                                                  ttCGCATCTGGGGGATTTCATCGAGGATGCCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCg                                                                    >  1:1014143/1‑149 (MQ=255)                                                                                   tCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCCCCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCGGGCCGG‑ACCGCCCGTGAAGCAAACGTTCTCCGCATGCCTTTCGGtgcgag                                                                   >  2:1342383/1‑144 (MQ=255)                                                                                    ggCATGTGGGTGATTTCACCGAGGAGACCCCCCTTGTGTTTCCGCTGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGat                                                                  <  2:1485224/148‑1 (MQ=255)                                                                                     gCATCTGGGGGATTTCATCGAGGATACCACCCTTGTTTTGTTGATGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGata                                                                 <  2:527414/149‑1 (MQ=255)                                                                                     gCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAGGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTCAAGGAAAAGATCTGCGTACACGGTTCGGTATCGaca                                                                 >  2:861086/1‑147 (MQ=255)                                                                                     gCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGata                                                                 >  1:1242411/1‑149 (MQ=255)                                                                                      cATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGatat                                                                <  1:847021/149‑1 (MQ=255)                                                                                       aTCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATg                                                               <  2:1517851/149‑1 (MQ=255)                                                                                       aTCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGAGTATGCGTTTCGGTATCGATATg                                                               <  1:819373/149‑1 (MQ=255)                                                                                            tgggTTTTCATCCATGATTCCACCCTCTTGTTTTTGAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac                                                          <  2:155767/148‑1 (MQ=255)                                                                                            ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac                                                          <  1:1173397/149‑1 (MQ=255)                                                                                             ggaaTTTCATCGAGGATACCACCCTCTATTTGCCTCTTGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc                                                         <  2:647217/146‑1 (MQ=255)                                                                                                aTTTCATCGAGGCTCCCCTCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTCCGTATGCGTTTCGGTATCGATATGAACACCGAc                                                      <  1:994261/149‑1 (MQ=255)                                                                                                  ttCATCGCGGACACCCCCCTCGACCTGCCGCTGGATTCTGCGACCACCCAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTa                                                    <  1:763693/149‑1 (MQ=255)                                                                                                  ttCAGCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCCCCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTa                                                    <  1:764055/149‑1 (MQ=255)                                                                                                     aTCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCCCCGAAAGCCTGCGGGGGGCACCGGGAGACGTGCTGACAGCCCGG‑CCCCCGCGAGAAGAAAGAGTACAGCGTATGACTTTCCGTCTCAATATGCACACCCACGacca                                                 >  2:764775/1‑147 (MQ=255)                                                                                                           gATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGga                                           <  1:931675/149‑1 (MQ=255)                                                                                                                caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCGGCGTGCGGCAACGCACGACGTGCTGGCTGGCCGG‑ACCGCGCGTGAAGCCAAAGGTCTGCGTAGGCCGTTCGGTAGCGATATGAACCCCAACTACCCCCCGgaaagag                                      >  2:716351/1‑145 (MQ=255)                                                                                                                 accCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                     <  1:521366/149‑1 (MQ=255)                                                                                                                 accCTCGAGCTGCCGCTGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                     <  1:616021/149‑1 (MQ=255)                                                                                                                   ccTCGGGCTGCCGCTGGATTCTGCGACCACCGAAAGCCGGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAATACCGACTACACGCTGGAAGAAGTgg                                   <  1:780518/149‑1 (MQ=255)                                                                                                                   ccTCGGGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                                   <  2:1402355/149‑1 (MQ=255)                                                                                                                   ccTCGCGCTGCCCCTGGATTCGGCGACCACCGAGAGCCGGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                                   <  2:576743/149‑1 (MQ=255)                                                                                                                   ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                                   >  1:383449/1‑149 (MQ=255)                                                                                                                   ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                                   <  2:1014143/149‑1 (MQ=255)                                                                                                                   ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCGGCGCTCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGAGCGTGGACCAAAAGTTGTGCGTAGGGGGTGAGGTATCGATATGAACACCGACCACAAACAGGAAGAAgggg                                   >  2:1437445/1‑149 (MQ=255)                                                                                                                    cTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTggg                                  >  1:652047/1‑149 (MQ=255)                                                                                                                         tgttccGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAAc                             <  2:1233293/145‑1 (MQ=255)                                                                                                                            gCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt                          >  1:1200428/1‑149 (MQ=255)                                                                                                                             ccGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGtt                         >  1:265778/1‑149 (MQ=255)                                                                                                                              cGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCaa                                                                                                                                       >  2:184544/1‑39 (MQ=255)                                                                                                                              cGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCaa                                                                                                                                       <  1:184544/39‑1 (MQ=255)                                                                                                                              cGCTGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTc                        <  2:1542952/149‑1 (MQ=255)                                                                                                                                cTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGGCGTGCTGGCCGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGa                      <  2:673922/149‑1 (MQ=255)                                                                                                                                  ggATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACg                    <  1:915884/149‑1 (MQ=255)                                                                                                                                  ggATGCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACg                    <  2:973330/149‑1 (MQ=255)                                                                                                                                   gATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCgccgg                   >  2:121983/1‑145 (MQ=255)                                                                                                                                   gATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGt                   >  2:95917/1‑149 (MQ=255)                                                                                                                                     ttCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACAAGCTGGAAGAAGTGGGTAAACAGTTCGACGTTa                 <  2:530944/149‑1 (MQ=255)                                                                                                                                     ttCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTa                 >  1:1388900/1‑149 (MQ=255)                                                                                                                                      tCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAc                >  1:515967/1‑149 (MQ=255)                                                                                                                                        tGCGACCACCGAAAGCCTGCGTGCGGCAGAAGTCGAGGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc              >  2:1107192/1‑149 (MQ=255)                                                                                                                                        tGCCACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc              <  2:691787/149‑1 (MQ=255)                                                                                                                                           gACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCcgcg           >  1:213147/1‑149 (MQ=255)                                                                                                                                             ccaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGaa         <  1:406761/149‑1 (MQ=255)                                                                                                                                                ccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATACGTTGCGGTCTCGATATGAACACGGACTACACGATGGTAGAACAGAGTAACCGGGTCGGCGTTACCagcggccgt      >  1:1486881/1‑144 (MQ=255)                                                                                                                                                 cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGGAAACAGTTCGACGTTACCCGCGAACGTa     >  1:1005817/1‑149 (MQ=255)                                                                                                                                                  gAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAt    >  1:283335/1‑149 (MQ=255)                                                                                                                                                   aaaGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATc   >  2:843552/1‑149 (MQ=255)                                                                                                                                                    aaGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATcc  <  2:1363611/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTG‑ACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCC  >  NC_000913/3214511‑3214806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑