Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	44.4%	98.3 / 20.4	63	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/28);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGC  >  NC_000913/1376361‑1376653                                                                                                                                                  |                                                                                                                                                     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                 >  1:322480/1‑150 (MQ=255)  cgccgATGACCTCGCCGCGTAGCCTGGACGGCGTCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                               <  2:352054/151‑1 (MQ=255)    ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                             <  1:836708/151‑1 (MQ=255)      gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                           >  2:434056/1‑151 (MQ=255)          aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                       >  1:398200/1‑151 (MQ=255)            cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                     >  1:187371/1‑151 (MQ=255)             tCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCTCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                    >  1:1015067/1‑151 (MQ=255)               gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                  >  1:765582/1‑151 (MQ=255)                     aagcggggACGGCGACCGAAAAAGGGTGAGTGTTCCCCTGCGCGTACGGAAACAGGTCGCCGCGCGTACGGGACCCGTGGTGTATCTCGACCTGTAAAC‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACa                                                                                                                              <  2:706987/143‑1 (MQ=255)                           ggACGGCGCCCGAAAAATGGTGGGTGATCCCCGGCGCGTACTGGAACGGGCCCCCGCGCGTCCCGGACCGGGGTTGTATCCCAACCCGTAAAC‑CCGCTACCAGGATAAAATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                      <  1:469677/151‑1 (MQ=255)                                gcgcCCGAAAAAGGGTGGGTGATCCCCTGCGCGTACGGGAAAGGGTCGCCGCGCGTACCGGAACCGGGGTGTATCTCGACCCTTAAAA‑CCGTTATAGGGATAATATGATCACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                   <  2:830762/149‑1 (MQ=255)                                    ccGAAAAATGGTGAGTGATCCCCTGCGCGTCCTGGAACAGGTCGCCGGGCGTCCCGGACCGGTGGTGTATCTCGACCCGTAAAC‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                             <  2:394200/151‑1 (MQ=255)                                        aaaaGGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGAACCTTAAAC‑CCGTTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                           <  2:348929/149‑1 (MQ=255)                                        aaaaGGGTGAGTGATCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCGGTGGTGTACCTCGACCTTTTAAC‑CCGCTACAAGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                         <  2:1172272/151‑1 (MQ=255)                                           aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                      >  1:1101178/1‑151 (MQ=255)                                                gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                   >  1:651376/1‑149 (MQ=255)                                                gagtgaTCCCCGGCGCGTACGGAAACAGGTCGCCGCGCGTACCGGAACCGGGGTGTATCTCAAACCTTTAAC‑CCGCTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                  <  2:5260/150‑1 (MQ=255)                                                   tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                               >  1:156544/1‑150 (MQ=255)                                                        cccGGCGCGTACGGGAACAGGCCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCGTTAAC‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                          <  2:322515/150‑1 (MQ=255)                                                         ccggcgcGTCCTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCAACCCGTTAAA‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                        <  2:567163/148‑1 (MQ=255)                                                              gcgTACGGGAACGGGTCGCCGCGCGTCCCGGACCGGGGGTGTATCTCGACCCTTTAAC‑CCTCTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                    <  2:1015170/150‑1 (MQ=255)                                                               cgTACTGGAACGGGTCGCCGCGCGTACCGGAACGGTGGTGTATCCCGACCCTTAAAC‑CCGCTATCAGAATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                       <  2:654745/146‑1 (MQ=255)                                                                 tACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                 >  2:1238229/1‑150 (MQ=255)                                                                   cggaaaCGGGTCGCCGCGTGTCCCGGAACGGTGGTGTATCTCGACCCGTAAAC‑CCGCTATCAGGATCATAGGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                 <  2:1194239/144‑1 (MQ=255)                                                                    tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                              >  2:568385/1‑150 (MQ=255)                                                                     ggagcgggccgccgcGCGTACCGAAACGGTGGTGTATCCGGCACCTTTAAA‑CCGTTAAGAGGATCAAATCATAAACCCCCTCGCGAATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                                <  2:1048095/138‑1 (MQ=255)                                                                     gaaaCAGGCCGCCGCGCGTACCGGACCGGTGGTTTATCCCGAACCTTAAAC‑CCGCTATCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                            <  1:305988/149‑1 (MQ=255)                                                                       aaCGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCCCGAACCGTTAAC‑CCGCTATCAGGATAATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                                <  2:1057785/145‑1 (MQ=255)                                                                       aaCGGGTCCCCGCGCTTACGGGACCGGTGGTTTATCCCGACCCTTAAAC‑CCGTTACAAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                           <  2:493969/150‑1 (MQ=255)                                                                       aaCAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                          >  1:843585/1‑151 (MQ=255)                                                                        aCAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                         >  2:970690/1‑151 (MQ=255)                                                                         cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                         >  2:1165801/1‑150 (MQ=255)                                                                         cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                         >  2:327182/1‑150 (MQ=255)                                                                           ggTCGCCGCGCGGCCGGGACCCGTGGTGTATCCCGAACCGTTAAC‑CCGCTATCAGGTTCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                       <  2:328006/150‑1 (MQ=255)                                                                              cgcggcgatacccGGAACCGTGTTGTACCTCGAACCGTAAAC‑CCGTTATAAGGATAATAGGATCAACACCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                     <  2:710177/138‑1 (MQ=255)                                                                              cgccggGCGCCCGGGACCCGTGTTTTATCCCGACCTGTAAAC‑CCATTATAAGAATCATATGATAACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                     <  1:183181/133‑1 (MQ=255)                                                                              cgccgcGCGTCCCGGAACCGTGGTGTATCTCGACCCTTTAAC‑CCGCTATCAGAATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                   <  2:765649/151‑1 (MQ=255)                                                                                   cgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                 >  1:571170/1‑148 (MQ=255)                                                                                   cccgTCCCGAACCCGTGTTTTATCTCGAACGGTTAAC‑CCGCTATCAGTATCATATGATCACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                               <  2:579280/148‑1 (MQ=255)                                                                                     cgcaCCGAACCCGTGGTTTACCCCGACCCGTAACC‑CCGCTATCAGGTTAAAAAAACCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                             <  1:231134/147‑1 (MQ=255)                                                                                      gTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                              >  1:270831/1‑148 (MQ=255)                                                                                           ggAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                      >  2:867833/1‑151 (MQ=255)                                                                                            gAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                     >  1:302648/1‑151 (MQ=255)                                                                                              aCCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                   >  1:377421/1‑151 (MQ=255)                                                                                                   ggttTATCTCAACCCTTTAAC‑CCCTTATAAGAATAATATGATCAACCCCCTCGCCGTTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                               <  2:675989/150‑1 (MQ=255)                                                                                                   ggtgTACCCCGACCCTTAAAC‑CCGATAAAAGGATCAAATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                               <  1:733825/150‑1 (MQ=255)                                                                                                     tgtATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                            >  2:169275/1‑151 (MQ=255)                                                                                                           ccaacccGTAAAC‑CCGTTAACAGAATAATATAACCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGGGCt                                         <  1:1034133/143‑1 (MQ=255)                                                                                                              aaCCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                    >  2:146759/1‑150 (MQ=255)                                                                                                               aCCCTTAAACCCGCAATCAGGATAATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                           <  2:618034/143‑1 (MQ=255)                                                                                                                    tAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATGGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                             >  1:808959/1‑151 (MQ=255)                                                                                                                     aaa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                            >  2:1249101/1‑151 (MQ=255)                                                                                                                         ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                         >  1:851593/1‑151 (MQ=255)                                                                                                                                caggttaatataacaaaCCCCCCCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                   <  2:398239/135‑1 (MQ=255)                                                                                                                                 agtaTCATATGACCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                        <  2:651434/141‑1 (MQ=255)                                                                                                                                  ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                >  1:319423/1‑151 (MQ=255)                                                                                                                                   gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc               >  1:1252923/1‑151 (MQ=255)                                                                                                                                     tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt             >  1:557805/1‑151 (MQ=255)                                                                                                                                        tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg          >  1:160414/1‑151 (MQ=255)                                                                                                                                        tatGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt             <  2:753723/148‑1 (MQ=255)                                                                                                                                            aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat       >  2:16782/1‑150 (MQ=255)                                                                                                                                             tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat     >  1:137483/1‑151 (MQ=255)                                                                                                                                               aaCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt           <  2:328561/143‑1 (MQ=255)                                                                                                                                                acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc  >  1:1196771/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGC  >  NC_000913/1376361‑1376653

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