Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	662,636	T→G	T3P (ACC→CCC)	ybeD ←	UPF0250 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	662,636	0	T	G	97.9%	342.8 / ‑1.9	95	T3P (ACC→CCC)	ybeD	UPF0250 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (45/48);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATT  >  NC_000913/662539‑662730                                                                                                  |                                                                                               aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgtt                                                                                              <  1:1057179/100‑1 (MQ=255)  ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgttt                                                                                             >  2:462752/1‑100 (MQ=255)  ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgttt                                                                                             <  1:1481757/100‑1 (MQ=255)   cACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTc                                                                                            <  1:690005/100‑1 (MQ=255)    aCTTCAACCACCTGATCAACCAGCTCAGGTCACACCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGGTTGGGTTTCa                                                                                           >  1:1279880/1‑100 (MQ=255)    aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                           <  2:189470/100‑1 (MQ=255)    aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                           >  2:1593817/1‑100 (MQ=255)     cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCAt                                                                                          >  2:1063283/1‑100 (MQ=255)      ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATg                                                                                         <  1:1688094/100‑1 (MQ=255)      ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATg                                                                                         <  1:148605/100‑1 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAg                                                                                     <  2:117260/100‑1 (MQ=255)           ccaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGc                                                                                    >  2:392199/1‑100 (MQ=255)           ccaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGc                                                                                    <  1:703135/100‑1 (MQ=255)            caccTGATCAACCAGCTCAGGTAACGCCTGCCCGATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTAGTTAAGTTTGGGTTTCATGt                                                                                        >  2:898273/1‑95 (MQ=255)            caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCt                                                                                   >  2:432007/1‑100 (MQ=255)              ccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTc                                                                                  <  2:822047/99‑1 (MQ=255)               cTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCg                                                                                >  2:877143/1‑100 (MQ=255)               cTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGGTCGTTAAGTTTGGGTTTCATGTCAGCTCCg                                                                                >  2:736382/1‑100 (MQ=255)                tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCgg                                                                               >  2:447347/1‑100 (MQ=255)                tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCgg                                                                               >  2:810705/1‑100 (MQ=255)                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGc                                                                              >  1:1608149/1‑100 (MQ=255)                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGc                                                                              <  1:618420/100‑1 (MQ=255)                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGc                                                                              <  2:329656/100‑1 (MQ=255)                  aTCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcg                                                                             <  1:1004524/100‑1 (MQ=255)                    cAACCAGCTCAGGTATCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgta                                                                           <  2:1917672/100‑1 (MQ=255)                    cAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgta                                                                           <  2:1016499/100‑1 (MQ=255)                    cAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgta                                                                           >  1:1917659/1‑100 (MQ=255)                      aCCAGTTCAGGTAACGCATGCCCTATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaac                                                                         <  1:444811/100‑1 (MQ=255)                      aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaac                                                                         <  2:1420644/100‑1 (MQ=255)                       ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacg                                                                        <  2:1800043/100‑1 (MQ=255)                       ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacg                                                                        >  1:812357/1‑100 (MQ=255)                       ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacg                                                                        >  1:994247/1‑100 (MQ=255)                           cTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAAt                                                                    >  1:247084/1‑100 (MQ=255)                            tCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAAtt                                                                   >  1:810682/1‑100 (MQ=255)                             cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                  <  1:1219913/100‑1 (MQ=255)                             cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                  >  2:1860965/1‑100 (MQ=255)                             cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                  <  2:350483/100‑1 (MQ=255)                               ggTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTaaa                                                                >  2:1673122/1‑100 (MQ=255)                                 tctcGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                              <  1:1886022/97‑1 (MQ=255)                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                              <  2:838029/100‑1 (MQ=255)                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                              <  2:1610325/100‑1 (MQ=255)                                  aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAg                                                             <  2:1105631/100‑1 (MQ=255)                                  aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAg                                                             <  1:1170130/100‑1 (MQ=255)                                   aCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGc                                                            <  1:1162881/100‑1 (MQ=255)                                    cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATaaga                                                           >  2:1886035/1‑97 (MQ=255)                                     gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCaa                                                          >  1:940249/1‑100 (MQ=255)                                      ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                         >  1:823810/1‑100 (MQ=255)                                       cTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACt                                                        <  2:1103327/100‑1 (MQ=255)                                        tGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                       >  1:391305/1‑100 (MQ=255)                                         gCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                      <  1:1897549/100‑1 (MQ=255)                                           cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCg                                                    >  2:882164/1‑100 (MQ=255)                                            ccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGc                                                   >  1:1220991/1‑100 (MQ=255)                                                  cTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAg                                             <  2:1220986/100‑1 (MQ=255)                                                   tttGTAAGTAAAAGGAGTAGGGAATTCAAGTAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAgg                                            >  2:1115038/1‑100 (MQ=255)                                                     tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                          <  2:382167/100‑1 (MQ=255)                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGccttctggcgg                                         <  1:643608/100‑10 (MQ=255)                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                         >  2:899389/1‑100 (MQ=255)                                                       taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                        <  2:1694808/100‑1 (MQ=255)                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGa                                       <  2:421441/100‑1 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                      <  1:293010/100‑1 (MQ=255)                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                      >  2:1896843/1‑99 (MQ=255)                                                           taaAAGGAGTAGGGAATTCAAGCAGTTCGTAAAATTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAGGGCGGGAGtt                                    >  1:1279540/1‑100 (MQ=255)                                                                ggAGTTGTGAATTCAAGCAGTTCATTAAGTTTTGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                <  2:1624564/99‑1 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                               <  2:106292/100‑1 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                               <  1:56573/100‑1 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                               >  2:752464/1‑100 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                               >  1:1943027/1‑100 (MQ=255)                                                                 gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCTACTCCCGCCAGAAGGCGGGAGTTTTTTAc                              <  1:736381/100‑1 (MQ=255)                                                                 gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                              <  2:1653001/100‑1 (MQ=255)                                                                 gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                              <  2:1486244/100‑1 (MQ=255)                                                                   gTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTg                            >  2:1077175/1‑100 (MQ=255)                                                                    tAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACtaa                           >  1:1179842/1‑98 (MQ=255)                                                                    tAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                           >  2:505157/1‑100 (MQ=255)                                                                    tAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                           <  2:537633/100‑1 (MQ=255)                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGTAGTTTTTTACTGAt                          >  1:1321937/1‑100 (MQ=255)                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                          >  2:1164223/1‑100 (MQ=255)                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                          >  1:1587478/1‑100 (MQ=255)                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTg                            >  1:265322/1‑97 (MQ=255)                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                         >  2:1874235/1‑100 (MQ=255)                                                                        gAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                       >  1:713569/1‑100 (MQ=255)                                                                        gAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                       <  2:1230621/100‑1 (MQ=255)                                                                         aaTTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGATTTTTACTGATGCtt                      >  1:1936916/1‑100 (MQ=255)                                                                           ttCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                    <  2:1169571/100‑1 (MQ=255)                                                                           ttCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                    <  1:1591154/100‑1 (MQ=255)                                                                            tCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTATATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAgn                   <  1:1080979/100‑2 (MQ=255)                                                                            tCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                   >  1:315506/1‑100 (MQ=255)                                                                            tCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                   >  1:315426/1‑100 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                  <  2:1137489/100‑1 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                  <  1:746034/100‑1 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                  <  2:810683/100‑1 (MQ=255)                                                                               aGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtata                >  2:927057/1‑100 (MQ=255)                                                                                   gTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATggg            <  2:1085832/100‑1 (MQ=255)                                                                                    ttCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg           <  1:1048808/100‑1 (MQ=255)                                                                                    ttCGTTAAGTTTGGGTTGCATGTCAGGTGGGGTGTAATGTAATTAAATAGCATCTCCCGCCAGAAGGCGGGAGGTTTTTACTGATGCTTAGTATATgggg           <  1:1254344/100‑1 (MQ=255)                                                                                      cGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc         >  2:560414/1‑100 (MQ=255)                                                                                       gTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACg        <  2:498697/100‑1 (MQ=255)                                                                                           attttGGGTTTCAAGTCAGCTCCGGCGTAACATAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa    >  2:1313154/3‑100 (MQ=255)                                                                                             tttGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAtt  <  2:979515/100‑1 (MQ=255)                                                                                                  |                                                                                               ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATT  >  NC_000913/662539‑662730

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