Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,196,608	T→C	5.5%	pseudogene (135/2847 nt)	yehH →	DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,196,608	0	T	C	5.5%	226.5 / 3.1	73	pseudogene (135/2847 nt)	yehH	DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (2/2);  total (40/33)
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.79e-01

CAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAGGAA  >  NC_000913/2196462‑2196755                                                                                                                                                   |                                                                                                                                                    cAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGAtgctg                                                                                                                                                    <  8:343719/148‑1 (MQ=255)    ggATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTgcgg                                                                                                                                                 >  4:334330/1‑148 (MQ=255)        gCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATATATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTAGAGATGCTgcggcagcg                                                                                                                                            >  2:144762/1‑149 (MQ=255)         cAGAATGAGACACTTTATCTCTCAGGACTAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTGGGCACCAAAGGGAAAAGCCAGATAAAAAGTTTTTCAGATGCTgcggcagcg                                                                                                                                            <  3:234632/148‑1 (MQ=255)            aaTGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTgcggcagcggca                                                                                                                                         >  1:63237/1‑148 (MQ=255)                 gaCACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCg                                                                                                                                   >  6:255106/1‑149 (MQ=255)                 gaCACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCg                                                                                                                                   >  5:277318/1‑149 (MQ=255)                   cacTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGACACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGa                                                                                                                                 <  5:292833/149‑1 (MQ=255)                   cacTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGa                                                                                                                                 <  5:255106/149‑1 (MQ=255)                    acTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAg                                                                                                                                <  6:270344/149‑1 (MQ=255)                          tctatcAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTaa                                                                                                                           <  3:334330/148‑1 (MQ=255)                           ctatcAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGc                                                                                                                         >  5:30709/1‑149 (MQ=255)                            tatcAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCt                                                                                                                        >  2:68641/1‑149 (MQ=255)                             atcAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCt                                                                                                                        >  6:188187/1‑148 (MQ=255)                              tcAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCt                                                                                                                        <  2:63237/147‑1 (MQ=255)                                  gACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGc                                                                                                                  >  1:306340/1‑149 (MQ=255)                                   aCGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCg                                                                                                                 <  7:389934/149‑1 (MQ=255)                                    cGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGTCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCgg                                                                                                                <  1:144762/149‑1 (MQ=255)                                        aaaTCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGaa                                                                                                            <  1:68641/149‑1 (MQ=255)                                          aTCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAgg                                                                                                          <  5:188187/149‑1 (MQ=255)                                           tCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGt                                                                                                         <  5:189240/149‑1 (MQ=255)                                               aTAAATTCAGGGCAGTTGAGCAACAGGGAGACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATTAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTgaag                                                                                                     <  5:170103/149‑1 (MQ=255)                                               aTAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTg                                                                                                        >  4:277715/1‑146 (MQ=255)                                                 aaaTTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGt                                                                                                         <  6:169620/143‑1 (MQ=255)                                                 aaaTTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGt                                                                                                         >  5:169620/1‑143 (MQ=255)                                                  aaTTCAGGGCAGGTGCGCATCAGGGAACCGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCGGATTGCGGAGAAGGTGaagaag                                                                                                  <  2:267672/149‑1 (MQ=255)                                                    ttCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATAACAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAggg                                                                                                >  8:218530/1‑149 (MQ=255)                                                              ttGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAg                                                                                                                                >  1:143814/1‑107 (MQ=255)                                                              ttGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAg                                                                                                                                <  2:143814/107‑1 (MQ=255)                                                                   cAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGc                                                                                 >  7:106629/1‑149 (MQ=255)                                                                    aaCAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGc                                                                                 >  4:196652/1‑148 (MQ=255)                                                                       aGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAg                                                                             >  8:333884/1‑149 (MQ=255)                                                                           aaaCGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGAt                                                                          >  1:175182/1‑148 (MQ=255)                                                                              cGACTTGCATATCACGTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAACAGTTTTGCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATaaat                                                                      <  6:155844/149‑3 (MQ=255)                                                                                  ttGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTg                                                                                                        >  2:329861/1‑111 (MQ=255)                                                                                  ttGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCttt                                                                  <  6:30709/149‑1 (MQ=255)                                                                                     cATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGc                                                                >  7:283656/1‑148 (MQ=255)                                                                                        atCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAAcc                                                            >  1:388437/1‑149 (MQ=255)                                                                                                 ggAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAAc                                                   <  2:175182/149‑1 (MQ=255)                                                                                                  gAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGACGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACg                                                  <  3:277715/149‑1 (MQ=255)                                                                                                     aaGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGt                                                >  5:260782/1‑148 (MQ=255)                                                                                                       gTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAAc                                             <  8:106629/149‑1 (MQ=255)                                                                                                       gTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAAc                                             <  2:388437/149‑1 (MQ=255)                                                                                                            cACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCCGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGct                                        >  3:295858/1‑149 (MQ=255)                                                                                                            cACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCCGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGct                                        >  3:290951/1‑149 (MQ=255)                                                                                                             aCCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAg                                                                                  <  4:232229/106‑1 (MQ=255)                                                                                                             aCCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGc                                                                <  1:154409/124‑1 (MQ=255)                                                                                                             aCCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGc                                                                >  2:154409/1‑124 (MQ=255)                                                                                                             aCCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGctc                                       >  7:61369/1‑149 (MQ=255)                                                                                                             aCCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGGATGTGGAGCAAg                                                                                  >  3:232229/1‑106 (MQ=255)                                                                                                              ccAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGctc                                       >  1:333021/1‑148 (MQ=255)                                                                                                               cAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGCGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGTTAACGGGctctc                                     >  5:280419/1‑149 (MQ=255)                                                                                                                 aaGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGctctct                                    >  2:336298/1‑148 (MQ=255)                                                                                                                  aGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCtctctc                                   >  8:43207/1‑148 (MQ=255)                                                                                                                     gCAAAGCCAGATAAAAAGTTTTTCAGATGCCGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAgg                               <  4:290951/149‑1 (MQ=255)                                                                                                                     gCAAAGCCAGATAAAAAGTTTTTCAGATGCCGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAgg                               <  4:295858/149‑1 (MQ=255)                                                                                                                       aaaGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAgg                               >  2:385713/1‑147 (MQ=255)                                                                                                                             aGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGTCTCTCTCAAGGTAGCGGAc                       <  2:333021/149‑1 (MQ=255)                                                                                                                             aGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGctct                                      <  4:248844/134‑1 (MQ=255)                                                                                                                             aGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGctct                                      >  3:248844/1‑134 (MQ=255)                                                                                                                             aGATAAAAAGGTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGAc                       <  5:295547/149‑1 (MQ=255)                                                                                                                               aTAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACtt                     <  2:306340/149‑1 (MQ=255)                                                                                                                               aTAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACt                      >  6:282388/1‑148 (MQ=255)                                                                                                                                tAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTa                    <  7:218530/149‑1 (MQ=255)                                                                                                                                   aaaGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCaa                                 <  6:13726/133‑1 (MQ=255)                                                                                                                                   aaaGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCaa                                 >  5:13726/1‑133 (MQ=255)                                                                                                                                   aaaGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATc                  >  1:141097/1‑148 (MQ=255)                                                                                                                                    aaGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATcca                <  8:61369/149‑1 (MQ=255)                                                                                                                                     aGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATccac               <  3:196652/149‑1 (MQ=255)                                                                                                                                      gTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATccac               >  6:347216/1‑148 (MQ=255)                                                                                                                                             aGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCattatt        >  3:253695/1‑148 (MQ=255)                                                                                                                                             aGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCattatt        >  4:203590/1‑148 (MQ=255)                                                                                                                                              gATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCattatt        >  2:302249/1‑147 (MQ=255)                                                                                                                                               aTGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAg     >  5:215028/1‑149 (MQ=255)                                                                                                                                                 gctgcGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAgg    >  8:167281/1‑148 (MQ=255)                                                                                                                                                  ctgcGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGCAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAGGaa  <  6:277318/149‑1 (MQ=255)                                                                                                                                                  ctgcGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAGGaa  <  7:333884/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGGGAAAAGTTGGCACCAAAGGGCAAAGCCAGATAAAAAGTTTTTCAGATGCTGCGGCAGCGGCAAAAGCGGAGCTTAAGCTGATTGCGGAGAAGGTGAAGAAGGGGTATGTGGAGCAAGCGAAGGATAATTCTTTGCAACCTTCCCAAACGGTAACGGGCTCTCTCAAGGTAGCGGACTTATCCACCATTATTCAGGAA  >  NC_000913/2196462‑2196755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑