Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	248,434	G→A	5.2%	pseudogene (1637/1713 nt)	lfhA ←	flagellar system protein, promoterless fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	248,434	0	G	A	5.2%	244.2 / 2.1	76	pseudogene (1637/1713 nt)	lfhA	flagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base A (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCGGCAC  >  NC_000913/248292‑248579                                                                                                                                               |                                                                                                                                                   gCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAaacaa                                                                                                                                             >  2:915320/1‑150 (MQ=255)   cTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAAcc                                                                                                                                           >  2:1067810/1‑150 (MQ=255)      aaCGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTg                                                                                                                                        <  1:1230722/150‑1 (MQ=255)       aCGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGc                                                                                                                                       <  1:667322/150‑1 (MQ=255)       aCGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGc                                                                                                                                       <  1:877491/150‑1 (MQ=255)           tAGGACGGATAAGGCGTTCACGACGCATACGGCAGTTGTGAACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAACTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCACAAACAACCTTGCATaa                                                                                                                                   >  2:313843/1‑150 (MQ=255)             ggACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAAcg                                                                                                                                 >  2:1369305/1‑150 (MQ=255)              gACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAAcgc                                                                                                                                <  2:327409/150‑1 (MQ=255)              gACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAAcgc                                                                                                                                >  2:1146295/1‑150 (MQ=255)              gACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAAcgc                                                                                                                                <  2:639217/150‑1 (MQ=255)                        gCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAg                                                                                                                      <  1:469295/150‑1 (MQ=255)                           ttCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGgc                                                                                                                   >  1:1297762/1‑150 (MQ=255)                             cACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGgcgc                                                                                                                 >  1:330977/1‑150 (MQ=255)                               cgccgcATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAg                                                                                                               >  2:235309/1‑150 (MQ=255)                                gccgcATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGc                                                                                                              <  1:980113/150‑1 (MQ=255)                                       ccGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACCTGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTg                                                                                                      >  2:30519/1‑150 (MQ=255)                                       ccGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTg                                                                                                      <  2:704358/150‑1 (MQ=255)                                          gCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCa                                                                                                   <  1:465456/150‑1 (MQ=255)                                          gCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCa                                                                                                   <  1:1471003/150‑1 (MQ=255)                                          gCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCa                                                                                                   <  1:1219517/150‑1 (MQ=255)                                           cAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCAc                                                                                                  <  2:69833/150‑1 (MQ=255)                                                gtgCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCACCAAAGGGCGCAGCTGC‑GGTGGCACcagca                                                                                             >  2:1402980/1‑150 (MQ=255)                                                  gCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCagcagc                                                                                           >  2:1035056/1‑150 (MQ=255)                                                   cACCAATCCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCTGGTGCCACc‑gcagct                                                                                          <  1:437056/150‑2 (MQ=255)                                                           cTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCt                                                                                  >  1:278487/1‑150 (MQ=255)                                                             ggATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCttt                                                                                >  1:12716/1‑150 (MQ=255)                                                                  cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCccg                                                                           >  2:1036896/1‑150 (MQ=255)                                                                  cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCccg                                                                           <  2:700373/150‑1 (MQ=255)                                                                  cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCccg                                                                           >  1:1098911/1‑150 (MQ=255)                                                                  cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCACAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCccg                                                                           >  2:195347/1‑150 (MQ=255)                                                                         cAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCttt                                                                    >  2:1233298/1‑150 (MQ=255)                                                                           aCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCa                                                                  <  1:1369305/150‑1 (MQ=255)                                                                              cccccATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCt                                                               >  1:1179288/1‑150 (MQ=255)                                                                                cccATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGc                                                             <  1:470947/150‑1 (MQ=255)                                                                                cccATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGc                                                             <  1:1402980/150‑1 (MQ=255)                                                                                 ccATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGct                                                            >  1:1157832/1‑150 (MQ=255)                                                                                   aTAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGctct                                                          <  1:1030113/150‑1 (MQ=255)                                                                                       tCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCa                                                      <  1:1067810/150‑1 (MQ=255)                                                                                       tCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCa                                                      <  1:778489/150‑1 (MQ=255)                                                                                        cTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCAc                                                     >  2:118538/1‑150 (MQ=255)                                                                                         ttCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGTATTTGGTGAAGTTGT‑GGTGGGACCAGCTTGAGCGTGTCTTTCTAAGCCGCTTTCATTTGCTGTTTCTcc                                                    <  2:1028726/150‑1 (MQ=255)                                                                                         ttCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCAcc                                                    >  1:264258/1‑150 (MQ=255)                                                                                         ttCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCAcc                                                    <  2:1297762/150‑1 (MQ=255)                                                                                          tCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCt                                                   >  2:298741/1‑150 (MQ=255)                                                                                            aaCTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCtgt                                                 <  2:73816/150‑1 (MQ=255)                                                                                             aCTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCtgtg                                                <  1:235309/150‑1 (MQ=255)                                                                                             aCTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCtgtg                                                <  1:609798/150‑1 (MQ=255)                                                                                             aCTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCtgtg                                                <  1:1146295/150‑1 (MQ=255)                                                                                              cTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTgg                                               <  2:318405/150‑1 (MQ=255)                                                                                                ccAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCa                                             <  1:1155813/150‑1 (MQ=255)                                                                                                   gCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCg                                          <  1:895077/150‑1 (MQ=255)                                                                                                        tCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTActc                                     >  1:377897/1‑150 (MQ=255)                                                                                                         ccGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTActct                                    >  2:540808/1‑150 (MQ=255)                                                                                                               cTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACt                              >  1:1347459/1‑150 (MQ=255)                                                                                                                  cGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt                           >  1:1161623/1‑150 (MQ=255)                                                                                                                      ataCGACAGCACATGCAGCCCCGGCACAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAg                       <  1:313843/150‑1 (MQ=255)                                                                                                                       taCGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGc                      >  2:1321438/1‑150 (MQ=255)                                                                                                                       taCGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGc                      >  1:1345226/1‑150 (MQ=255)                                                                                                                       taCGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGc                      <  2:278487/150‑1 (MQ=255)                                                                                                                        aCGACAGCACATGCAGCCCCGGCACAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCa                     <  1:195347/150‑1 (MQ=255)                                                                                                                         cGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCat                    >  2:1016441/1‑150 (MQ=255)                                                                                                                          gACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCata                   >  1:905547/1‑150 (MQ=255)                                                                                                                             aGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTc                <  2:1347459/150‑1 (MQ=255)                                                                                                                             aGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTc                >  2:621319/1‑150 (MQ=255)                                                                                                                              gCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCg               <  1:915320/150‑1 (MQ=255)                                                                                                                              gCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCg               <  2:1199350/150‑1 (MQ=255)                                                                                                                                 caTGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGAt            >  2:106140/1‑150 (MQ=255)                                                                                                                                  aTGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATc           <  2:1482962/150‑1 (MQ=255)                                                                                                                                  aTGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATc           <  1:298741/150‑1 (MQ=255)                                                                                                                                   tGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCg          <  1:1053833/150‑1 (MQ=255)                                                                                                                                   tGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGCGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCAGATTCGGATCg          >  2:1263132/1‑150 (MQ=255)                                                                                                                                       gCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCg      >  1:812864/1‑150 (MQ=255)                                                                                                                                       gCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCg      >  1:433541/1‑150 (MQ=255)                                                                                                                                       gCCCCGGCGCAAACAACCTTGCATAACGCGCAAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCCACTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTAGGATCGACCg      >  2:145414/1‑150 (MQ=255)                                                                                                                                          ccGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCGGCa   >  2:1024474/1‑150 (MQ=255)                                                                                                                                           cGGCGCAAACAACCTTGCATAACTCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTTGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCGGCAc  <  2:1389275/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGC‑GGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATATTCGGATCGACCGGCAC  >  NC_000913/248292‑248579

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