Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,173	0	T	G	21.7%	135.9 / 12.0	60	G31G (GGT→GGG)	ymdF	conserved protein YmdF
Reads supporting (aligned to +/- strand):  ref base T (17/30);  new base G (13/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 4.64e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACC  >  NC_000913/1068034‑1068311                                                                                                                                            |                                                                                                                                           cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                                            <  2:519025/140‑1 (MQ=255) cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                                            <  2:294391/140‑1 (MQ=255)             gCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAAt                                                                                                                                <  1:871167/140‑1 (MQ=255)                 tatGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGAcc                                                                                                                            >  2:596265/1‑140 (MQ=255)                 tatGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTACGGGCATTTTGTAGGAGACCGCGAAAGAGGATCAAGAGCACGTAAAAAAAGGGGAGAAAACAACAGGGGTCAGTTTAAAAATGAcc                                                                                                                            >  2:565289/1‑140 (MQ=255)                      tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTTCGGCAATTTTTCAGAAGACCGCGAGAAAGCATCAGAACCAGGTAAAAAAGGTGGAAAACACAGCGGGGGGAATTTCAAGAATTACCCGCAg                                                                                                                       >  2:1070920/1‑140 (MQ=255)                      tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAg                                                                                                                       >  2:1199580/1‑140 (MQ=255)                        tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAATACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcg                                                                                                                     <  2:1074884/140‑1 (MQ=255)                        tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcg                                                                                                                     <  2:828650/140‑1 (MQ=255)                         cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCCGAAGACCGCGAAAAAGCATCAGAAGCAGGGAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAgcgc                                                                                                                    >  2:227074/1‑140 (MQ=255)                         cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgc                                                                                                                    <  2:435852/140‑1 (MQ=255)                         cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgc                                                                                                                    <  1:1199580/140‑1 (MQ=255)                           ataCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAgcgcgc                                                                                                                  >  2:458278/1‑140 (MQ=255)                           ataCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAAAGACCCGCAgcgcgc                                                                                                                  >  2:27552/1‑140 (MQ=255)                            taCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa                                                                                                                 >  1:1130312/1‑140 (MQ=255)                            taCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCCCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAACACCGCGGGGGTAAATTCAAAAATAACCCTCCgcgcgcc                                                                                                                 >  2:791111/1‑139 (MQ=255)                             aCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTTCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCAt                                                                                                                >  2:946537/1‑140 (MQ=255)                               ccAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCt                                                                                                              >  2:671640/1‑140 (MQ=255)                                  gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                           >  1:187919/1‑140 (MQ=255)                                  gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                           <  1:565289/140‑1 (MQ=255)                                  gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                           >  1:434278/1‑140 (MQ=255)                                  gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                           >  2:1158534/1‑140 (MQ=255)                                  gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTCGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                           >  1:812782/1‑140 (MQ=255)                                   aCGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgcg                                                                                                                   <  1:246891/131‑1 (MQ=255)                                   aCGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAg                                                                                                          <  1:867410/140‑1 (MQ=255)                                   aCGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAgcgcg                                                                                                                   >  2:246891/1‑131 (MQ=255)                                     ggAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCa                                                                                                        <  2:1139955/140‑1 (MQ=255)                                      gAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAg                                                                                                       >  2:63162/1‑140 (MQ=255)                                        ggTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGGAAAAAAAGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCCCAGCGCGCATCTGAAGCAggg                                                                                                     >  1:612570/1‑139 (MQ=255)                                         gTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTTCAGAAAACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAAGGGGGCAGCACAGCGGGGGGAAATTTAAGAAAGACCCCCCGCGCGCATCTGAAGCAggga                                                                                                    >  2:612640/1‑138 (MQ=255)                                         gTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGGAAGAGCATCAGAAGCAGGGAAAAAAGGGGGACAGCACAACGGGGGGAATTTCAAAAATGACCCGCAGCGCGCATCTGAAGCAGGTa                                                                                                    >  2:867410/1‑140 (MQ=255)                                         gTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAggga                                                                                                    >  1:1016532/1‑138 (MQ=255)                                             gTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaaa                                                                                                <  1:1165563/140‑1 (MQ=255)                                              tAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaaaa                                                                                               <  2:39099/140‑1 (MQ=255)                                                  ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggtg                                                                                           <  1:612640/140‑1 (MQ=255)                                                   gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggggg                                                                                          >  1:450218/1‑140 (MQ=255)                                                    cAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggtggt                                                                                         <  1:671640/140‑1 (MQ=255)                                                        ccATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAaga                                                                                     >  1:396195/1‑140 (MQ=255)                                                         cATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAagag                                                                                    <  1:508696/140‑1 (MQ=255)                                                         cATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAagag                                                                                    <  1:1161841/140‑1 (MQ=255)                                                         cATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAagag                                                                                    <  2:434278/140‑1 (MQ=255)                                                             gAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTa                                                                                        <  2:201185/132‑1 (MQ=255)                                                             gAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGt                                                                                <  1:596265/140‑1 (MQ=255)                                                             gAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTa                                                                                        >  1:201185/1‑132 (MQ=255)                                                                gCGGTTCCGGCAATTTTTCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGGGGTAAGAGCAGTCAc                                                                             >  1:863516/1‑140 (MQ=255)                                                                    ttCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCa                                                                         <  1:458278/140‑1 (MQ=255)                                                                         gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGc                                                                    <  2:396195/140‑1 (MQ=255)                                                                                  cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                                           <  2:612570/140‑1 (MQ=255)                                                                                   agaagaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGc                                                          <  1:63162/140‑1 (MQ=255)                                                                                       gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                                      <  2:18275/140‑1 (MQ=255)                                                                                            cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                                 >  2:989752/1‑140 (MQ=255)                                                                                                     aTCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa                                        <  2:187919/140‑1 (MQ=255)                                                                                                       cAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAAt                                      <  2:1130312/140‑1 (MQ=255)                                                                                                       cAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAAt                                      <  1:366211/140‑1 (MQ=255)                                                                                                               ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGACAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                              >  1:662221/1‑140 (MQ=255)                                                                                                                    aaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTg                         <  1:1070920/140‑1 (MQ=255)                                                                                                                          tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCt                    <  2:317726/139‑1 (MQ=255)                                                                                                                          tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCt                    >  1:317726/1‑139 (MQ=255)                                                                                                                          tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc                   <  1:791111/140‑1 (MQ=255)                                                                                                                          tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACAAATGCATTTTTGTCTGTAGCTc                   >  1:871718/1‑140 (MQ=255)                                                                                                                           ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCg                  >  1:1217304/1‑140 (MQ=255)                                                                                                                           ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCg                  >  2:947454/1‑140 (MQ=255)                                                                                                                                      cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCa       >  1:1020757/1‑140 (MQ=255)                                                                                                                                      cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCa       >  2:262987/1‑140 (MQ=255)                                                                                                                                         gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCa    >  2:1072544/1‑140 (MQ=255)                                                                                                                                          ggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAc   <  2:863516/140‑1 (MQ=255)                                                                                                                                           gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc  >  1:512274/1‑140 (MQ=255)                                                                                                                                           gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc  >  1:198582/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                           CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACC  >  NC_000913/1068034‑1068311

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑