Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,254	A→T	9.7%	intergenic (+32/‑115)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,254	0	A	T	9.7%	171.1 / 9.2	66	intergenic (+32/‑115)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base A (39/20);  new base T (2/5);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 9.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCA  >  NC_000913/781112‑781394                                                                                                                                               |                                                                                                                                               ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                                         >  2:194935/1‑150 (MQ=255)   tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTATCACAGTTGGTATAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAAAACAGCACGACCCACCACTTTTAAAATGTAAAAAAGAAGATAACCGTGAAGGTTAACg                                                                                                                                       >  1:573292/1‑150 (MQ=11)   tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCACTTTTAAAATCAAAATCAGATGTTAAGCGTGAAGGATAACg                                                                                                                                       >  1:557466/1‑150 (MQ=255)   tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTGAAGCATAACg                                                                                                                                       >  1:1287265/1‑150 (MQ=255)   tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAAATGGTAGCAGGTTCGAATCCTGCACGACCCACCAATTTTAACATCAAAATAAGATTTTAAACGTGAAGGATAACg                                                                                                                                       >  1:557472/1‑150 (MQ=14)     ggcgtTTCTATCCCGTCATCACCAACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                                    <  2:305957/150‑1 (MQ=255)     ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTCTAGCTCAGATGTTAATCGTGAAGGTTAACGtt                                                                                                                                     >  2:657010/1‑150 (MQ=14)        ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                                  <  2:1039742/150‑1 (MQ=255)         gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcgc                                                                                                                                 <  2:351395/150‑1 (MQ=255)         gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcgc                                                                                                                                 >  1:905998/1‑150 (MQ=255)          ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTGAAGCGTGAAGGATAACGTTgcgac                                                                                                                                >  2:918733/1‑148 (MQ=255)             gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                             >  1:1026542/1‑150 (MQ=255)               tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCaa                                                                                                                           >  1:219707/1‑150 (MQ=255)               tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCaa                                                                                                                           >  1:219649/1‑150 (MQ=255)                   gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGc                                                                                                                       <  2:832829/150‑1 (MQ=255)                   gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCCCCACTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAAAGGTGCGCCAGCAACCGc                                                                                                                       >  1:603896/1‑150 (MQ=14)                            caccacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATCAATTTGTCGCAGGTTCGAATACTGCACGACCCACCACATTTAAAATAAAAAAAAGTAGTTAAGCGGTAAGGTTAACGTTTCGCAAGCAAACGCCCggatcgc                                                                                                              >  1:1172403/1‑144 (MQ=9)                                       tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAACCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag‑cga                                                                                                 <  2:1178726/150‑1 (MQ=255)                                             gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTc                                                                                              <  2:561637/147‑1 (MQ=255)                                             gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTc                                                                                              >  1:561637/1‑147 (MQ=37)                                               tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCt                                                                                         >  1:210453/1‑150 (MQ=255)                                                         gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                              >  1:222084/1‑89 (MQ=255)                                                         gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                              <  2:222084/89‑1 (MQ=255)                                                            cAGTTTACTTTTTATCAATTGGGCGCAGGTTCGAATCCCGGACGACCCCCCCGTTTTAACATCCAACTCAGATGTTTAGCGTGAAGGATAACGTTGCGTTAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACccac                                                                              <  2:1007022/150‑1 (MQ=11)                                                              cttGTCTTTTTATAAATTGGTCGCATGTTCGAATCCGGCACGTACCAACAGTTTTAACTTCACACTCAGTTGTTAATCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACccacca                                                                            <  2:557466/149‑1 (MQ=14)                                                                   cTTTGAAACAAGTGGGCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCTTTATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTaa                                                                     <  1:997388/150‑1 (MQ=14)                                                                     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                                   >  1:976413/1‑150 (MQ=255)                                                                       tAATCAATTGGTCGCAGGTTCGAATCGTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                                 <  2:1296420/150‑1 (MQ=255)                                                                        aaTCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                                <  2:1352754/150‑1 (MQ=255)                                                                         aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                               >  2:645642/1‑150 (MQ=255)                                                                                ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                        <  2:402955/150‑1 (MQ=11)                                                                                    gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTATGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCCCGACCCACCAATCTTAAAGAGTGGCCCCGAGTaa                                                    >  1:1207476/1‑150 (MQ=255)                                                                                    gCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTAGCGCCAGCAACCGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAGGGGCCCCGAGTat                                                    >  1:1207481/1‑149 (MQ=17)                                                                                      aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                  >  1:1235245/1‑150 (MQ=255)                                                                                          tCGACTCCTGCACGACCCACCAGTTTGAACATCTAACTCAGAGGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                              <  2:780447/150‑1 (MQ=18)                                                                                          tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                              <  2:536353/150‑1 (MQ=255)                                                                                          tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCAAGAACGACACACCAATCTGAGAGAGTGGCCCCGAGGAAAAAACt                                              >  1:1254157/1‑150 (MQ=17)                                                                                          tCGAATCCTGCACGACCCACCAGTTTTAACAACAAACTCAGATGATAACCGAGAAGGATACCTATGAGCAAGCAACGGCCCGTAGGCAGAG‑CGAAG‑CGAGTCATCCTGATCGACCCACCAATCTACAAGATTGGCCCCGAGTagacctct                                              >  1:1254044/1‑143 (MQ=17)                                                                                              aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                          <  1:1403671/150‑1 (MQ=255)                                                                                               tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                         >  1:733336/1‑150 (MQ=255)                                                                                               tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                         <  2:855446/150‑1 (MQ=255)                                                                                                 cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                       >  1:224533/1‑150 (MQ=255)                                                                                                  tGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                      >  1:48242/1‑150 (MQ=17)                                                                                                  tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                      >  2:1219647/1‑150 (MQ=255)                                                                                                   gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                     <  1:288118/150‑1 (MQ=255)                                                                                                    cACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                                    >  1:1227426/1‑150 (MQ=255)                                                                                                     ccGACCCACCAGTTTTAACATCATACACAGTTGATAATCGTGATGGATAACTCTGCGCCAGAACAGGCCCGTATGGCGAG‑CTAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAaca                                   <  2:186767/149‑1 (MQ=21)                                                                                                     aCGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAaca                                   <  2:39430/150‑1 (MQ=255)                                                                                                      cGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                                  <  2:48242/150‑1 (MQ=17)                                                                                                             ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                           >  1:1225429/1‑150 (MQ=255)                                                                                                             ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                           >  1:1225408/1‑150 (MQ=255)                                                                                                             ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGATg                           >  1:1225432/1‑150 (MQ=255)                                                                                                                   ataccATGAACCTCAGATGTGAAGCGTGTAGGATCACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAACGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                     <  2:1391280/146‑1 (MQ=21)                                                                                                                        aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCATCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCACGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                <  2:905998/150‑1 (MQ=255)                                                                                                                        aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                >  1:111394/1‑150 (MQ=255)                                                                                                                         tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc               >  2:430867/1‑150 (MQ=255)                                                                                                                               tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg         >  1:1072673/1‑150 (MQ=255)                                                                                                                                 aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGTCCCCGAGTAAAAAACTTTCAGGTAAAAACCGAATGGGTAGTTAGATCAGTTGGTa       >  2:758054/1‑150 (MQ=21)                                                                                                                                 aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGtc       >  1:1390977/1‑149 (MQ=255)                                                                                                                                 aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa       >  1:1391062/1‑150 (MQ=255)                                                                                                                                  gATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                                 <  2:996734/91‑1 (MQ=255)                                                                                                                                  gATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                                 >  1:996734/1‑91 (MQ=255)                                                                                                                                     gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc   >  1:85365/1‑150 (MQ=255)                                                                                                                                     gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGATCAGTTGGTAGAGc   >  1:150889/1‑150 (MQ=255)                                                                                                                                     gTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCAACAACCACGACACACCAATCTTAAAGATTGGCCCCGAGTAAAAAACTTGCAGGTAACACCCGAATGGGTCGTAAGCGCAGTTGGTAGAGc   >  1:85429/1‑150 (MQ=25)                                                                                                                                      ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa  <  2:697047/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               GGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCA  >  NC_000913/781112‑781394

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