Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,768,131	A→T	pseudogene (911/945 nt)	rhsJ →	Rhs‑family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,768,131	0	A	T	97.8%	309.1 / NA	92	pseudogene (911/945 nt)	rhsJ	Rhs‑family protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (49/40);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTAC  >  NC_000913/3767984‑3768277                                                                                                                                                    |                                                                                                                                                   gAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAg                                                                                                                                                  <  2:844205/150‑1 (MQ=255) gAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAg                                                                                                                                                  <  2:890838/150‑1 (MQ=255)  aGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAgg                                                                                                                                                 <  1:785912/150‑1 (MQ=255)  aGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAgg                                                                                                                                                 <  2:924401/150‑1 (MQ=255)         aaTGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCa                                                                                                                                          >  4:26719/1‑150 (MQ=255)           tGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGa                                                                                                                                        >  4:175774/1‑150 (MQ=255)            ggTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaa                                                                                                                                       >  4:162015/1‑150 (MQ=255)            ggTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaa                                                                                                                                       >  1:1106659/1‑150 (MQ=255)             gTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaa                                                                                                                                      >  2:513964/1‑150 (MQ=255)              tAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGATCAGTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaaa                                                                                                                                     >  1:315436/1‑150 (MQ=255)                 gTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAcc                                                                                                                                  <  4:388031/150‑1 (MQ=255)                 gTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAcc                                                                                                                                  <  2:283695/150‑1 (MQ=255)                  ttAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAccc                                                                                                                                 >  2:939534/1‑150 (MQ=255)                   tAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCa                                                                                                                                <  2:266945/150‑1 (MQ=255)                    aaaaaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAg                                                                                                                               <  2:698044/150‑1 (MQ=255)                       aaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGa                                                                                                                            <  1:939534/150‑1 (MQ=255)                       aaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGa                                                                                                                            <  1:530703/150‑1 (MQ=255)                       aaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGa                                                                                                                            <  1:430246/150‑1 (MQ=255)                        atagATACAGTGCTAACAGCACCAACGAAACTTAATTCAAAAAAGGTTCTAATGGCTATAATAACAAACGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGat                                                                                                                           <  4:50346/147‑1 (MQ=255)                             tACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAg                                                                                                                      <  1:240175/150‑1 (MQ=255)                             tACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAg                                                                                                                      >  4:180574/1‑150 (MQ=255)                                aGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATa                                                                                                                   >  3:448265/1‑150 (MQ=255)                                aGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATa                                                                                                                   >  1:806274/1‑150 (MQ=255)                                   gCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGt                                                                                                                <  2:771637/150‑1 (MQ=255)                                    cTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATCACAACAAAAGACTTTGGTAATTGGATAACATATGTGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCAAGTGATATAAGATAGAtt                                                                                                               <  4:21006/150‑1 (MQ=255)                                     tAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                              <  1:513964/150‑1 (MQ=255)                                      aaCAGCACCAACGAACCTTAATTCAAATAAGGATCTACTGGCTATAACAACAAAAGACCTCGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGACAGGTCGt                                                                                                             >  2:728700/1‑150 (MQ=255)                                      aaCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCCCCCCATTTCTAATAGGCCAGGTACGCCACAACAACCCAGTGATATAAGCTAgctcgc                                                                                                             >  2:728784/1‑145 (MQ=255)                                       aCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTa                                                                                                            >  1:881884/1‑150 (MQ=255)                                       aCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTa                                                                                                            >  2:737123/1‑150 (MQ=255)                                       aCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTa                                                                                                            >  4:337367/1‑150 (MQ=255)                                           cACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATg                                                                                                        >  1:999134/1‑150 (MQ=255)                                               aaCGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATg                                                                                                    >  2:425044/1‑150 (MQ=255)                                                  gAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGaaa                                                                                                 <  3:318349/150‑1 (MQ=255)                                                      cTTAATTCAAATAGGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCAt                                                                                             <  2:11231/150‑1 (MQ=255)                                                         aaTACAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCt                                                                                          <  1:66170/150‑1 (MQ=255)                                                          aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                         >  1:112701/1‑150 (MQ=255)                                                          aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                         >  2:852118/1‑150 (MQ=255)                                                          aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                         >  2:26901/1‑150 (MQ=255)                                                           ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTg                                                                                        >  1:294098/1‑150 (MQ=255)                                                           ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTg                                                                                        >  3:511261/1‑150 (MQ=255)                                                            tCAAATAAGGATCTAATGGCTATAAAAACAAAAAACTTTGGAAAGTGGATGAGATAAGGGACATCTGCGCACCCCATTTCTTTTATTCAAAAAACGAAAAAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGt                                                                                       <  4:75778/150‑1 (MQ=255)                                                              aaaTAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTa                                                                                     >  4:67004/1‑150 (MQ=255)                                                                aTAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACt                                                                                   <  4:55467/150‑1 (MQ=255)                                                                aTAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATAAGGGCCATCTGGGCACCCCATTTCTAATATGCAAAAAAAAAAAAAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACt                                                                                   <  4:55458/150‑1 (MQ=255)                                                                  aaGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                         >  1:1037537/1‑142 (MQ=255)                                                                  aaGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                         <  2:1037537/142‑1 (MQ=255)                                                                  aaGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAt                                                                                 >  1:1090465/1‑150 (MQ=255)                                                                      aTCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCGGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttt                                                                             <  1:356666/150‑1 (MQ=255)                                                                          aaTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaaaa                                                                                                                                    >  3:174634/1‑91 (MQ=255)                                                                          aaTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaaaa                                                                                                                                    <  4:174634/91‑1 (MQ=255)                                                                           aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAtt                                                                        >  3:244241/1‑150 (MQ=255)                                                                            tGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAttt                                                                       <  2:66038/150‑1 (MQ=255)                                                                            tGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAAGAAAAAAAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAttt                                                                       <  2:315436/150‑1 (MQ=255)                                                                                   aacaacAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtttt                                                                >  2:1152665/1‑150 (MQ=255)                                                                                   aacaacAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtttt                                                                >  2:952572/1‑150 (MQ=255)                                                                                     caacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAt                                                              >  4:147444/1‑150 (MQ=255)                                                                                      aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATACGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAtt                                                             >  2:104750/1‑150 (MQ=255)                                                                                       acaaAAGACCTTTGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAttt                                                            >  2:371817/1‑150 (MQ=255)                                                                                        caaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTAt                                                                                                     >  4:270264/1‑108 (MQ=255)                                                                                        caaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTAt                                                                                                     <  3:270264/108‑1 (MQ=255)                                                                                         aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAttt                                                                       >  1:627182/1‑137 (MQ=255)                                                                                         aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAttt                                                                       <  2:627182/137‑1 (MQ=255)                                                                                         aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGt                                                          <  3:175774/150‑1 (MQ=255)                                                                                         aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGt                                                          >  2:837501/1‑150 (MQ=255)                                                                                            aGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCaa                                                       >  2:781964/1‑150 (MQ=255)                                                                                            aGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCaa                                                       >  2:997983/1‑150 (MQ=255)                                                                                             gACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTTTCTTTTATTTGTCAAt                                                      >  4:407361/1‑150 (MQ=255)                                                                                                cTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGc                                                   >  2:2651/1‑150 (MQ=255)                                                                                                 ttGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTgaa                                                  >  4:305384/1‑148 (MQ=255)                                                                                                 ttGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCa                                                  >  3:190114/1‑150 (MQ=255)                                                                                                 ttGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCa                                                  >  1:926779/1‑150 (MQ=255)                                                                                                   ggTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATg                                                <  3:469132/150‑1 (MQ=255)                                                                                                       aTTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGt                                                                                       <  2:800821/107‑1 (MQ=255)                                                                                                       aTTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGt                                                                                       >  1:800821/1‑107 (MQ=255)                                                                                                                    tGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATa                               <  3:67004/150‑1 (MQ=255)                                                                                                                       gCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaa                            <  2:1123105/150‑1 (MQ=255)                                                                                                                       gCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCGTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaa                            <  2:527776/150‑1 (MQ=255)                                                                                                                         cATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGAGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaaca                          >  1:443089/1‑150 (MQ=255)                                                                                                                           tCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAc                                                                                                                                   >  4:449931/1‑43 (MQ=255)                                                                                                                           tCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAc                                                                                                                                   <  3:449931/43‑1 (MQ=255)                                                                                                                            cTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaacaaaa                       <  4:202362/150‑1 (MQ=255)                                                                                                                            cTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaacaaaa                       <  2:999134/150‑1 (MQ=255)                                                                                                                             tGGGCACCCCATTTCTAATATGCTCGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaacaaaac                      >  4:118981/1‑150 (MQ=255)                                                                                                                             tGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAAATTCTTGGTGCTGATGTGAAAATaaaacaaaac                      >  1:833925/1‑150 (MQ=255)                                                                                                                                 cACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGaaa                  <  2:806274/150‑1 (MQ=255)                                                                                                                                          tCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCGTGCTTGTTACTGTTTTAGTAGTGCTCTTTTATTTTGCAATTGCATTGTGCTGATGGGAAAATAAAACAAAACGAAATGATGATgg         >  1:581767/1‑150 (MQ=255)                                                                                                                                           cTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATggg        <  1:93399/150‑1 (MQ=255)                                                                                                                                            tAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAttt                                                            >  2:202362/1‑97 (MQ=255)                                                                                                                                            tAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAttt                                                            <  1:202362/97‑1 (MQ=255)                                                                                                                                                atGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTa   <  3:26719/150‑1 (MQ=255)                                                                                                                                                 tGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTAc  >  3:291418/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAGGAACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTAC  >  NC_000913/3767984‑3768277

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