Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,386,365	T→A	20.5%	F169Y (TTC→TAC)	yfbL →	putative M28A family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,365	0	T	A	20.5%	181.5 / 43.5	84	F169Y (TTC→TAC)	yfbL	putative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base T (39/27);  new base A (10/7);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACG  >  NC_000913/2386266‑2386460                                                                                                    |                                                                                                gCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCctt                                                                                                 <  1:1419361/100‑1 (MQ=255)  cGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCcttc                                                                                                <  1:1066708/100‑1 (MQ=255)     tgtgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCttt                                                                                             <  1:720115/100‑1 (MQ=255)       tgacAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAATTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCg                                                                                           <  2:71492/97‑1 (MQ=255)       tgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCg                                                                                           >  1:413375/1‑100 (MQ=255)        ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGt                                                                                          >  1:1389923/1‑100 (MQ=255)        ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                          >  1:1231198/1‑100 (MQ=255)        ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                          <  2:253765/100‑1 (MQ=255)        ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                          >  2:10770/1‑100 (MQ=255)        ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                          >  1:2154338/1‑100 (MQ=255)         gCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTa                                                                                         >  2:1222176/1‑100 (MQ=255)            ggATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCg                                                                                      >  2:676937/1‑100 (MQ=255)            ggATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCg                                                                                      >  1:1319501/1‑100 (MQ=255)                tACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGa                                                                                  <  1:2212706/100‑1 (MQ=255)                 aCTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGaa                                                                                 >  1:2397833/1‑100 (MQ=255)                 aCTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGa                                                                                  >  2:892982/1‑99 (MQ=255)                   tCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAAt                                                                               <  1:2313551/100‑1 (MQ=255)                    cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATg                                                                              >  2:1679559/1‑100 (MQ=255)                     gAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATgg                                                                             <  2:1366464/100‑1 (MQ=255)                     gAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATgg                                                                             <  1:768623/100‑1 (MQ=255)                      aaCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATggg                                                                            <  1:1761136/100‑1 (MQ=255)                      aaCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATggg                                                                            <  2:1749688/100‑1 (MQ=255)                         tGGCATGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAg                                                                         <  1:1973612/100‑1 (MQ=255)                          ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAgc                                                                        >  2:1352290/1‑100 (MQ=255)                          ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAgc                                                                        >  2:627176/1‑100 (MQ=255)                          ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAg                                                                         >  1:1921398/1‑99 (MQ=255)                               gTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGt                                                                   >  2:984514/1‑100 (MQ=255)                                tttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtg                                                                  <  1:945623/100‑1 (MQ=255)                                tttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGt                                                                   >  1:578203/1‑99 (MQ=255)                                   gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                               >  2:401365/1‑100 (MQ=255)                                   gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                               >  2:2019642/1‑100 (MQ=255)                                   gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                               >  1:1908334/1‑100 (MQ=255)                                   gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                               >  2:1229662/1‑100 (MQ=255)                                    ttACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGGGCGCGGtgcatg                                                              <  1:2283174/100‑1 (MQ=255)                                     tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                             >  2:1625781/1‑100 (MQ=255)                                     tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                             >  2:1193710/1‑100 (MQ=255)                                       cATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAg                                                           >  1:2375424/1‑100 (MQ=255)                                        aTCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGc                                                          <  1:2365967/100‑1 (MQ=255)                                          cAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCt                                                         >  2:2385632/1‑99 (MQ=255)                                           aGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTc                                                       >  1:998180/1‑100 (MQ=255)                                             cAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGc                                                     >  1:1021116/1‑100 (MQ=255)                                              aaGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCt                                                    <  2:2311306/100‑1 (MQ=255)                                                gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                  <  2:1098710/100‑1 (MQ=255)                                                gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                  <  2:420392/100‑1 (MQ=255)                                                gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                  <  1:1658316/100‑1 (MQ=255)                                                 tACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGa                                                 <  2:1312457/100‑1 (MQ=255)                                                 naCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCATTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCt                                                    >  2:1363570/2‑97 (MQ=255)                                                  nccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                >  2:2280519/2‑100 (MQ=255)                                                  aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                <  1:1208217/100‑1 (MQ=255)                                                  aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                <  2:276478/100‑1 (MQ=255)                                                  aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                >  1:226594/1‑100 (MQ=255)                                                   ccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGc                                               >  2:1205226/1‑100 (MQ=255)                                                   ccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGc                                               >  1:538200/1‑100 (MQ=255)                                                    cGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTCAGCg                                              >  2:526440/1‑100 (MQ=255)                                                      aaaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGACGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTc                                            <  2:228966/100‑1 (MQ=255)                                                       aaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTcc                                           <  2:457095/100‑1 (MQ=255)                                                        cataGGCGACCAGCTGCACGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCa                                          >  2:874069/4‑100 (MQ=255)                                                        aaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCa                                          <  2:2375396/100‑1 (MQ=255)                                                          caagcattccgCTGGTCCGCTATGCGTCGGAAGAACCGACATTGTTTAGTAGCGAGGAACAGAGCAGCGCGGTGAATGGAGCTTAGCTTGAGCGAACAGt                                        >  2:2281641/11‑100 (MQ=255)                                                            ggCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGa                                      >  1:1589069/1‑100 (MQ=255)                                                            ggCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGa                                      >  2:1074107/1‑100 (MQ=255)                                                              cGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaaa                                    >  2:1200755/1‑100 (MQ=255)                                                               gTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAt                                   >  2:2336589/1‑100 (MQ=255)                                                               gTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAt                                   >  1:145020/1‑100 (MQ=255)                                                                tGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAtt                                  >  1:1265731/1‑100 (MQ=255)                                                                   aGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAAt                               <  2:112561/100‑1 (MQ=255)                                                                    gCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                              >  2:923913/1‑100 (MQ=255)                                                                     cTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGa                             >  1:1264986/1‑100 (MQ=255)                                                                      tGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGAt                            <  1:2300508/100‑1 (MQ=255)                                                                      tGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGAt                            <  1:2396610/100‑1 (MQ=255)                                                                        gTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAg                          <  2:413372/100‑1 (MQ=255)                                                                          cGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                        <  1:338050/100‑1 (MQ=255)                                                                          cGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                        >  1:1677357/1‑100 (MQ=255)                                                                             cTATGCGTCGGAAGAACCGCCCTCCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATACCAccg                     <  2:1304114/100‑3 (MQ=255)                                                                                tGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGga                   >  2:1860011/1‑99 (MQ=255)                                                                                  cGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgat                >  2:1446604/1‑100 (MQ=255)                                                                                      ggAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGAtt            <  2:617329/100‑1 (MQ=255)                                                                                       gaagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTg           <  2:2199451/100‑1 (MQ=255)                                                                                        aagaaCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGAtt            >  2:1427275/1‑98 (MQ=255)                                                                                          gaaCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCt        <  1:384420/100‑1 (MQ=255)                                                                                           aaCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTa       <  2:1346775/100‑1 (MQ=255)                                                                                            aCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAt      <  1:417729/100‑1 (MQ=255)                                                                                             ccGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAtt     >  1:2096397/1‑100 (MQ=255)                                                                                               gaccTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGACTATTAc   >  1:1619306/3‑100 (MQ=255)                                                                                               gCCCTTCTTTCGTAGCAATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTa    >  2:520431/1‑99 (MQ=255)                                                                                                cccTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACg  <  2:2397805/100‑1 (MQ=255)                                                                                                    |                                                                                                GCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACG  >  NC_000913/2386266‑2386460

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