Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	973,599	C→A	T21T (ACC→ACA)	cmoM →	tRNA cmo(5)U34 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	973,599	0	C	A	98.8%	301.9 / ‑3.2	86	T21T (ACC→ACA)	cmoM	tRNA cmo(5)U34 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (41/44);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

TTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCA  >  NC_000913/973464‑973736                                                                                                                                        |                                                                                                                                           ttCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAG‑TTTCCCGTAACATTTACGGCACCACAAAAgg                                                                                                                                       <  1:816499/140‑1 (MQ=255)          gCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTc                                                                                                                               >  2:951933/1‑140 (MQ=255)           ccGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCg                                                                                                                              <  1:444522/140‑1 (MQ=255)             gTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGAc                                                                                                                            <  1:1420313/140‑1 (MQ=255)              tGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACa                                                                                                                           <  2:1013205/140‑1 (MQ=255)              tGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACa                                                                                                                           <  1:1544044/140‑1 (MQ=255)                gCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCTAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAgg                                                                                                                         <  2:1488959/140‑1 (MQ=255)                gCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAgg                                                                                                                         >  1:1736190/1‑140 (MQ=255)                 cTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGc                                                                                                                        <  2:1736190/140‑1 (MQ=255)                        aaaaTAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCtg                                                                                                                 <  1:363950/140‑1 (MQ=255)                          aaTAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCtgtg                                                                                                               <  2:150883/140‑1 (MQ=255)                          aaTAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCtgtg                                                                                                               >  1:1258260/1‑140 (MQ=255)                           aTAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCA‑CAAAAGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  1:1764779/140‑1 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  2:180359/140‑1 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             >  1:18751/1‑140 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  2:1658667/140‑1 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  2:244904/140‑1 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  2:1258260/140‑1 (MQ=255)                            tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATGTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGc                                                                                                             <  1:1006286/140‑1 (MQ=255)                               cGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑gg                                                                                                         >  1:482404/1‑140 (MQ=255)                               cGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑gg                                                                                                         >  2:62312/1‑140 (MQ=255)                                 tCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGAt                                                                                                       >  1:1376303/1‑140 (MQ=255)                                     tATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCg                                                                                                   >  2:1042174/1‑140 (MQ=255)                                      aTCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGa                                                                                                  >  2:344687/1‑140 (MQ=255)                                              aCGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGc                                                                                          <  2:320814/140‑1 (MQ=255)                                                  cTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATTGCGTGCTGGc                                                                                      >  2:1699075/1‑140 (MQ=255)                                                  cTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGc                                                                                      <  2:817092/140‑1 (MQ=255)                                                     aGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGa                                                                                   >  2:1835588/1‑140 (MQ=255)                                                       ccAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGaaa                                                                                 <  1:72368/140‑1 (MQ=255)                                                          aCCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATgg                                                                              >  2:1066975/1‑140 (MQ=255)                                                          aCCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCATCGACAGGCTATTCTGTGGCA‑GGACCTCGGTCGCGGGCTTGCGGAAAttg                                                                              >  2:249593/1‑138 (MQ=255)                                                           ccGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATggg                                                                             >  2:517291/1‑140 (MQ=255)                                                           ccGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATggg                                                                             >  1:683923/1‑140 (MQ=255)                                                           ccGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATggg                                                                             >  1:1259213/1‑140 (MQ=255)                                                           ccGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATggg                                                                             <  1:344687/140‑1 (MQ=255)                                                             ggTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGcc                                                                           <  2:935873/140‑1 (MQ=255)                                                               tGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGGCACTCGATCGCGTGCTGGCGGAAAGGGGccc                                                                          >  2:713342/1‑140 (MQ=255)                                                                       gCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACt                                                                 <  1:1847286/140‑1 (MQ=255)                                                                       gCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACt                                                                 >  2:1366183/1‑140 (MQ=255)                                                                       gCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACt                                                                 >  1:215709/1‑140 (MQ=255)                                                                         aaTGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGc                                                               >  2:975120/1‑140 (MQ=255)                                                                         aaTGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGc                                                               >  1:58095/1‑140 (MQ=255)                                                                           tGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCgt                                                             >  1:1201663/1‑140 (MQ=255)                                                                           tGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCgt                                                             <  2:18751/140‑1 (MQ=255)                                                                              aGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGtgtg                                                          <  1:975120/140‑1 (MQ=255)                                                                                gATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCt                                                        >  1:994884/1‑140 (MQ=255)                                                                                gATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCt                                                        >  1:584461/1‑140 (MQ=255)                                                                                 aTCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTg                                                       <  2:683923/140‑1 (MQ=255)                                                                                     cAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATg                                                   >  2:846118/1‑140 (MQ=255)                                                                                     cAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATg                                                   >  2:757988/1‑140 (MQ=255)                                                                                     cAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATg                                                   >  2:17512/1‑140 (MQ=255)                                                                                       aTTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCt                                                 <  1:17512/140‑1 (MQ=255)                                                                                       aTTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCt                                                 >  1:1166848/1‑140 (MQ=255)                                                                                       aTTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCt                                                 <  2:58095/140‑1 (MQ=255)                                                                                        ttttGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTg                                                >  1:1249321/1‑140 (MQ=255)                                                                                         tttGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTgg                                               <  2:215709/140‑1 (MQ=255)                                                                                         tttGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTgg                                               >  2:1006955/1‑140 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAGCATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCgg                                            >  2:545470/1‑140 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCtggtgg                                            <  1:356826/140‑1 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCgg                                            <  1:3904/140‑1 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCgg                                            <  2:353606/140‑1 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCgg                                            <  2:1904750/140‑1 (MQ=255)                                                                                            gatgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCgg                                            <  1:1066975/140‑1 (MQ=255)                                                                                              tgatATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTg                                          <  1:517291/140‑1 (MQ=255)                                                                                                 tatTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgag                                       <  2:1201663/140‑1 (MQ=255)                                                                                                 tatTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgag                                       <  1:1294034/140‑1 (MQ=255)                                                                                                  atTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGTTGGCGGTGgaga                                      >  2:1622493/1‑140 (MQ=255)                                                                                                  atTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgaga                                      >  1:1384652/1‑140 (MQ=255)                                                                                                  atTGCGGAAAAGCTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGTCA‑GGATCTCGATCGCGTGCTGGAGGAAAATGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgaga                                      >  2:1152089/1‑140 (MQ=255)                                                                                                    tGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAg                                    >  2:1788849/1‑140 (MQ=255)                                                                                                      cGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAggg                                  >  2:1219004/1‑140 (MQ=255)                                                                                                      cGGAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAggg                                  >  2:1794479/1‑140 (MQ=255)                                                                                                       ggAAAAGTTTTCCCGTAGCATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGc                                 <  1:545470/140‑1 (MQ=255)                                                                                                       ggAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgaga                                      <  1:567330/135‑1 (MQ=255)                                                                                                       ggAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGgaga                                      >  2:567330/1‑135 (MQ=255)                                                                                                       ggAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGc                                 <  2:1588087/140‑1 (MQ=255)                                                                                                        gAAAAGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCa                                <  1:62312/140‑1 (MQ=255)                                                                                                            aGTTTTCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGAcc                            <  2:969660/140‑1 (MQ=255)                                                                                                               tttCCCGTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCa                         <  1:380782/140‑1 (MQ=255)                                                                                                                     gTAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATTGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaaa                   <  1:1699075/140‑1 (MQ=255)                                                                                                                      tAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCg                           >  1:1353788/1‑131 (MQ=255)                                                                                                                      tAACATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCg                           <  2:1353788/131‑1 (MQ=255)                                                                                                                       aaCATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAAACCGCAATCAAAAt                 >  2:327241/1‑140 (MQ=255)                                                                                                                         cATTTACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATgg               <  2:834699/140‑1 (MQ=255)                                                                                                                            ttACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCg            >  2:1586066/1‑140 (MQ=255)                                                                                                                             tACGGCACCACAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCCGGCCGCAATCAAAATGGCCgg           >  2:1002758/1‑139 (MQ=255)                                                                                                                                  caccacAAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTg      <  2:1793759/140‑1 (MQ=255)                                                                                                                                   accacaAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTgg     >  1:1938080/1‑140 (MQ=255)                                                                                                                                   accacaAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTgg     >  2:449739/1‑140 (MQ=255)                                                                                                                                   accacaAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTgg     >  1:1088355/1‑140 (MQ=255)                                                                                                                                   accacaAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GG‑TCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGGGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  1:1515883/1‑140 (MQ=255)                                                                                                                                    ccacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  2:856025/1‑140 (MQ=255)                                                                                                                                    ccacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  2:36201/1‑140 (MQ=255)                                                                                                                                    ccacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  1:56551/1‑140 (MQ=255)                                                                                                                                    ccacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  1:929203/1‑140 (MQ=255)                                                                                                                                    ccacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg    >  1:101232/1‑140 (MQ=255)                                                                                                                                     cacaaAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GG‑TCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTCGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCa  <  2:1515883/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                           TTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCA‑GGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCA  >  NC_000913/973464‑973736

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