Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,388 A→T 6.0% intergenic (‑24/‑180) moeA ← / → iaaA molybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913866,3880AT 6.0% 456.7 / 14.9 133intergenic (‑24/‑180)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (59/66);  new base T (4/4);  total (63/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 139 total aligned reads displayed.

TCGCCAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAA  >  NC_000913/866241‑866536
                                                                                                                                                   |                                                                                                                                                    
tcgcCAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAAt                                                                                                                                                     >  3:1205628/1‑149 (MQ=255)
    gaGAATACGACCAACACACTGTACCAGTGGCATCGTTTTCTGGTGGGTCTGTGGCGATACGCGAGAACGCATCTCATTAAGCGAGGTGTCGAGCGACATTAATCCGGTGGTAAATTCCATGCAAACACTACCGCGGAGGCAAAATCGaa                                                                                                                                                 <  4:410863/148‑1 (MQ=255)
    cAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                                                                 <  3:356965/149‑1 (MQ=255)
    cAGAATACGACCAAAACACGGCACCAGTGGCAGCGTTTCCTGGGCGGTCAGGGGGGTGACGCGAGAAAGCATCTCATTAAGCGCAGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGac                                                                                                                                                 <  2:1013179/149‑2 (MQ=255)
      gAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAtt                                                                                                                                               >  4:1211682/1‑149 (MQ=255)
      gAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAtt                                                                                                                                               <  2:190872/149‑1 (MQ=255)
         tACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGc                                                                                                                                            <  3:252320/149‑1 (MQ=255)
           cGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCt                                                                                                                                          <  1:1480549/149‑1 (MQ=255)
            gACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtc                                                                                                                                         <  4:137434/149‑2 (MQ=255)
            gACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                                                                         >  4:160773/1‑149 (MQ=255)
             aCCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtat                                                                                                                                        <  2:249772/149‑1 (MQ=255)
               cAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatta                                                                                                                                      >  1:1231287/1‑149 (MQ=255)
                  aaCACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGt                                                                                                                                   >  4:178651/1‑149 (MQ=255)
                    cacTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                                                                 <  4:275944/149‑1 (MQ=255)
                         tACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaaa                                                                                                                            >  3:944515/1‑149 (MQ=255)
                           ccAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                                                                          <  2:572337/149‑1 (MQ=255)
                            cAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACg                                                                                                                         >  4:712670/1‑149 (MQ=255)
                             aGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                                                                        >  1:677738/1‑149 (MQ=255)
                              gTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGcc                                                                                                                       <  3:881188/149‑1 (MQ=255)
                                   aGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                                                                             <  1:1438440/138‑1 (MQ=255)
                                   aGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                                                                             >  2:1438440/1‑138 (MQ=255)
                                      gTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                                                                               >  1:380825/1‑149 (MQ=255)
                                       tttCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                                                                 >  2:1154285/1‑130 (MQ=255)
                                         tCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTAt                                                                                                            <  2:677738/149‑1 (MQ=255)
                                             ggacGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGCATGCCa                                                                                                        >  3:1045393/4‑149 (MQ=255)
                                              ggCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCAc                                                                                                       >  3:939855/1‑149 (MQ=255)
                                                   tCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                                                                  >  4:240214/1‑149 (MQ=255)
                                                        ggggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTa                                                                                             <  2:411139/149‑1 (MQ=255)
                                                        ggggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTa                                                                                             <  1:1045464/149‑1 (MQ=255)
                                                         gggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAg                                                                                            >  4:705237/1‑149 (MQ=255)
                                                           gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                            <  3:194665/99‑1 (MQ=255)
                                                               cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                                                                      <  4:939855/149‑1 (MQ=255)
                                                               cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                                                                      <  1:1371615/149‑1 (MQ=255)
                                                                  gagaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                                                                   <  2:143641/149‑1 (MQ=255)
                                                                   agaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                                                                  <  3:160773/149‑1 (MQ=255)
                                                                         cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGACtt                                                                            <  4:647587/149‑1 (MQ=255)
                                                                           tctcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                                                                          <  3:873514/149‑1 (MQ=255)
                                                                             tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                                                                        <  4:382529/149‑1 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                            <  2:283706/79‑1 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                            >  1:283706/1‑79 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATTCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                                                                        <  3:1058812/147‑1 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATTCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                                                                        >  4:1058812/1‑147 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGTCTTTACATg                                                                      >  3:1184844/1‑149 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCttt                                                                           <  4:686237/144‑1 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCttt                                                                           >  3:686237/1‑144 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                                                                      >  3:778109/1‑149 (MQ=255)
                                                                               aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                                                                      >  1:760594/1‑149 (MQ=255)
                                                                                    gcgcGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCata                                                                 >  2:940744/1‑149 (MQ=255)
                                                                                    gcgcGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCata                                                                 <  1:520374/149‑1 (MQ=255)
                                                                                     cgcgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                                                                <  1:394956/149‑1 (MQ=255)
                                                                                         gtgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCAt                                                            <  1:929646/149‑1 (MQ=255)
                                                                                              gAGCGACATCATTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc                                                              >  3:1234442/1‑142 (MQ=255)
                                                                                              gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc                                                              <  4:1234442/142‑1 (MQ=255)
                                                                                              gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc                                                       >  4:769470/1‑149 (MQ=255)
                                                                                               aGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCt                                                      >  4:494269/1‑149 (MQ=255)
                                                                                                     aTCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTc                                                >  1:378670/1‑149 (MQ=255)
                                                                                                       cAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGACACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaa                                              >  1:1136232/1‑149 (MQ=255)
                                                                                                        aaTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                                                                   >  4:1164269/1‑111 (MQ=255)
                                                                                                        aaTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                                                                   <  3:1164269/111‑1 (MQ=255)
                                                                                                        aaTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaa                                             >  1:586217/1‑149 (MQ=255)
                                                                                                         aTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaa                                            >  3:652329/1‑149 (MQ=255)
                                                                                                         aTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaa                                            >  1:480943/1‑149 (MQ=255)
                                                                                                          tCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaaa                                           >  2:202699/1‑149 (MQ=255)
                                                                                                          tCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaaa                                           <  2:1354122/149‑1 (MQ=255)
                                                                                                            cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAAtt                                         <  2:384786/149‑1 (MQ=255)
                                                                                                              gtggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGa                                       <  3:391926/149‑1 (MQ=255)
                                                                                                               tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGaa                                      <  4:928133/149‑1 (MQ=255)
                                                                                                               tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGaa                                      <  3:715688/149‑1 (MQ=255)
                                                                                                               tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGaa                                      >  2:764301/1‑149 (MQ=255)
                                                                                                                 gtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATg                                    >  2:558420/1‑149 (MQ=255)
                                                                                                                  tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                   <  1:1191080/149‑1 (MQ=255)
                                                                                                                  tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                   <  1:1429380/149‑1 (MQ=255)
                                                                                                                  tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                   <  1:632545/149‑1 (MQ=255)
                                                                                                                    aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGt                                 <  3:33648/149‑1 (MQ=255)
                                                                                                                      ttCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               <  3:769470/149‑1 (MQ=255)
                                                                                                                      ttCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               <  1:929858/149‑1 (MQ=255)
                                                                                                                         cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc                                                       <  4:625274/122‑1 (MQ=255)
                                                                                                                          aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCa                           >  3:926416/1‑149 (MQ=255)
                                                                                                                           tGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                                                                  <  2:991152/77‑1 (MQ=255)
                                                                                                                               aaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATaaaa                      >  4:891047/1‑149 (MQ=255)
                                                                                                                                   cTCCTGCGGAGGCAAAATCGAATTAGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTc                  <  4:778109/149‑1 (MQ=255)
                                                                                                                                    tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGTCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCt                 <  4:1184844/149‑1 (MQ=255)
                                                                                                                                    tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                                                                <  2:635404/102‑1 (MQ=255)
                                                                                                                                    tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                                                                >  1:635404/1‑102 (MQ=255)
                                                                                                                                    tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCt                 >  1:187021/1‑149 (MQ=255)
                                                                                                                                    tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCt                 <  3:208873/149‑1 (MQ=255)
                                                                                                                                     acataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTg                >  2:695838/8‑149 (MQ=255)
                                                                                                                                       tGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGat              >  3:375939/1‑149 (MQ=255)
                                                                                                                                        taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               >  2:331823/5‑131 (MQ=255)
                                                                                                                                        taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               >  2:392806/5‑131 (MQ=255)
                                                                                                                                        taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               <  1:392806/127‑1 (MQ=255)
                                                                                                                                        taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                               <  1:331823/127‑1 (MQ=255)
                                                                                                                                          ggAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATAtt           <  3:178651/149‑1 (MQ=255)
                                                                                                                                             ggCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTtat        >  1:279537/1‑149 (MQ=255)
                                                                                                                                             ggCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTtat        >  1:1256668/1‑149 (MQ=255)
                                                                                                                                                aaaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGc     <  2:760594/149‑1 (MQ=255)
                                                                                                                                                 aaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGTCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCa    >  1:1041074/1‑149 (MQ=255)
                                                                                                                                                 aaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCa    >  1:718333/1‑149 (MQ=255)
                                                                                                                                                  aaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCaa   >  3:316927/1‑149 (MQ=255)
                                                                                                                                                   aTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCaaa  <  1:960304/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGCCAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAA  >  NC_000913/866241‑866536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: