Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,393 A→T 6.0% intergenic (‑29/‑175) 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,3930AT 6.0% 457.3 / 13.8 134intergenic (‑29/‑175)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (62/63);  new base T (4/4);  total (66/68)
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.66e-01

Only 100 of 143 total aligned reads displayed.

CAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGG  >  NC_000913/866245‑866541
                                                                                                                                                    |                                                                                                                                                    
gaGAATACGACCAACACACTGTACCAGTGGCATCGTTTTCTGGTGGGTCTGTGGCGATACGCGAGAACGCATCTCATTAAGCGAGGTGTCGAGCGACATTAATCCGGTGGTAAATTCCATGCAAACACTACCGCGGAGGCAAAATCGaa                                                                                                                                                      <  4:410863/148‑1 (MQ=255)
cAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                                                                      <  3:192476/149‑1 (MQ=255)
cAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                                                                      <  3:356965/149‑1 (MQ=255)
cAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                                                                      <  1:315832/149‑1 (MQ=255)
     tACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGc                                                                                                                                                 <  3:252320/149‑1 (MQ=255)
     tACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                                 <  3:298650/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)
           cAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatta                                                                                                                                           >  1:1231287/1‑149 (MQ=255)
            aaaaCACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtattat                                                                                                                                          >  1:143641/1‑149 (MQ=255)
              aaCACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGt                                                                                                                                        >  4:178651/1‑149 (MQ=255)
                cacTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                                                                      <  4:275944/149‑1 (MQ=255)
                       ccAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                                                                               >  3:847746/1‑149 (MQ=255)
                         aGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                                                                             >  1:677738/1‑149 (MQ=255)
                          gTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGcc                                                                                                                            <  3:881188/149‑1 (MQ=255)
                               aGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                                                                                  >  2:1438440/1‑138 (MQ=255)
                               aGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                                                                                  <  1:1438440/138‑1 (MQ=255)
                                 cGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACt                                                                                                                     <  2:351502/149‑1 (MQ=255)
                                  gTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                                                                                    >  1:380825/1‑149 (MQ=255)
                                   tttCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                                                                      >  2:1154285/1‑130 (MQ=255)
                                    ttCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTa                                                                                                                  <  2:616844/149‑1 (MQ=255)
                                     tCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTAt                                                                                                                 <  2:677738/149‑1 (MQ=255)
                                       cTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGACAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                                                                               >  1:463144/1‑133 (MQ=255)
                                       cTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGc                                                                                                               >  1:384786/1‑149 (MQ=255)
                                         ggacGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGCATGCCa                                                                                                             >  3:1045393/4‑149 (MQ=255)
                                          ggCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCAc                                                                                                            >  3:939855/1‑149 (MQ=255)
                                               tCAGTGGGGTGACGCGAGAAAGCATTTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATCTGCCTATTATGTCAGAAAAACGCCACAGACTGTa                                                                                                                  <  3:1046668/138‑1 (MQ=255)
                                               tCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTa                                                                                                                  >  4:1046668/1‑138 (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)
                                                           cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                                                                           <  3:761630/149‑1 (MQ=255)
                                                                     cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGACtt                                                                                 <  4:647587/149‑1 (MQ=255)
                                                                      aTCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatt                                                                                                                                            >  2:2323/1‑89 (MQ=255)
                                                                      aTCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatt                                                                                                                                            <  1:2323/89‑1 (MQ=255)
                                                                       tctcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                                                                               <  3:873514/149‑1 (MQ=255)
                                                                         tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACt                                                                                                                     >  4:926951/1‑109 (MQ=255)
                                                                         tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACt                                                                                                                     <  3:926951/109‑1 (MQ=255)
                                                                         tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                                                                             <  4:382529/149‑1 (MQ=255)
                                                                           aTTATGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                                                                           >  1:1354122/1‑149 (MQ=255)
                                                                           aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                                 <  2:283706/79‑1 (MQ=255)
                                                                           aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATTCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                                                                             <  3:1058812/147‑1 (MQ=255)
                                                                           aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGTCTTTACATg                                                                           >  3:1184844/1‑149 (MQ=255)
                                                                           aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCttt                                                                                <  4:686237/144‑1 (MQ=255)
                                                                           aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCttt                                                                                >  3:686237/1‑144 (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                                                     >  2:1500914/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)
                                                                                                    aaTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaa                                                  >  1:586217/1‑149 (MQ=255)
                                                                                                     aTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaa                                                 >  3:652329/1‑149 (MQ=255)
                                                                                                      tCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaaa                                                >  2:202699/1‑149 (MQ=255)
                                                                                                      tCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCaaaaa                                                <  2:1354122/149‑1 (MQ=255)
                                                                                                        cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAAtt                                              <  2:738229/149‑1 (MQ=255)
                                                                                                          gtggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGa                                            <  3:391926/149‑1 (MQ=255)
                                                                                                           tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGaa                                           <  3:1511842/149‑1 (MQ=255)
                                                                                                             gtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATg                                         >  2:558420/1‑149 (MQ=255)
                                                                                                              tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGACTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCAGCTTTCTATATTCAAAAATTGAATGa                                        <  4:744013/149‑1 (MQ=255)
                                                                                                              tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                        <  1:632545/149‑1 (MQ=255)
                                                                                                              tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                        <  1:1191080/149‑1 (MQ=255)
                                                                                                              tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                        <  2:853114/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)
                                                                                                                       tGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                                                                       >  1:991152/1‑77 (MQ=255)
                                                                                                                       tGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCAt                               >  4:697561/1‑149 (MQ=255)
                                                                                                                           aaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATaaaa                           >  4:891047/1‑149 (MQ=255)
                                                                                                                                tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGTCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCt                      <  4:1184844/149‑1 (MQ=255)
                                                                                                                                tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                                                                     <  2:635404/102‑1 (MQ=255)
                                                                                                                                tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                                                                     >  1:635404/1‑102 (MQ=255)
                                                                                                                                   tGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGat                   >  3:375939/1‑149 (MQ=255)
                                                                                                                                    taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                                    >  2:392806/5‑131 (MQ=255)
                                                                                                                                    taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTaa                                    >  2:331823/5‑131 (MQ=255)
                                                                                                                                      ggAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATAtt                <  3:178651/149‑1 (MQ=255)
                                                                                                                                        aGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTta              >  3:1178330/1‑149 (MQ=255)
                                                                                                                                        aGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTta              >  4:1190559/1‑149 (MQ=255)
                                                                                                                                         ggCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTtat             >  1:279537/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)
                                                                                                                                                tCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCaaaa      <  3:468649/149‑1 (MQ=255)
                                                                                                                                                tCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCaaaa      >  1:315732/1‑149 (MQ=255)
                                                                                                                                                    aTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                        >  3:952534/1‑111 (MQ=255)
                                                                                                                                                    aTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGa                                        <  4:952534/111‑1 (MQ=255)
                                                                                                                                                    aTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTgg  >  2:554905/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGG  >  NC_000913/866245‑866541

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: