Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,889,005	T→C	G92G (GGT→GGC)	tnaA →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,889,005	0	T	C	98.7%	275.2 / ‑4.5	79	G92G (GGT→GGC)	tnaA	tryptophanase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/48);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATA  >  NC_000913/3888869‑3889141                                                                                                                                         |                                                                                                                                         ttATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTAt                                                                                                                                       <  1:1089677/140‑1 (MQ=255) ttATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTAt                                                                                                                                       <  2:408871/140‑1 (MQ=255)  tATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATc                                                                                                                                      <  2:1083205/140‑1 (MQ=255)   aTCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCa                                                                                                                                     <  2:464380/140‑1 (MQ=255)    tCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCaa                                                                                                                                    <  1:346820/140‑1 (MQ=255)             tGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAACGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCAtt                                                                                                                           <  1:590440/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGa                                                                                                                       <  1:673211/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGa                                                                                                                       <  2:360523/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGa                                                                                                                       <  2:567229/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGa                                                                                                                       <  1:689818/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGa                                                                                                                       <  1:1143989/140‑1 (MQ=255)                 cGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGa                                                                                                                       <  2:167610/140‑1 (MQ=255)                  gACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGAc                                                                                                                      >  2:893725/1‑140 (MQ=255)                    cAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTc                                                                                                                    >  1:24806/1‑140 (MQ=255)                      gCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCAc                                                                                                                  <  2:941009/140‑1 (MQ=255)                       cGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCAcc                                                                                                                 >  1:937919/1‑140 (MQ=255)                          cACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAgg                                                                                                              >  2:219369/1‑140 (MQ=255)                              ggggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCg                                                                                                          <  1:219369/140‑1 (MQ=255)                               gggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGt                                                                                                         >  2:482518/1‑140 (MQ=255)                                    gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAAGGCCGTGgcgc                                                                                                    >  2:105280/1‑140 (MQ=255)                                        cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:305380/140‑1 (MQ=255)                                        cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:1161510/140‑1 (MQ=255)                                        cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:61972/140‑1 (MQ=255)                                          cAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                >  2:1069978/1‑138 (MQ=255)                                          cAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:1069978/138‑1 (MQ=255)                                          cAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                >  2:1103220/1‑138 (MQ=255)                                          cAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:1103220/138‑1 (MQ=255)                                          cAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCa                                                                                              <  1:354061/140‑1 (MQ=255)                                           agagCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaa                                                                                             >  2:337865/1‑140 (MQ=255)                                             agCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAt                                                                                           <  2:1077658/140‑1 (MQ=255)                                                atacAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta                                                                                        <  2:420760/137‑1 (MQ=255)                                                aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta                                                                                        <  2:585988/140‑1 (MQ=255)                                                aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta                                                                                        <  1:831875/140‑1 (MQ=255)                                                aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta                                                                                        <  1:964683/140‑1 (MQ=255)                                                       cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                  <  2:593075/123‑1 (MQ=255)                                                       cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                  >  1:191063/1‑123 (MQ=255)                                                       cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                  <  2:191063/123‑1 (MQ=255)                                                       cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                  >  1:593075/1‑123 (MQ=255)                                                       cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCg                                                                                 >  1:1060260/1‑140 (MQ=255)                                                        tGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                <  1:465248/124‑1 (MQ=255)                                                        tGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                >  2:465248/1‑124 (MQ=255)                                                        tGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCgg                                                                                >  1:307734/1‑140 (MQ=255)                                                         gCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt                                                                               <  2:937919/140‑1 (MQ=255)                                                           gatgatGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTAc                                                                             <  2:117235/140‑1 (MQ=255)                                                                tgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAtt                                                                        >  2:494132/1‑140 (MQ=255)                                                                tgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAtt                                                                        >  2:361557/1‑140 (MQ=255)                                                                  cgcgGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaa                                                                      >  2:834922/1‑140 (MQ=255)                                                                      gcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTaaaaaa                                                                  >  1:534112/1‑140 (MQ=255)                                                                        gacgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGTAAATCTATATTCCGGTACTGATTAAAAAAcg                                                                <  1:2761/140‑1 (MQ=255)                                                                          cgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAAcgcg                                                              >  1:854264/1‑140 (MQ=255)                                                                             aGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGc                                                           <  1:337865/140‑1 (MQ=255)                                                                              gCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCa                                                          <  1:856638/140‑1 (MQ=255)                                                                                 tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGa                                                       >  1:952668/1‑140 (MQ=255)                                                                                   cAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaa                                                     <  2:24806/140‑1 (MQ=255)                                                                                   cAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaa                                                     <  2:106657/140‑1 (MQ=255)                                                                                   cAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaa                                                     <  1:356685/140‑1 (MQ=255)                                                                                   cAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaa                                                     <  2:474953/140‑1 (MQ=255)                                                                                         cAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAGAAACGCGAGCAGGAAAAAGGcc                                               >  2:63750/1‑140 (MQ=255)                                                                                              gTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGAt                                          <  1:910509/140‑1 (MQ=255)                                                                                              gTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGAt                                          <  1:639962/140‑1 (MQ=255)                                                                                                     ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAGAAACGCGAGCAGGAAAAAGGCCTGGATCgcagca                                   <  1:63750/140‑1 (MQ=255)                                                                                                             tAGCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggtg                           >  2:674848/1‑140 (MQ=255)                                                                                                               gCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGc                         >  2:989819/1‑140 (MQ=255)                                                                                                               gCCGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGc                         >  1:235365/1‑140 (MQ=255)                                                                                                                ccGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCg                        <  1:893725/140‑1 (MQ=255)                                                                                                                ccGAGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCg                        <  1:482518/140‑1 (MQ=255)                                                                                                                   aGTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtc                     >  1:120072/1‑140 (MQ=255)                                                                                                                    gTCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                    >  1:650400/1‑140 (MQ=255)                                                                                                                     tCAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtctc                   >  1:535208/1‑140 (MQ=255)                                                                                                                      cAGTGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTctct                  <  2:536341/140‑1 (MQ=255)                                                                                                                         tGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                    >  2:778586/1‑135 (MQ=255)                                                                                                                         tGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                    <  1:778586/135‑1 (MQ=255)                                                                                                                         tGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc               >  2:10980/1‑140 (MQ=255)                                                                                                                         tGAAAAATATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc               >  2:197115/1‑140 (MQ=255)                                                                                                                            aaaaTATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAt            <  1:494132/140‑1 (MQ=255)                                                                                                                            aaaaTATCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAt            <  1:361557/140‑1 (MQ=255)                                                                                                                                 atCTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt       <  2:449772/140‑1 (MQ=255)                                                                                                                                   cTTTGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTTTTTg     >  1:428983/1‑140 (MQ=255)                                                                                                                                    tttGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa    <  2:235365/140‑1 (MQ=255)                                                                                                                                      tGGCTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa  <  2:997351/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATA  >  NC_000913/3888869‑3889141

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