Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,758	T→C	intergenic (+129/‑7)	lysT → / → valT	tRNA‑Lys/tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,758	0	T	C	97.9%	312.2 / ‑2.9	93	intergenic (+129/‑7)	lysT/valT	tRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (43/48);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCA  >  NC_000913/780663‑780837                                                                                                |                                                                                ggATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgcagta                                                                             >  2:797022/1‑100 (MQ=255)   aTAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATc                                                                           <  1:1069208/100‑1 (MQ=255)   aTAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATc                                                                           <  1:1067584/100‑1 (MQ=255)    tAACGTTGCATGACCAACGGCCCGAAGGGCGAGACTAAGTCGAGTCATCCTGCACCACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCg                                                                          <  1:1166865/100‑1 (MQ=25)     aaCGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCgg                                                                         <  1:379813/100‑1 (MQ=255)      aCGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCggg                                                                        >  1:2201102/1‑100 (MQ=255)       cGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGt                                                                       <  1:2585580/100‑1 (MQ=255)       cGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGt                                                                       >  1:736451/1‑100 (MQ=255)       cGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGt                                                                       <  2:1805092/100‑1 (MQ=255)       cGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGt                                                                       <  2:2401601/100‑1 (MQ=255)       cGTTGCATGAGCAACGGCCCGAACGGCGAGATGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCCGTATCGGGt                                                                       <  1:2761583/100‑1 (MQ=25)         ttGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGa                                                                     <  2:2386664/100‑1 (MQ=255)         ttGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGa                                                                     <  1:452987/100‑1 (MQ=255)          tGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGAt                                                                    <  1:1714152/100‑1 (MQ=255)          tGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGAt                                                                    >  1:2050020/1‑100 (MQ=255)           gCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGAtt                                                                   <  1:1767239/100‑1 (MQ=255)           gCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGAtt                                                                   >  2:752717/1‑100 (MQ=255)            ctcatgCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTa                                                                  >  1:461013/6‑100 (MQ=255)             aTGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAg                                                                 >  2:1996014/1‑100 (MQ=255)              tGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGc                                                                >  1:1540523/1‑100 (MQ=255)              catgCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGc                                                                >  2:2135214/4‑100 (MQ=255)               gAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCt                                                               <  1:2336385/100‑1 (MQ=255)                aGCAACGGCCGAAAGCGCGAAACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGc                                                                <  1:1144409/98‑1 (MQ=25)                 gCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCa                                                             >  2:2484508/1‑100 (MQ=255)                 gCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCa                                                             >  2:774027/1‑100 (MQ=255)                      ggCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTgg                                                        <  2:474286/100‑1 (MQ=255)                      ggCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTgg                                                        <  1:1090223/100‑1 (MQ=255)                      ggCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTgg                                                        >  2:2202584/1‑100 (MQ=255)                       gCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTggg                                                       <  2:703828/100‑1 (MQ=255)                       gCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTggg                                                       <  1:362046/100‑1 (MQ=255)                        cccGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGga                                                      >  2:758580/1‑100 (MQ=255)                         cttccgatctagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGgag                                                     <  1:1930875/90‑1 (MQ=25)                          cGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGgaga                                                    <  2:1302460/100‑1 (MQ=255)                           gAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGagag                                                   <  1:453924/100‑1 (MQ=255)                            aaGGGCGAGACGAAGTCGAGTCATCCTGCAAGACACACCAATAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGc                                                  <  1:964426/100‑1 (MQ=255)                             aGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCa                                                 >  2:1398133/1‑100 (MQ=255)                             aGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCa                                                 <  2:857536/100‑1 (MQ=255)                              gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                <  1:1513856/100‑1 (MQ=255)                              gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                >  2:2331988/1‑100 (MQ=255)                              gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                <  2:736456/100‑1 (MQ=255)                              gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                <  2:686514/100‑1 (MQ=255)                              gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                >  2:2827287/1‑100 (MQ=255)                               ggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAcc                                               >  2:1739554/1‑100 (MQ=255)                                gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCt                                              >  2:870568/1‑100 (MQ=255)                                 cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCt                                              >  1:1663058/1‑99 (MQ=255)                                  gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                                            >  2:2339146/1‑100 (MQ=255)                                  gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                                            >  2:858102/1‑100 (MQ=255)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGagagatcggaag                                           >  2:1930903/1‑92 (MQ=25)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                           <  1:875737/100‑1 (MQ=255)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                           <  2:2574977/100‑1 (MQ=255)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                           <  1:1344314/100‑1 (MQ=255)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                           <  2:2356582/100‑1 (MQ=255)                                   agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                           <  2:99136/100‑1 (MQ=255)                                    gaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt                                          >  2:753834/1‑100 (MQ=255)                                      cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa                                        <  2:2333131/100‑1 (MQ=255)                                      cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa                                        <  2:1626339/100‑1 (MQ=255)                                      cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa                                        <  1:884334/100‑1 (MQ=255)                                         aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACaa                                     <  1:656428/100‑1 (MQ=255)                                            cGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                  >  2:111585/1‑100 (MQ=255)                                            cGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATCAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                  >  1:1222577/1‑100 (MQ=255)                                             gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggag                                 >  1:455335/1‑100 (MQ=255)                                             gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga                                  >  2:2661484/1‑99 (MQ=255)                                              aGTCATCCTGCACGACCCACCACTCACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggagg                                <  1:1162441/100‑1 (MQ=255)                                              aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggagg                                >  2:1331972/1‑100 (MQ=255)                                              aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggagg                                >  2:1331959/1‑100 (MQ=255)                                                  aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggggg                             <  1:1344978/99‑1 (MQ=255)                                                  aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGt                            >  2:2814512/1‑100 (MQ=255)                                                   tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTc                           <  1:2178695/100‑1 (MQ=255)                                                     cTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcgg                         <  2:324842/100‑1 (MQ=255)                                                     cTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcgg                         >  1:2293955/1‑100 (MQ=255)                                                      tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTACCTTACAAGGAGGGGGTcggc                        <  1:1279521/100‑1 (MQ=255)                                                         aCGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGt                     <  2:847667/100‑1 (MQ=255)                                                           gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTc                   <  2:902633/100‑1 (MQ=255)                                                             cccACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGa                 >  2:2525007/1‑100 (MQ=255)                                                             cccACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGa                 >  1:179942/1‑100 (MQ=255)                                                               caccacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATc               >  2:259489/1‑100 (MQ=255)                                                                accacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATcc              <  1:1165373/100‑1 (MQ=255)                                                                 ccacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATccc             <  2:1256341/100‑1 (MQ=255)                                                                 ccacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATccc             <  1:1250640/100‑1 (MQ=255)                                                                  cacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCg            <  2:707639/100‑1 (MQ=255)                                                                  cacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCg            >  2:2464878/1‑100 (MQ=255)                                                                    cTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtc          >  1:2016364/1‑100 (MQ=255)                                                                     tAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtca         <  2:378383/100‑1 (MQ=255)                                                                      aaCATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat        <  1:1083943/100‑1 (MQ=255)                                                                      aaCATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat        >  2:1389296/1‑100 (MQ=255)                                                                      aaCATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat        >  1:1153719/1‑100 (MQ=255)                                                                      aaCATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTAATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcat        <  1:2578331/100‑1 (MQ=255)                                                                       aCATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcatc       >  1:432349/1‑100 (MQ=25)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca  <  1:2638904/100‑1 (MQ=11)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca  >  1:1633234/1‑100 (MQ=11)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca  >  2:1529156/1‑100 (MQ=11)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca  >  2:1842876/1‑100 (MQ=11)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca  >  2:1842822/1‑100 (MQ=11)                                                                            gttagttGTAGTATCCAGCGCAGTATCGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGTGGTCGGCGGTTCGATCCCGTCATcaccca  >  1:1743397/1‑100 (MQ=255)                                                                                                |                                                                                GGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCA  >  NC_000913/780663‑780837

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