Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,820	G→A	56.9%	R269C (CGC→TGC)	purU ←	formyltetrahydrofolate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,820	0	G	A	56.9%	25.9 / 101.5	86	R269C (CGC→TGC)	purU	formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (26/23);  total (45/41)
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.46e-01

ACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACA  >  NC_000913/1287722‑1287917                                                                                                   |                                                                                                  aCTGGTTGCGTAATTTTCATCCATAACGGATAAAAGGAAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACgcgc                                                                                                  <  1:1749871/100‑1 (MQ=255)     gTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCtgtg                                                                                              <  2:1497278/100‑1 (MQ=255)       tGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGcc                                                                                            >  1:1007428/1‑100 (MQ=255)       tGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGcc                                                                                            <  1:1039740/100‑1 (MQ=255)        gCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCa                                                                                           >  2:588260/1‑100 (MQ=255)         cGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAg                                                                                          <  2:1675193/100‑1 (MQ=255)          gTAATTTTCATCCGTAACGGATTAAAGGTAAGCAGATATTTTTGCTGGCGATTGGATAATAATAGCTGGATTACCGGAACCAAAGAACCGGtgtgaaa                                                                                           >  1:2121554/1‑95 (MQ=255)          gTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAg                                                                                          <  1:1353963/99‑1 (MQ=255)           tAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTTTTTTTTTTTTTAAAATAAACGTTCGATTACCGTATACAAAGTCGCACtgtgtccgt                                                                                         >  1:185148/1‑94 (MQ=255)            aaTTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTAc                                                                                       >  2:1500191/1‑100 (MQ=255)              ttttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACtt                                                                                     >  2:344943/1‑100 (MQ=255)                ttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTg                                                                                   <  2:868840/100‑1 (MQ=255)                ttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCACTACTTTg                                                                                   <  2:1357587/100‑1 (MQ=255)                ttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCCATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTg                                                                                   <  1:1383006/100‑1 (MQ=255)                 tCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGt                                                                                  <  2:2270542/100‑1 (MQ=255)                  cATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGta                                                                                 <  1:14279/100‑1 (MQ=255)                  cATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGta                                                                                 >  2:1379143/1‑100 (MQ=255)                    tCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGtata                                                                               >  1:338525/1‑100 (MQ=255)                     ccGTAACGGATTAAAGGTAACCAGTTATTTTTTTTTTTTTTTTTAAAATAAACGTGCGATTACCGTATACCACGACGCACTGTGCCAGTACTTTGta                                                                                 >  1:523832/1‑97 (MQ=255)                      cGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGt                                                                             >  1:2250708/1‑100 (MQ=255)                       gTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTg                                                                            >  1:1484102/1‑100 (MQ=255)                       gTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGAGTAAAGAATAATCATTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTTTATAGt                                                                             <  1:844814/99‑1 (MQ=255)                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGc                                                                           >  1:1358445/1‑100 (MQ=255)                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGc                                                                           >  2:1225452/1‑100 (MQ=255)                         aaCGGATTAGAGGTAACCAGGTATTTTGGCTGGCGAGTAAAGAATAATCGTATGATTACGGTGAAGAAAGATGCAATGTGCCAGTACTTTGTATAGTGc                                                                           >  1:1782360/1‑99 (MQ=255)                          aCGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCAc                                                                         <  2:1358453/100‑1 (MQ=255)                           cGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACg                                                                        <  2:2179599/100‑1 (MQ=255)                              aTTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACt                                                                     >  2:1315460/1‑100 (MQ=255)                              aTTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACt                                                                     >  2:1497109/1‑100 (MQ=255)                               ttAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACtt                                                                    <  2:1078037/100‑1 (MQ=255)                                 aaaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGGCGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTaa                                                                  >  2:1786874/1‑100 (MQ=255)                                  aaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAg                                                                 >  2:2019626/1‑100 (MQ=255)                                  aaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAg                                                                 >  2:1549600/1‑100 (MQ=255)                                     gTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACg                                                              <  2:2035965/100‑1 (MQ=255)                                        aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGttt                                                           <  1:60593/100‑1 (MQ=255)                                           aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                        >  2:1714272/1‑100 (MQ=255)                                            gTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtct                                                       <  1:2585899/100‑1 (MQ=255)                                             ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGTTTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtctc                                                      <  1:1361837/100‑1 (MQ=255)                                             ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtctc                                                      <  1:1304211/100‑1 (MQ=255)                                              tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCg                                                     <  1:347795/100‑1 (MQ=255)                                              tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCg                                                     <  2:1718259/100‑1 (MQ=255)                                              tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGGACGACTTAAGACGTTTTTCTCg                                                     <  1:1500182/100‑1 (MQ=255)                                              tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCg                                                     >  1:1733991/1‑100 (MQ=255)                                              tATTTTGGCTGGCGATTAAAGGATTATCGTTCGATTGCCGTAAACAAAGACGCACTGTGCCAGGACTTTGTATAGTGCACGACTTGAGACGTTTTTCTCg                                                     <  1:415322/100‑1 (MQ=255)                                               aTTTTTGGTGGCGATTAAAGAATAATCGTTCGATTATCGTAAACAAAGACGAGCTGTGCGTGTACTTTGTATAGTGCACGACTTAAGACGttttt                                                         >  1:2127673/1‑95 (MQ=255)                                                tttttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGAc                                                   <  1:1949960/100‑1 (MQ=255)                                                 ttttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACg                                                  >  2:1675537/1‑100 (MQ=255)                                                   ttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTc                                                <  2:976185/100‑1 (MQ=255)                                                   ttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTc                                                >  1:2012282/1‑100 (MQ=255)                                                     gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCAc                                              >  2:2531599/1‑100 (MQ=255)                                                     gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCAc                                              >  2:2346269/1‑100 (MQ=255)                                                      cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACg                                             <  1:1241454/100‑1 (MQ=255)                                                      cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACg                                             <  1:1496428/100‑1 (MQ=255)                                                        ggCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATTGTGCACGACTTAAGACGTCTTTCTCGACGTCACGAc                                           >  2:731250/1‑100 (MQ=255)                                                          agATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCt                                         >  2:876490/2‑100 (MQ=255)                                                               aaaGAATAATCGTTCGATTACCGTAAACAAAGGCGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                                    <  1:1786865/100‑1 (MQ=255)                                                               aaaGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                                    <  1:1840654/100‑1 (MQ=255)                                                                 acggtaatCGAACGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGca                                  >  2:183991/5‑100 (MQ=255)                                                                 aGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGca                                  >  1:736609/1‑100 (MQ=255)                                                                  gAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcat                                 >  1:493742/1‑100 (MQ=255)                                                                   aataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                                <  1:2317206/100‑1 (MQ=255)                                                                    ataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatca                               >  1:1268923/1‑100 (MQ=255)                                                                    ataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatca                               <  2:2563561/100‑1 (MQ=255)                                                                      aatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTAGACGTCACGACCTGCGCGCATCata                             >  2:544275/1‑100 (MQ=255)                                                                        tCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATc                           <  2:1438023/100‑1 (MQ=255)                                                                         cGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCt                          >  2:1429284/1‑100 (MQ=255)                                                                           ttCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTc                        <  2:1359365/100‑1 (MQ=255)                                                                            tCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTAGAGTGCACGACTTAAGACGTTTTTCTCAACGTCACGACCTGCGCGCATCATATCTTc                        <  1:1844205/99‑1 (MQ=255)                                                                               aTTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                    >  2:1882617/1‑100 (MQ=255)                                                                               aTTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                    >  1:1841325/1‑100 (MQ=255)                                                                               aTTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTGTCTCGACGGCACGACCTGCGAGCATCACATc                           <  1:1396448/93‑1 (MQ=255)                                                                               aTTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                    <  2:2250723/100‑1 (MQ=255)                                                                                ttACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtg                   <  2:2399907/100‑1 (MQ=255)                                                                                   ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt                >  2:96549/1‑100 (MQ=255)                                                                                   ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt                >  2:849801/1‑100 (MQ=255)                                                                                     gTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg              <  1:2408468/100‑1 (MQ=255)                                                                                     gTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg              >  1:2522120/1‑100 (MQ=255)                                                                                      tAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAgg             >  2:2148277/1‑100 (MQ=255)                                                                                       aaaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGt            >  1:2436924/1‑100 (MQ=255)                                                                                        aaCAATGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACTTCACGACCTGCGCGCATCATATCTTGAGCGCTGTAGGTa           >  1:1287292/1‑100 (MQ=255)                                                                                        aaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTa           >  1:795209/1‑100 (MQ=255)                                                                                         aCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTAt          <  2:1914984/100‑1 (MQ=255)                                                                                          cAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATg         >  1:1947706/1‑100 (MQ=255)                                                                                          cAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATg         <  2:1316334/100‑1 (MQ=255)                                                                                          cAAAGACGCACTGAGCCAGTACTTTGTATAGTGCACGAATTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATg         >  1:1600426/1‑100 (MQ=255)                                                                                           aaaGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa        <  2:609849/100‑1 (MQ=255)                                                                                              gACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCg     <  1:2132268/100‑1 (MQ=255)                                                                                                 gcgcTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACa  <  2:1530303/100‑1 (MQ=255)                                                                                                   |                                                                                                  ACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACA  >  NC_000913/1287722‑1287917

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: ‑