Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,900,481	G→T	intergenic (‑208/‑111)	ygcW ← / → yqcE	putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,900,481	0	G	T	90.6%	258.8 / 16.6	86	intergenic (‑208/‑111)	ygcW/yqcE	putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (6/2);  new base T (33/44);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTCTATTAC  >  NC_000913/2900384‑2900576                                                                                                   |                                                                                                ggCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCtcttct                                                                                               <  1:793441/100‑1 (MQ=255) ggCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCtcttct                                                                                               <  1:1243243/100‑1 (MQ=255)    aTAATTGCTAATGTTATCGTAATCATAGTGTGATATTCGTATTGTCAACTAATTTTAAATAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTc                                                                                            <  1:297831/100‑1 (MQ=255)    aTAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTc                                                                                            <  1:845633/100‑1 (MQ=255)      ctcTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCa                                                                                          <  2:752028/98‑1 (MQ=255)      aaCTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCa                                                                                          <  1:1299988/100‑1 (MQ=255)       aCTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCAt                                                                                         >  2:108816/1‑100 (MQ=255)        cTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCaag                                                                                        >  1:752027/1‑98 (MQ=255)        cTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATg                                                                                        >  2:1201931/1‑100 (MQ=255)        cTGCTAATGTGATCGTAATAACAGTGTGATATTCGTATTGTCTAATAATTATAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCCTCATTCCATg                                                                                        <  1:1363172/100‑1 (MQ=255)         tGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGt                                                                                       >  1:149341/1‑100 (MQ=255)          gCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGtt                                                                                      >  2:270490/1‑100 (MQ=255)             aaTGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGttatt                                                                                   >  1:1050616/1‑100 (MQ=255)              aTGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGttattt                                                                                  <  2:14272/100‑1 (MQ=255)               tgtgATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGttattta                                                                                 >  1:1286431/1‑100 (MQ=255)               tgtgATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGttattta                                                                                 >  1:312779/1‑100 (MQ=255)                  gATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATttatta                                                                              >  2:1419015/1‑100 (MQ=255)                  gATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATttatta                                                                              >  2:452521/1‑100 (MQ=255)                   aTCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTtattat                                                                             <  2:201068/100‑1 (MQ=255)                   aTCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTtattat                                                                             >  2:855680/1‑100 (MQ=255)                    tCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTattatt                                                                            >  2:1740353/1‑100 (MQ=255)                    tCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTattatt                                                                            <  1:919243/100‑1 (MQ=255)                       tAATCACAGTGTGGTATTCGTATTGTCAAATAATTTTAGAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCa                                                                         >  1:1201879/1‑100 (MQ=255)                       tAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAACGCAGAAGCGTGTTTACTCTCACA‑CAACAACCCTTCTTTCCATGTTATTTATTATTTCa                                                                         <  1:89858/100‑1 (MQ=255)                       tAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCa                                                                         >  2:1801788/1‑100 (MQ=255)                          tCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAAtt                                                                      <  1:728295/100‑1 (MQ=255)                           cacaGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTa                                                                     >  1:57927/1‑100 (MQ=255)                            acaGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTAc                                                                    <  1:1659369/100‑1 (MQ=255)                             caGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACt                                                                   <  2:980516/100‑1 (MQ=255)                              aGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTa                                                                  >  1:500071/1‑100 (MQ=255)                              aGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTa                                                                  <  1:1589399/100‑1 (MQ=255)                                tgtgATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTaga                                                                >  1:254665/1‑100 (MQ=255)                                 gtgATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTagag                                                               <  1:1615665/100‑1 (MQ=255)                                    atatTCGTATTGTCACATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGAcc                                                            <  2:1916269/100‑1 (MQ=255)                                      atTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTg                                                          <  2:892713/100‑1 (MQ=255)                                       ttCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTga                                                         >  1:1178501/1‑100 (MQ=255)                                       ttCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTga                                                         <  1:108816/100‑1 (MQ=255)                                        tCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTgag                                                        <  2:1178503/100‑1 (MQ=255)                                         cGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTgaga                                                       <  1:1511638/100‑1 (MQ=255)                                             ttGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAg                                                   <  2:1422517/100‑1 (MQ=255)                                             ttGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAg                                                   <  2:320613/100‑1 (MQ=255)                                              tGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGt                                                  <  2:1070625/100‑1 (MQ=255)                                              tGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGt                                                  <  1:921297/100‑1 (MQ=255)                                               gTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTa                                                 >  1:1078751/1‑100 (MQ=255)                                               gTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTa                                                 >  2:468807/1‑100 (MQ=255)                                                  aaaTAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATc                                              >  2:1544051/1‑100 (MQ=255)                                                  aaaTAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATc                                              >  2:1543916/1‑100 (MQ=255)                                                    ataatTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATcc                                             >  1:1385057/1‑99 (MQ=255)                                                     taatTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCAc                                           >  2:682770/1‑100 (MQ=255)                                                       atTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCt                                         >  2:1892140/1‑100 (MQ=255)                                                       atTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCt                                         <  1:1768476/100‑1 (MQ=255)                                                        ttttAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTg                                        <  1:1265764/100‑1 (MQ=255)                                                        ttttAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTg                                        >  1:809718/1‑100 (MQ=255)                                                         tttAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGa                                       <  1:1391626/100‑1 (MQ=255)                                                          ttAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGaa                                      >  2:1110172/1‑100 (MQ=255)                                                           tAAAAAGCAGAAGCCGTATTACTTTCACACCACCACCCCT‑CCTTCCATGTTATTTTTTTTTTCCAATACTTGAGACCTTAGAAAAGTTATTCACCTTaaa                                     >  1:1711385/1‑100 (MQ=255)                                                            aaaaaGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAg                                    >  1:1940846/1‑100 (MQ=255)                                                            aaaaaGCAGAAGCGTGATTACTCTCACA‑CAACACCCCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGaaa                                     >  2:1469409/1‑99 (MQ=255)                                                             aaaaGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAgg                                   <  2:91377/100‑1 (MQ=255)                                                              aaaGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAggg                                  <  2:500073/100‑1 (MQ=255)                                                              aaaGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAggg                                  >  1:123445/1‑100 (MQ=255)                                                               aaGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGggc                                 <  2:57927/100‑1 (MQ=255)                                                                 gCAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGggcgg                               >  2:478296/1‑100 (MQ=255)                                                                  cAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGggcggc                              <  1:274829/100‑1 (MQ=255)                                                                  cAGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGggcggc                              <  2:1731925/100‑1 (MQ=255)                                                                   aGAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCt                             <  1:1740342/100‑1 (MQ=255)                                                                    gAAGCGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCtt                            <  2:532128/100‑1 (MQ=255)                                                                        cGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACtt                        <  1:73553/100‑1 (MQ=255)                                                                        cGTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACtt                        >  1:522468/1‑100 (MQ=255)                                                                         gTGATTACTCTCCCA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACttt                       <  1:625178/100‑1 (MQ=255)                                                                         gTGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACttt                       <  1:1003096/100‑1 (MQ=255)                                                                          tGATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACtttt                      <  2:1233091/100‑1 (MQ=255)                                                                           gATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTtc                     >  2:889854/1‑100 (MQ=255)                                                                           gATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTtc                     >  1:814737/1‑100 (MQ=255)                                                                           gATTACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTtc                     <  2:760486/100‑1 (MQ=255)                                                                             ttACTCTCACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTtct                    <  1:1595696/99‑1 (MQ=255)                                                                              tACTCTCACA‑CAACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTctct                  >  1:400390/1‑100 (MQ=255)                                                                                  ctcACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCttttt              >  2:1195068/1‑100 (MQ=255)                                                                                   tcACA‑CAACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCtttttt             <  2:1337159/100‑1 (MQ=255)                                                                                     aca‑caACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCtttttttn           <  2:278686/100‑2 (MQ=255)                                                                                     aca‑caACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCtttttttt           >  1:1164381/1‑100 (MQ=255)                                                                                      ca‑caACACCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTg          <  2:1098109/100‑1 (MQ=255)                                                                                         caacaCCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTc        <  1:1419007/100‑1 (MQ=255)                                                                                         caacaCCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTc        <  2:839926/100‑1 (MQ=255)                                                                                         caacaCCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTc        <  1:855679/100‑1 (MQ=255)                                                                                         caacaCCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTc        <  2:1229189/100‑1 (MQ=255)                                                                                          aacaCCTCTTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAACGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTCt       >  1:271134/1‑100 (MQ=255)                                                                                               ctctTCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTCTATTAc  >  1:177237/1‑100 (MQ=255)                                                                                                   |                                                                                                GGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACA‑CAACACCTCTGCTTTCCATGTTATTTATTATTTCAATTACTAGAGACCTGAGAAAAGTAATCCACCTGAAAGGGCGGCTTACTTTTCTCTTTTTTTTGTCTATTAC  >  NC_000913/2900384‑2900576

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