Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	662,636	T→G	100%	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.7%	296.4 / ‑2.9	87	T3P (ACC→CCC)	ybeD	UPF0250 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (43/42);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTAC  >  NC_000913/662539‑662732                                                                                                  |                                                                                                 aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgtt                                                                                                <  1:1849863/100‑1 (MQ=255)    aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                             >  2:2362548/1‑100 (MQ=255)     cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCAt                                                                                            >  2:1916287/1‑100 (MQ=255)      ttCAACCACCTGATCAACCAGCTCAGGTAACTCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATg                                                                                           <  2:1725086/100‑1 (MQ=255)      ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATg                                                                                           <  2:1857174/100‑1 (MQ=255)       tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGt                                                                                          <  1:1519340/100‑1 (MQ=255)        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTc                                                                                         <  2:1859978/100‑1 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTctg                                                                                       <  1:1724434/100‑3 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAg                                                                                       >  1:2302443/1‑100 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAg                                                                                       <  2:2306338/100‑1 (MQ=255)            caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCt                                                                                     >  2:712910/1‑100 (MQ=255)             accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTCTCATGTCAGCTc                                                                                    <  1:108194/100‑1 (MQ=255)                tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCgg                                                                                 >  2:106336/1‑100 (MQ=255)                tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCgg                                                                                 >  1:1242217/1‑100 (MQ=255)                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGGTAAGTATGGGTTTCATGTCAGCTCCg                                                                                  >  2:2286248/1‑98 (MQ=255)                   tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcg                                                                               <  2:1921618/99‑1 (MQ=255)                   tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcg                                                                               >  1:458391/1‑99 (MQ=255)                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaa                                                                            >  1:775787/1‑100 (MQ=255)                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaa                                                                            >  2:688002/1‑100 (MQ=255)                      ggttgatCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaac                                                                           >  1:1429323/7‑100 (MQ=255)                      aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCACGCAGTTCGTTAAGTTTGGGGTTCATGTCAGCTCCGGcgtaac                                                                           >  2:1658412/1‑100 (MQ=255)                      aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCTGTTCGTCAAGTTTGGGTTTCATGTCAGCTCCGGcgtaac                                                                           >  2:1801199/1‑100 (MQ=255)                      aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaac                                                                           >  1:148739/1‑100 (MQ=255)                       ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacg                                                                          <  1:852382/100‑1 (MQ=255)                        cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacgt                                                                         <  2:1644692/100‑1 (MQ=255)                           cTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAAt                                                                      >  2:469275/1‑100 (MQ=255)                            tCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAAtt                                                                     <  1:2442174/100‑1 (MQ=255)                             cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                    <  1:758180/100‑1 (MQ=255)                             cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                    >  1:1497381/1‑100 (MQ=255)                              aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTaa                                                                   >  1:668304/1‑100 (MQ=255)                                gTAACGCCTGCCCCGTAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAAt                                                                 <  2:967945/100‑1 (MQ=255)                                gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAAt                                                                 <  2:198915/100‑1 (MQ=255)                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                                <  1:1180074/100‑1 (MQ=255)                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                                <  2:1228926/100‑1 (MQ=255)                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCAGGCGTAACGTAATTAAATa                                                                >  2:1525752/1‑100 (MQ=255)                                  aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAGATAg                                                               >  1:1583475/1‑100 (MQ=255)                                   aCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGc                                                              <  1:342561/100‑1 (MQ=255)                                     gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCaa                                                            <  1:1916301/100‑1 (MQ=255)                                      ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                           >  2:133527/1‑100 (MQ=255)                                        tGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                         >  2:1286383/1‑100 (MQ=255)                                         gCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                        <  1:243253/100‑1 (MQ=255)                                           cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCg                                                      >  2:1540014/1‑100 (MQ=255)                                           cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCg                                                      >  2:1557871/1‑100 (MQ=255)                                                     tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                            <  1:2081937/100‑1 (MQ=255)                                                     tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                            <  2:1317566/100‑1 (MQ=255)                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTACCGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                           >  2:282265/1‑100 (MQ=255)                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                           <  1:694454/100‑1 (MQ=255)                                                       taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGccttctggcggg                                          <  1:200923/100‑11 (MQ=255)                                                       taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                          >  1:1524419/1‑100 (MQ=255)                                                       taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                          <  2:1353643/100‑1 (MQ=255)                                                       taagtaaAACGAGTAGGGAATTCAAGCAGTCCATTCAGTTTGGGCTTCATGTCAGCACCGGCGTACCACAATCACATAGCCActccctcccgcagacg                                            >  2:1416366/1‑90 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                        <  1:538417/100‑1 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                        >  1:1310003/1‑100 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                        <  2:2044147/100‑1 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTAGGTTTCATGTCAGCTCAGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                        <  2:1394733/100‑1 (MQ=255)                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGGTCCGGCGTACCGTAAGTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                       >  2:1730740/1‑100 (MQ=255)                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                       >  2:1069474/1‑100 (MQ=255)                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                       <  2:39857/100‑1 (MQ=255)                                                               cgatctAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                  <  1:1208081/95‑1 (MQ=255)                                                                ggAGTGGGGAATCCAATCAGTTCGTTAAGTTTGGGTTTCATGGCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                  <  2:445413/99‑1 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                  <  2:2457889/99‑1 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                  >  1:1260887/1‑99 (MQ=255)                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                 >  2:1808768/1‑100 (MQ=255)                                                                 gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACGCCCGCCAGAAGGCGGGAGTTTTTTAc                                <  2:60302/100‑1 (MQ=255)                                                                  aGTAGGGAATTCAAGCAGTTCGTTAAATTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACt                               <  2:141264/100‑1 (MQ=255)                                                                   gTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAACTTCGACGTAACATAATTAAATAGCAACTCCCGACAGAAAGCAGGAATTTTTTACt                               >  2:921098/1‑99 (MQ=255)                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTagatcg                            >  2:1208073/1‑95 (MQ=255)                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGGGTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                             >  2:2413791/1‑99 (MQ=255)                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                           >  2:1922355/1‑100 (MQ=255)                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                           <  2:192863/100‑1 (MQ=255)                                                                        gAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                         >  1:542681/1‑100 (MQ=255)                                                                          aTTCAAGCAGTTCGTTAAGTTTGGGTTTTATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                       <  1:269071/100‑1 (MQ=255)                                                                           ttCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                      >  2:1456176/1‑100 (MQ=255)                                                                           ttCAAGCAGTTCATTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGAAATTAAATGGCAACTGCCGCCAGCAGGCGGGTGTTTTTTACTGCTGCTTa                       >  2:690318/1‑99 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                    <  2:1062810/100‑1 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                    >  2:1886047/1‑100 (MQ=255)                                                                              aaGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtat                   <  1:165469/100‑1 (MQ=255)                                                                              aaGCAGTTCGTTAAGTTTGAGTTTCATGTCAGCTCCGGCGTAACGTAATTAACTAGCAACTCCCGCCACAAGGCGGGAGTTTTTTACCGATGCTCAGtat                   <  2:115416/100‑1 (MQ=255)                                                                                   gTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATggg              <  1:935647/100‑1 (MQ=255)                                                                                     tCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGa            >  2:2381101/1‑100 (MQ=255)                                                                                      cGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc           >  1:1608866/1‑100 (MQ=255)                                                                                      cGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc           >  1:1008077/1‑100 (MQ=255)                                                                                      cGTTAAGTTTGGGTTTCATGTCAGCAGCGGCGTAGCGTAATTAAATGGCATCTGCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATg                <  2:1974929/95‑1 (MQ=255)                                                                                       gTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACg          <  2:2233260/100‑1 (MQ=255)                                                                                        ttAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACgg         <  1:2266867/100‑1 (MQ=255)                                                                                          aaGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaa       <  1:43012/100‑1 (MQ=255)                                                                                           aGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa      <  1:1564036/100‑1 (MQ=255)                                                                                               tGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTac  >  1:1697000/1‑100 (MQ=255)                                                                                                  |                                                                                                 ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTAC  >  NC_000913/662539‑662732

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