Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,751	0	A	C	40.5%	122.7 / 33.3	85	N59K (AAT→AAG)	ralR	Rac prophage; restriction alleviation protein
Reads supporting (aligned to +/- strand):  ref base A (20/30);  new base C (34/0);  total (55/30)
Fisher's exact test for biased strand distribution p-value = 9.50e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGG  >  NC_000913/1413602‑1413900                                                                                                                                                          |                                                                                                                                                       gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCat                                                                                                                                                       <  2:1106679/151‑1 (MQ=255)        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCatt‑attat                                                                                                                                                <  1:623545/150‑1 (MQ=255)        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCatt‑atta                                                                                                                                                 <  2:506406/149‑1 (MQ=255)         gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTAA‑ATTAAAA‑CCT‑CCCTT‑Atttttt                                                                                                                                              >  1:208747/1‑151 (MQ=255)            aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCa                                                                                                                                           <  1:1293552/151‑1 (MQ=255)               ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAT‑ATTAAAA‑CCT‑CCCTT‑TTTATTTTCAATg                                                                                                                                        >  2:990423/1‑151 (MQ=255)                    ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTAA‑ATAAACA‑CAT‑CCCTT‑TTTATTTTCCGTGGTTCg                                                                                                                                   >  2:1380333/1‑151 (MQ=255)                     tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTc                                                                                                                                    <  2:1153064/149‑1 (MQ=255)                      cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGTTCGtt                                                                                                                                 >  2:407511/1‑151 (MQ=255)                        tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGttt                                                                                                                                <  1:719945/150‑1 (MQ=255)                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCC‑CCATT‑ATTATTTTCAGTGGGTCGTTTAt                                                                                                                              >  2:1232987/1‑151 (MQ=255)                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTa                                                                                                                               <  2:1266808/150‑1 (MQ=255)                          tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTAt                                                                                                                              <  2:117135/150‑1 (MQ=255)                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCT‑CCATT‑ATTTTTTTCATTgggttttttttttt                                                                                                                           >  1:1286723/1‑139 (MQ=255)                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTc                                                                                                                            <  1:741955/150‑1 (MQ=255)                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTcc                                                                                                                           <  2:1396639/150‑1 (MQ=255)                              ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATTAACA‑CCC‑CCATT‑TTTTTTTTCCGTGGGTCGttttttttat                                                                                                                         >  1:656525/1‑147 (MQ=255)                                  agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGTTTGGTTTTTTCATCttt                                                                                                                     >  2:174525/1‑151 (MQ=255)                                  agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                                     <  1:284792/151‑1 (MQ=255)                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                                   <  1:929533/151‑1 (MQ=255)                                     gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                                                                                                                  <  1:260012/151‑1 (MQ=255)                                         tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAAAA‑CCC‑CCCCT‑TTTATTTCCCGTGGGTCGTTTTTTcccttctttttttgg                                                                                                              >  2:1037464/1‑138 (MQ=255)                                          ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTTTTTTCAGTGGTTCGTTTTTTTCATCTTTTTAGTGct                                                                                                             >  1:1167906/1‑151 (MQ=255)                                           tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAAAA‑ACT‑CCATT‑TTTTTTTTCAGTGGTTCGTTTTTTTCATCTTTTTAGTGctt                                                                                                            >  1:630877/1‑151 (MQ=255)                                            ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATAAAAA‑ACC‑CCCTT‑ATTATTTTCCGTGGTTCGTTTATTTCCTCTTTCGAGTGcttc                                                                                                           >  2:274845/1‑151 (MQ=255)                                            ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAT‑AATAAAA‑CCT‑CCCTT‑ATTATTTTTTGTGGGTCGGTTTTTCCTTCTTTCggattcttc                                                                                                           >  2:494305/1‑143 (MQ=255)                                              cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑ACC‑CCCTT‑TTTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGcttttt                                                                                                         >  1:578789/1‑151 (MQ=255)                                              cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑AACAAAA‑CCTCCCATA‑TTTTTTTCCAAGGGTT‑GTTTTTTTCAACTTTctgggggttttt                                                                                                         >  2:128922/1‑140 (MQ=255)                                               gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                                        <  2:173677/151‑1 (MQ=255)                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAAAA‑CCT‑CCCTT‑ATTATTTTCCGTGGGTCGTTTTTTTCCTCTTTTTAGTGCtttttttt                                                                                                      >  2:182064/1‑151 (MQ=255)                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                                      <  1:159647/151‑1 (MQ=255)                                                   tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TCATTTTTTTA‑ATGAAAA‑CCC‑CCCCT‑TTTTTTTTTTGTGGTTCGGTTATTTCCTCTTTTGAGTGCCTCttttttt                                                                                                    >  1:1120081/1‑149 (MQ=255)                                                    cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTAATTTTTAA‑ATGAACA‑CCT‑TCATT‑ATTATTTTCAGGGGTTTTTTTTTTTCACCTTTTCCGGGTTTCTTTTTCAc                                                                                                   >  2:660427/1‑151 (MQ=255)                                                     tGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑AAGAAAA‑ACC‑CCCTT‑TTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTGAGTGGCTCttttttaat                                                                                                  >  1:1226529/1‑147 (MQ=255)                                                      gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTTC‑ATGAACA‑CCT‑CCCTT‑TTTTTTTTCACAGGTGTGTTTTTTTCATCTTTCCAGTTCTTCTTTTTTACtt                                                                                                 >  1:123877/1‑151 (MQ=255)                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATTAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGTTCGTTTTTTTCCTCTTTCTAGGGCTTTTTTTTTACTTccacc                                                                                            >  2:1369874/1‑151 (MQ=255)                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATTT‑AATAACA‑ACC‑CCCTT‑ATTATTTCCAGTGGGTCGTTTTTTTCTTCTTTTCGGGGGTTCTTTTTTACTTccacc                                                                                            >  2:1211761/1‑151 (MQ=255)                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCC‑CCCTT‑ATTATTTTCCGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTTTTTTTTACTTccccc                                                                                            >  2:751608/1‑151 (MQ=255)                                                            gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAACA‑ACT‑CCCTT‑TTTTTTTTCAGTGGTTTGTTTTTTTCCTTTTTTGGGTGCCTCTTTTTTTCTTtcacc                                                                                            >  1:1359334/1‑150 (MQ=255)                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGATAATTTT‑TAA‑AATAACA‑CCT‑CCAAT‑TTTTTTTTCAGGGGTTGGTGTTTTATATTCTTTGTAAGCTTCTTTTTCCCTTcctccaca                                                                                         >  2:317675/1‑148 (MQ=255)                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTATTTTCAGTGGGTTGTTTATTTCCTCTTTTGAGTGGTTTTTTTTTTCTTCCACCATa                                                                                         >  2:1167004/1‑150 (MQ=255)                                                                 tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑GAAATTTTTAC‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCACTGGTTCGTTTTTTTCCTCTTTTGAGTGGTTTTTTTTTACTTTCACCCTaaac                                                                                      >  2:773314/1‑149 (MQ=255)                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAACA‑CCC‑CCCCT‑ATTTTTTCCAGCGGGTTGTTTATTCCACCTTTTCAGTGCTTCTTTTTTACTTCCCCCAtcaacc                                                                                     >  2:182430/1‑146 (MQ=255)                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc                                                                                      <  1:848891/150‑1 (MQ=255)                                                                    aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCgg                                                                                    <  1:1098802/150‑1 (MQ=255)                                                                    aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGt                                                                                   <  1:1118772/151‑1 (MQ=255)                                                                     gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAAAA‑CCC‑CCCAT‑ATTATTTTCAATGGGTTGTTTTTTTCATTTTTTCAGGGCTTTTTTTTTTCTTCCCCCCCAAcccggt                                                                                  >  1:1203874/1‑147 (MQ=255)                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTTTCTTTCCCCCTAAACGGGTCgggg                                                                             >  1:454672/1‑150 (MQ=255)                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAACA‑CCC‑CCATT‑ATTATTTTCAATGGTTTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTACTTTCCCCATAAACGGGTTgggg                                                                             >  2:814336/1‑150 (MQ=255)                                                                            aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑AATAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCACCCTAACCGGGTCGGGGtt                                                                           >  2:1001154/1‑151 (MQ=255)                                                                            aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTATTA‑ATGAACA‑CCC‑CCCTT‑TTTATTTTCAGTGGTTTGTTTTTTTCCTCTTTTGAGGGCTTCTTTTTTACTTTCACCATAACCGGTTTgggggt                                                                           >  1:332997/1‑148 (MQ=255)                                                                               atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑AAGAACACCCC‑CCTTT‑TTTTTTTCCCGGGGTCCG‑TTATTCCATTTTTTGAGGGGTTCTTTTTTACTTCCACCATAAACGGCTCGGGGTTTGc                                                                        >  2:1419609/1‑151 (MQ=255)                                                                                 cTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCC‑CCATT‑TTTATTTTCAGTGGTTTGTTTTTTTCATCTTTTGAGTGCTTTTTTTTTTCTTCCACCCTAAACCGGTCGGGGTTCGCAt                                                                      >  2:218038/1‑151 (MQ=255)                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAAAA‑ACT‑CCCTTAATTTTTTTCAG‑GGTGCGTTTTTTTCTTCTTTTTAGTGCTTCTTTTTTACTTTCACCATCAACGCGTTGGGGTTCGCCTcc                                                                    >  1:18429/1‑151 (MQ=255)                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTTTTTCCCGTGGGTTGTTTTTTTCCTCTTTTTAGTGGTTCTTTTTTTCTTTCCCCATAACCCGGTTGGGGTTTGCAttc                                                                    >  1:72475/1‑149 (MQ=255)                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGc                                                                        <  2:1340501/147‑1 (MQ=255)                                                                                    gTCCAGTNAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTANNTNNNGTGGTTCGTTTATTNNNNNNNTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCg                                                                   <  2:105962/151‑1 (MQ=255)                                                                                      ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGGTTGTTTTTTTCATTTTTTTAGTGCTTCTTTTTTTCTTTCACCATAACCGGTTTGGGGTTCGCCTCCGtt                                                                 >  1:751654/1‑151 (MQ=255)                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTTACTTCCCCCCTAAACGGGTCGGGGTTTCCATCCCTTac                                                               >  2:1411708/1‑150 (MQ=255)                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAT‑AAGGAAA‑CCC‑CCCTT‑TTTTTTTTCCCTGGTTTGTTTTTTTCCTCTTTTTTGGGCTTCTTTTTTTCTTCCCCCCTAAACCGTTCGGGGTTTGCCTCCCtttcc                                                              >  1:736404/1‑148 (MQ=255)                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTaca                                                              <  1:1410352/151‑1 (MQ=255)                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                          <  1:925649/151‑1 (MQ=255)                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                          <  1:209045/151‑1 (MQ=255)                                                                                                  aGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCg                                                     <  1:407519/151‑1 (MQ=255)                                                                                                   gTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCCCCCTAAACGGGTTGGGGTTCGCATCCCGTACACCTCGCTCgg                                                    >  2:51893/1‑151 (MQ=255)                                                                                                     gCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAACA‑CCC‑CCCCT‑TTTTTTTTCAGTGGTTTGTTTTTTTCATCTTTCGCGGGCTTCTTTTTTTCTTCCACCATAACCGGTTTCGGGTTCGCCTCCGTTACACCTCGCTCggg                                                   >  2:700145/1‑149 (MQ=255)                                                                                                      ccACCTTCCTTTTCAATAGTGGCGG‑GAATTTTTTACAATTAACC‑CCC‑CCCTT‑TTTTTTTCCA‑TGGTTCGTTTATTTCTTCTTTCGTGTGCTTCTTTTTTTCTTTCTCCATTACCGGGTTGGGGTTGTCTTCCGTTACACTTCGCTCGTTa                                                  >  1:182619/1‑150 (MQ=255)                                                                                                             ccttTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTCGTTTATTTCATCTTTTGAGTGCTTCTTTTTTACTTCCCCCATAAACGGTTTGGGGTTTGCATCCGTTACACTTCGCTCGGGAATATCCTg                                          >  2:303610/1‑151 (MQ=255)                                                                                                                   cAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATgg                                    <  1:1186173/151‑1 (MQ=255)                                                                                                                       aGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                                 <  2:124788/150‑1 (MQ=255)                                                                                                                        gTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                               <  2:1008895/151‑1 (MQ=255)                                                                                                                        gTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGGTCGTTTATTTCATCTTTTGAGTGCTTCTTTTTTACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtttt                               >  1:134264/1‑150 (MQ=255)                                                                                                                           gCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc                            >  2:1381507/1‑151 (MQ=255)                                                                                                                              g‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt                         >  1:1043222/1‑151 (MQ=255)                                                                                                                                 aaTTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc                       >  1:476681/1‑151 (MQ=255)                                                                                                                                     ttATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt                   >  1:386179/1‑151 (MQ=255)                                                                                                                                     ttATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt                   >  1:1373104/1‑151 (MQ=255)                                                                                                                                      tataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt                   <  1:139761/150‑1 (MQ=255)                                                                                                                                      tataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg                  >  1:323660/1‑151 (MQ=255)                                                                                                                                       ataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg                 <  1:1437566/151‑1 (MQ=255)                                                                                                                                       ataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg                 <  2:1202449/151‑1 (MQ=255)                                                                                                                                         aC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGa                <  2:881938/150‑1 (MQ=255)                                                                                                                                               aaCA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAcc           >  1:963378/1‑150 (MQ=255)                                                                                                                                                aca‑cCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACcac         <  1:1037460/151‑1 (MQ=255)                                                                                                                                                 ca‑cCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACcaca        <  1:494312/151‑1 (MQ=255)                                                                                                                                                     cT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACaaa      >  2:748175/1‑150 (MQ=255)                                                                                                                                                        ccATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATg   >  2:165476/1‑151 (MQ=255)                                                                                                                                                         cATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATgg  >  1:941855/1‑151 (MQ=255)                                                                                                                                                          att‑attattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATgg  >  1:48290/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                       GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGG  >  NC_000913/1413602‑1413900

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑