Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,197,365	G→A	100%	pseudogene (892/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,197,365	0	G	A	98.0%	310.1 / ‑3.0	98	pseudogene (892/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 G (1/1);  new base A (57/39);  total (58/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAC  >  NC_000913/2197269‑2197461                                                                                                 |                                                                                                 tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                               >  1:1854282/1‑100 (MQ=255) tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                               <  1:179519/100‑1 (MQ=255)  agttGGAGCAATATTACAGATTTGAATGACGCCTACGTAAACAAATATCACTGACGGAAGAGTCTCTCTGGAAAAAAAAAAAAAAAAAACTTATTACCGc                                                                                              <  2:1471218/98‑1 (MQ=255)  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAGCTTATTACCGc                                                                                              <  2:1309747/100‑1 (MQ=255)  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                               >  2:2418919/1‑99 (MQ=255)  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                              >  2:1388666/1‑100 (MQ=255)  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                              <  2:814609/100‑1 (MQ=255)  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                              >  1:671868/1‑100 (MQ=255)   cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCa                                                                                             >  2:528402/1‑100 (MQ=255)     tGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAAAATTTATCACTGACGGAAGAGTCTCTCTAGCAAAAATGTGCACAAAAACTTATTACCGCAAt                                                                                           >  2:1803623/1‑100 (MQ=255)     gggAGCAATATTCCAGATTTGAATGACGCCTACGTAACAATTTATCACTGACGGAAGAGTATCTCTGGCAAAAATGTGCACAAAAACATATTACCGCaa                                                                                            >  2:637647/2‑99 (MQ=255)      ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                          <  1:851157/100‑1 (MQ=255)       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                         >  1:461392/1‑100 (MQ=255)       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                         >  2:1051959/1‑100 (MQ=255)       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                         >  1:877282/1‑100 (MQ=255)       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGAAATTc                                                                                         >  1:473171/1‑100 (MQ=255)         gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCca                                                                                       <  1:753791/100‑1 (MQ=255)          cAATATTTCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGTAAAAATGTGCACAAAAACTTATTACCTCAATTcc                                                                                        <  2:2128508/98‑1 (MQ=255)          cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCCCTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcac                                                                                      <  2:702454/100‑1 (MQ=255)            atatTCCAGATTTGGATGAAGCCTACGCAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCAGAAAACCTTATTACCACAATTCCACat                                                                                    >  2:1023602/1‑100 (MQ=255)              atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACatat                                                                                  <  1:621990/100‑1 (MQ=255)               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                                 >  1:1293612/1‑100 (MQ=255)               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                                 >  1:1273974/1‑100 (MQ=255)               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                                 >  1:1008681/1‑100 (MQ=255)                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                               >  1:1578561/1‑100 (MQ=255)                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                              >  2:1107464/1‑100 (MQ=255)                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                              >  1:338367/1‑100 (MQ=255)                    gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGa                                                                            >  2:2051543/1‑100 (MQ=255)                    gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGa                                                                            >  1:595729/1‑100 (MQ=255)                    gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGa                                                                            >  1:531165/1‑100 (MQ=255)                     aTTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGaa                                                                           >  2:1757544/1‑100 (MQ=255)                         gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                       <  2:1008673/100‑1 (MQ=255)                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                      <  1:671419/100‑1 (MQ=255)                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                       >  1:680827/1‑99 (MQ=255)                             gACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCa                                                                   >  2:169877/1‑100 (MQ=255)                               cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                 >  1:878758/1‑100 (MQ=255)                               cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                 >  2:568435/1‑100 (MQ=255)                               agCCTACGGAAACAATTATCACTGACGGAAGAGTCTCTCTGGGAAACATGTGCACAAAAATTTATTAGCGCAATTCAACATATTTCAGAATgacggaca                                                                  >  2:618437/2‑96 (MQ=255)                                gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                                <  1:1052904/100‑1 (MQ=255)                                gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                                >  1:2315269/1‑100 (MQ=255)                                gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                                <  2:146369/100‑1 (MQ=255)                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                               <  2:1331238/100‑1 (MQ=255)                                      gTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                          >  1:251208/1‑100 (MQ=255)                                      gTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                          <  2:1446869/100‑1 (MQ=255)                                        aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                        <  1:528409/100‑1 (MQ=255)                                        aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                        >  1:529420/1‑100 (MQ=255)                                        aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                        >  1:1364343/1‑100 (MQ=255)                                         aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACACAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                       >  2:1418479/1‑100 (MQ=255)                                         aaCATTTATCACTGACGGAAGAGTCGCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                       <  2:159571/100‑1 (MQ=255)                                             tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                   <  2:1265023/100‑1 (MQ=255)                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                                 >  2:864585/1‑100 (MQ=255)                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                                 >  1:406489/1‑100 (MQ=255)                                                aTCACTGACGGAAGTGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                                <  2:330081/100‑1 (MQ=255)                                                aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                                <  2:1168453/100‑1 (MQ=255)                                                  cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCAACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                              >  1:2228525/1‑100 (MQ=255)                                                   aCTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatta                                             <  2:2253537/100‑1 (MQ=255)                                                    cTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                            <  2:448585/100‑1 (MQ=255)                                                      gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                          <  1:2089670/100‑1 (MQ=255)                                                      gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                           >  1:323749/1‑99 (MQ=255)                                                       aCGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCg                                         <  2:682967/100‑1 (MQ=255)                                                        cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGa                                        >  1:1519718/1‑100 (MQ=255)                                                        cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGa                                        <  1:611568/100‑1 (MQ=255)                                                        cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGa                                        <  1:1909118/100‑1 (MQ=255)                                                        cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGa                                        <  1:215246/100‑1 (MQ=255)                                                         ggAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATGACCCGaa                                       >  1:423421/1‑100 (MQ=255)                                                           aaGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaa                                     <  2:262052/100‑1 (MQ=255)                                                              agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                  <  2:1930707/100‑1 (MQ=255)                                                              agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                  <  1:239205/100‑1 (MQ=255)                                                               gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCa                                 >  1:726452/1‑100 (MQ=255)                                                               gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCa                                 >  1:1582712/1‑100 (MQ=255)                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                               >  1:1934590/1‑100 (MQ=255)                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                               >  2:686878/1‑100 (MQ=255)                                                                  tctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                              <  1:2060265/100‑1 (MQ=255)                                                                     ctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAtt                           >  2:382869/1‑100 (MQ=255)                                                                       ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                         >  1:2075589/1‑100 (MQ=255)                                                                       ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                         <  1:167471/100‑1 (MQ=255)                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        <  2:1273966/100‑1 (MQ=255)                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        >  1:514375/1‑100 (MQ=255)                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                       >  2:2553821/1‑100 (MQ=255)                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                       <  2:1809057/100‑1 (MQ=255)                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                       >  1:1831290/1‑100 (MQ=255)                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTAATTACCCGAAAAACCAGAAATTGca                        >  1:1671728/1‑99 (MQ=255)                                                                          aaaaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcaca                      <  2:15507/100‑1 (MQ=255)                                                                           aaaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACAt                     >  1:1212284/1‑100 (MQ=255)                                                                           aaaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACAt                     <  2:296909/100‑1 (MQ=255)                                                                              aTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGaa                  >  1:405257/1‑100 (MQ=255)                                                                                gtgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAt                >  1:710092/1‑100 (MQ=255)                                                                                 tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt               >  1:2205816/1‑100 (MQ=255)                                                                                 tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt               >  1:1311940/1‑100 (MQ=255)                                                                                  gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg              >  1:1707467/1‑100 (MQ=255)                                                                                   cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             >  1:1945414/1‑100 (MQ=255)                                                                                   cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             <  1:768330/100‑1 (MQ=255)                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTTTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             <  2:301462/99‑1 (MQ=255)                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc            <  1:2494643/100‑1 (MQ=255)                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc            >  1:67286/1‑100 (MQ=255)                                                                                    acaAAAACGTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc            >  1:2165730/1‑100 (MQ=255)                                                                                           cTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc     >  2:803014/1‑100 (MQ=255)                                                                                            ttattaCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCt    <  2:1016255/100‑1 (MQ=255)                                                                                            ttattaCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCt    <  1:568442/100‑1 (MQ=255)                                                                                              aTTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAc  <  1:1621022/100‑1 (MQ=255)                                                                                                 |                                                                                                 TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAC  >  NC_000913/2197269‑2197461

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