Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,627,837	A→T	100%	I175F (ATC→TTC)	nadR →	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,837	0	A	T	96.7%	300.1 / 0.5	91	I175F (ATC→TTC)	nadR	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (40/48);  total (41/50)
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.63e-01

CAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATC  >  NC_000913/4627740‑4627932                                                                                                  |                                                                                                cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgttc                                                                                               <  2:2412049/100‑1 (MQ=255)   aCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTa                                                                                             <  2:1793008/100‑1 (MQ=255)     gAAGAGGGGATGGAGCCGTATCCTCACGGCTGGGATGTGGTGAGTAACGTCATCAAAAATTTTATGGCTGAAAAAGGGATCCAGCCGgatctgttct                                                                                              <  2:1650020/97‑1 (MQ=255)         agGGCATGGAGCCGTATCCGCACGGCTGGGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTc                                                                                       <  2:2101474/100‑1 (MQ=255)            atctGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGga                                                                                    <  1:1107659/97‑1 (MQ=255)              aTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaag                                                                                  <  1:1182177/100‑1 (MQ=255)              aTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaag                                                                                  <  1:2285870/100‑1 (MQ=255)               tGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaan                                                                                  <  2:1514401/99‑2 (MQ=255)               tGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaaga                                                                                 <  1:1961264/100‑1 (MQ=255)                ggAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaagat                                                                                >  2:1107651/1‑99 (MQ=255)                ggAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaagaa                                                                                >  1:2017090/1‑100 (MQ=255)                  aGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGc                                                                              <  1:2266302/100‑1 (MQ=255)                   gCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGcc                                                                             >  2:57540/1‑100 (MQ=255)                    ccgtttcttcACGGCTGGGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATGCAGCCGGATCTGTTCTACACCTCGGAAGAAGcc                                                                             <  1:370929/90‑1 (MQ=255)                    ccGTATCCGTACGACTGGGTTGTGTGGAGCATCGGCATCAAAAAGTTTATGGCTGAAAATGGGAGCCAGCCGGATCTGTTCTACACCTCGTAAGAAGcc                                                                             <  2:998339/99‑1 (MQ=255)                     cGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGa                                                                           >  1:1934374/1‑100 (MQ=255)                      gTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGAt                                                                          <  1:501506/100‑1 (MQ=255)                           cGCACGGCTGGGATGTGTGGAGCAACTGCATCAAAAAGTTTATGTCTGAAAAAGGGATCCAGCCCGATCTGTTCTACACCTCGGAAGAAGCCGATGCGcc                                                                     <  2:1262290/100‑1 (MQ=255)                             cACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCcac                                                                   <  1:938844/100‑1 (MQ=255)                               cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAg                                                                 <  1:2048550/100‑1 (MQ=255)                                 tcTGGGATGTGTGAAACACCGATATCAAAAAGTTTATGGCCGAAAAAGGCATCCAGCCGGATCTGTTCTACACCCCTGAAGAAGCCGATGCGACACAGta                                                               <  2:1395954/99‑1 (MQ=255)                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGta                                                               <  2:782449/100‑1 (MQ=255)                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGta                                                               <  1:676424/100‑1 (MQ=255)                                  cTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGtat                                                              <  1:692744/100‑1 (MQ=255)                                  cTGGGATGTGTGGAGCAACGGCATCAAAAAGAAAAAAAAAAAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGta                                                               <  2:2019092/99‑1 (MQ=255)                                   tGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGtata                                                             <  2:1061334/100‑1 (MQ=255)                                     ggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCCCAGTATATg                                                           >  2:1453760/1‑100 (MQ=255)                                     ggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                           <  2:2159038/100‑1 (MQ=255)                                     ggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                           >  2:1453680/1‑100 (MQ=255)                                      gATGTGTGGAGCAACGGCATCAAAAAGTTTATTGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                           <  1:2305321/99‑1 (MQ=255)                                      gATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATgg                                                          >  2:2477605/1‑100 (MQ=255)                                       attcgtcgAGCATAGGCACCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                           <  1:2305278/91‑1 (MQ=255)                                       aTGTGTGGAGCAACGGCATCAAAAGCTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCGCGGAAGAAGCCGATGCGCCACAGTATATGGa                                                         <  2:1067019/100‑1 (MQ=255)                                       aTGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGa                                                         <  2:1191110/100‑1 (MQ=255)                                         gtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAc                                                       >  1:317985/1‑100 (MQ=255)                                          tgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACa                                                      >  1:636415/1‑100 (MQ=255)                                          tgtgGAGCAACGGCATCAAAAAGTTTAGGGCGGAAAAAGGGATCCAGCCGGAGCGGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAc                                                       <  2:1482000/99‑1 (MQ=255)                                           gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACAt                                                     >  1:2228468/1‑100 (MQ=255)                                             ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                   >  2:151398/1‑100 (MQ=255)                                             ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                   >  1:70633/1‑100 (MQ=255)                                               aGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgg                                                 <  2:1441782/100‑1 (MQ=255)                                                gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                <  2:2269624/100‑1 (MQ=255)                                                gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                <  1:2073914/100‑1 (MQ=255)                                                    cGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGTAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                            <  1:2032427/100‑1 (MQ=255)                                                    cGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                            <  1:693282/100‑1 (MQ=255)                                                     ggCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                           <  1:708212/100‑1 (MQ=255)                                                     ggCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                           <  1:2457193/100‑1 (MQ=255)                                                     ggCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                           >  2:1572397/1‑100 (MQ=255)                                                     ggCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                           <  2:1545883/100‑1 (MQ=255)                                                       cATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgan                                         <  2:1415913/100‑2 (MQ=255)                                                        aTCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgaga                                        <  1:1507884/100‑1 (MQ=255)                                                         tCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGAc                                       >  1:2116213/1‑100 (MQ=255)                                                          cAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACCTCTGGGGATCGAGACg                                      >  1:1492374/1‑100 (MQ=255)                                                          cAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACg                                      >  1:888363/1‑100 (MQ=255)                                                           aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACgg                                     >  1:2405985/1‑100 (MQ=255)                                                           aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACgg                                     >  2:2184130/1‑100 (MQ=255)                                                           aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAACCCGATCCGCCACAGTATATGGAACATCTGGGGATCGAGACg                                      >  2:2102241/1‑99 (MQ=255)                                                           aaaaaGTGTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACgg                                     >  2:2510619/1‑100 (MQ=255)                                                            aaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCGACAGTATATGGAACATCTGGGGAGCGAGACGGt                                    >  2:922997/1‑100 (MQ=255)                                                             aaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGAGATCggaagagcg                                   >  2:1267124/1‑92 (MQ=255)                                                              aaGTTTATGGTTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                  <  2:1683196/100‑1 (MQ=255)                                                              aaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                  <  1:2102762/100‑1 (MQ=255)                                                               aGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCt                                 <  1:1263985/100‑1 (MQ=255)                                                                  ttATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGt                              >  1:2116900/1‑100 (MQ=255)                                                                   tATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                             <  2:274029/100‑1 (MQ=255)                                                                     tGGCTGAAAAAGGGATCCCGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGa                           <  2:1289390/100‑1 (MQ=255)                                                                     tGGCTGAAAAAGGGATCCAGCGGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATTGAACATCAGGGGATCGAGACGGTGATGATCGa                           >  2:621653/1‑100 (MQ=255)                                                                     tGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGa                           >  2:382651/1‑100 (MQ=255)                                                                      ggCTGAAAAAGGGATCCATCCGGATCTGTTCTACAGCGCGGAAGAAGCCGACGCGCCACAGTATACGGAACATCCGGGGATCGAGACGGCGCTGATCg                            >  2:570027/1‑98 (MQ=255)                                                                       gCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATc                         >  1:2262821/1‑100 (MQ=255)                                                                       gCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATc                         >  2:2187439/1‑100 (MQ=255)                                                                       gCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATc                         >  2:989226/1‑100 (MQ=255)                                                                       gCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGCGCTGGTCg                            >  2:1407163/1‑97 (MQ=255)                                                                        cTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                        >  2:374001/1‑100 (MQ=255)                                                                        cTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                        <  1:1250977/100‑1 (MQ=255)                                                                         tGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCg                       >  1:580187/1‑100 (MQ=255)                                                                         tGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCg                       <  2:2211890/100‑1 (MQ=255)                                                                          gAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGa                      >  1:1029731/1‑100 (MQ=255)                                                                           aaaaaGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaa                     <  1:2442538/100‑1 (MQ=255)                                                                              aaGGGATCCAGCCGGATCTATTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGCTCGAGACGGTGCTGGTCGATCCGAAAc                   >  2:2325980/1‑99 (MQ=255)                                                                                gggATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTa                >  2:1777751/1‑100 (MQ=255)                                                                                 ggATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTAc               >  2:1932794/1‑100 (MQ=255)                                                                                   aTCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCt             <  1:373736/100‑1 (MQ=255)                                                                                   aTCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCt             >  1:2421674/1‑100 (MQ=255)                                                                                     ccAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCttt           >  1:1828710/1‑100 (MQ=255)                                                                                      cAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTa          >  2:2104216/1‑100 (MQ=255)                                                                                      cAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTa          >  1:759768/1‑100 (MQ=255)                                                                                        gCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATg        >  2:2406877/1‑100 (MQ=255)                                                                                         ccGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGa       <  1:1441375/100‑1 (MQ=255)                                                                                         ccGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGa       <  1:661933/100‑1 (MQ=255)                                                                                          cGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg      <  1:1932808/100‑1 (MQ=255)                                                                                          cGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg      <  1:472106/100‑1 (MQ=255)                                                                                           ggATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGt     <  2:1300521/100‑1 (MQ=255)                                                                                              tctgttctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATc  <  2:2066511/100‑1 (MQ=255)                                                                                                  |                                                                                                CAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATC  >  NC_000913/4627740‑4627932

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