Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,415,525 G→T 9.7% E417D (GAG→GAT) ‡ aidB → DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,415,5250GT 9.7% 356.0 / 25.3 131E417D (GAG→GAT) ‡aidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (56/62);  new base T (8/5);  total (64/67)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 140 total aligned reads displayed.

GCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCGGT  >  NC_000913/4415376‑4415674
                                                                                                                                                     |                                                                                                                                                     
gCCAAATACTCCCTGAGGGCGCTATTACTACCGCCTGCCGCAAAATTTGTGATCAGCAAAAGAGGTATGCCGTTTCTGGCCGTAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGTGGATTGTGAGCTGCCGCGGCTTTACCGGgag                                                                                                                                                       <  2:334819/150‑1 (MQ=255)
gCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgag                                                                                                                                                       <  1:569900/150‑1 (MQ=255)
 ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTATGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                      >  2:1062154/1‑150 (MQ=255)
  cAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGAt                                                                                                                                                     >  1:1018328/1‑150 (MQ=255)
  cAAAGAAGCCCTGCGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGAt                                                                                                                                                     <  1:1307652/150‑1 (MQ=255)
    aagaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGc                                                                                                                                                   >  1:746294/1‑150 (MQ=255)
    aagaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGc                                                                                                                                                   >  1:168102/1‑150 (MQ=255)
       aagCCCTGTGGGCGCGTTTATTTACGCCGGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCgg                                                                                                                                                <  1:857017/150‑1 (MQ=255)
        atcccTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGGGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGt                                                                                                                                               <  1:925215/148‑1 (MQ=255)
         gCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTa                                                                                                                                              >  2:690989/1‑150 (MQ=255)
         gCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTa                                                                                                                                              <  2:1758996/150‑1 (MQ=255)
            cTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaag                                                                                                                                           <  1:1414547/150‑2 (MQ=255)
             tgtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                          <  1:1642673/150‑3 (MQ=255)
             tgtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGGTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                          <  1:667606/150‑3 (MQ=255)
             tgtgGGCGCGTCTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGGGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                          <  1:450358/150‑3 (MQ=255)
              gtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAg                                                                                                                                         <  1:1024666/150‑1 (MQ=255)
                gggCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCTGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTa                                                                                                                                       >  2:956551/1‑150 (MQ=255)
                gggCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTTCCGGGAGATGCCGGTAAACAGTa                                                                                                                                       >  2:341192/1‑150 (MQ=255)
                   cgcgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTAttt                                                                                                                                    <  1:382672/150‑1 (MQ=255)
                     cgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGCATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTgg                                                                                                                                  <  1:1251831/150‑1 (MQ=255)
                       tttatttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGCAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGa                                                                                                                                >  1:1668279/1‑150 (MQ=255)
                       tttatttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGa                                                                                                                                >  2:49263/1‑150 (MQ=255)
                       attatttaCGCCTGCGGCAAAATTTGTGATCTGCAATAGACATATGCACTTTGAGGACCAAGCGATGGAGGTGCTGGGCGGCATTCGTTTTTGCGAGGAGAGACAGCTGCCGCGACTTTACCGGGAGATGCCGGTAAACAGTATTTgtca                                                                                                                                >  1:320967/2‑147 (MQ=255)
                          atttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAgg                                                                                                                             <  1:1875660/150‑1 (MQ=255)
                          atttaCGCCTGCGGCAAAATTTCTGATCTCCAAAAGAGGTATGCCGCTTGTGGCCGAAGCGAGGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAgg                                                                                                                             <  2:252138/150‑1 (MQ=255)
                               cGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAACCAGGATTTGGGATGGTTCCg                                                                                                                        >  2:86968/1‑150 (MQ=255)
                                cccTGCGGCAAAATTTGTGATCTTCAAAAGCGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCGGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCgg                                                                                                                       <  2:648397/149‑1 (MQ=255)
                                cacTGCGGAAAAATTTTGTATCTGCAAAAGAGGTATGACGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGTGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGTGAAGGTTCCgg                                                                                                                       <  2:655053/148‑1 (MQ=255)
                                    gCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAat                                                                                                                   <  2:1120096/150‑1 (MQ=255)
                                    gCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAat                                                                                                                   <  1:1286694/150‑1 (MQ=255)
                                     cGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAata                                                                                                                  >  2:389269/1‑150 (MQ=255)
                                        cAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGAAGGAGGAGCTGGGCGCCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCTGTAAACAGTATTTGGGAAGGTTCCGGCAAtatta                                                                                                               >  1:1213999/1‑150 (MQ=255)
                                           aaTTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                            >  2:1565215/1‑150 (MQ=255)
                                           aaTTTGTGATATGCAAAAGAGGTATGCCGTTTGTGGCCGTAGCGTTGGAGGTGCTGTGCTGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGTCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                            <  1:322077/150‑1 (MQ=255)
                                             tttGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGc                                                                                                          <  2:1455280/150‑1 (MQ=255)
                                                gtgATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTg                                                                                                       >  1:398977/1‑150 (MQ=255)
                                                  gATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGTGAGTGAGCTGCAGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGAAATATTATGTGCCTGGa                                                                                                     >  1:81442/1‑150 (MQ=255)
                                                    tCTGCAAAAGAGGGATGCCGTTTGTGGCCGAAGCGATGGCGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtg                                                                                                   <  2:1286897/150‑1 (MQ=255)
                                                      tGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAATCAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtgtg                                                                                                 >  2:186066/1‑150 (MQ=255)
                                                      tGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtgtg                                                                                                 >  1:1593015/1‑150 (MQ=255)
                                                         aaaaGAGGTATGCCGTTTGTGGCCGAAGCCATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACAGGGAGATGCCGGTAAACAGTATTTGGGATGGTTCCGGCAATATTATGTGCCTGGATGTGTTg                                                                                              >  1:150803/1‑150 (MQ=255)
                                                          aaaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCTAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTgc                                                                                             >  2:96208/1‑150 (MQ=255)
                                                           aaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTgcg                                                                                            <  1:389269/150‑1 (MQ=255)
                                                             gagGTATGCCGTTTGAGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGGGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTAATTGGGAAGGTTCCGGCAAAATTATGTGCCTGGATGTGTTGcgcg                                                                                          <  2:1643992/150‑1 (MQ=255)
                                                                gTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtc                                                                                       <  2:973888/150‑1 (MQ=255)
                                                                 tATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtct                                                                                      >  1:1186/1‑150 (MQ=255)
                                                                  ctGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtctc                                                                                     <  2:1838279/149‑1 (MQ=255)
                                                                   tgtcGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                    <  1:14252/147‑1 (MQ=255)
                                                                   tGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTATACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                    <  1:1045643/150‑1 (MQ=255)
                                                                   tGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                    >  1:825660/1‑150 (MQ=255)
                                                                     ccGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGCGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAt                                                                                  >  2:1190293/1‑150 (MQ=255)
                                                                         ttGTGGCCGCTGCGATCGAGGTGCTGGGCGGGATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCAGGTAAACATTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagc                                                                              <  2:1251760/150‑1 (MQ=255)
                                                                          tgtgGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                             >  2:1722523/1‑150 (MQ=255)
                                                                             ggCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagcaagc                                                                          >  2:1646332/1‑150 (MQ=255)
                                                                              gCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCg                                                                         <  1:1483098/150‑1 (MQ=255)
                                                                               ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATACGCAAGCgg                                                                        >  1:344564/1‑150 (MQ=255)
                                                                               ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                        >  1:601879/1‑150 (MQ=255)
                                                                               ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                        >  1:1269667/1‑150 (MQ=255)
                                                                                 gAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGGGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGc                                                                      >  2:666371/1‑150 (MQ=255)
                                                                                 gAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAGAAGCAAGCGGGc                                                                      >  1:1169340/1‑150 (MQ=255)
                                                                                  aaGCGAGGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGGCTCAATAAGCAAGCGGGCg                                                                     >  2:1479964/1‑150 (MQ=255)
                                                                                       aTGGAGGTGCTGGGCGGCATTGGGTATTGCGAGGAGAGGGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGGTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATAc                                                                <  2:755070/150‑1 (MQ=255)
                                                                                         ggAGGTGCAGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGa                                                              >  2:1853149/1‑150 (MQ=255)
                                                                                            ggTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACtt                                                           <  2:932799/150‑1 (MQ=255)
                                                                                             gTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTa                                                          <  1:607299/150‑1 (MQ=255)
                                                                                              tGCTGGGCGGCATTGGTTATTGCGAGGAGTGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTAt                                                         <  2:1513151/150‑1 (MQ=255)
                                                                                               gCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTAtt                                                        >  2:1074577/1‑150 (MQ=255)
                                                                                                ccgggCGGCATTGGTTATTGCGAGGAGAGTGAGCAGCCGCGGCTTTACCTTTATATGCCGGTAAACAGTATTTGGGAATGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTg                                                       <  2:38295/148‑1 (MQ=255)
                                                                                                   ggcggcATTGGATATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGCTGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCg                                                    <  2:148519/150‑1 (MQ=255)
                                                                                                               tATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTtgtg                                        <  1:1070799/150‑1 (MQ=255)
                                                                                                                   gCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAg                                    >  1:1773168/1‑150 (MQ=255)
                                                                                                                   gCGAGCAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAg                                    <  1:540521/150‑1 (MQ=255)
                                                                                                                     gaggagAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTAGGGAAGGTTCCGGCAATATTATGTGCGTGGATGTGTTGCGCGTTCTCAATCAGCAAGCGGGCGTATACGCCTTATTGTCGGAAGCATTTGTGGAAGTg                                  >  2:1131425/1‑150 (MQ=255)
                                                                                                                      aggagAGTGAGCTGCCGCGGCTTTACCGGGGGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGa                                 >  1:737037/1‑150 (MQ=255)
                                                                                                                         agagTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAg                              <  2:573481/150‑1 (MQ=255)
                                                                                                                           agTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAggg                            >  1:546762/1‑150 (MQ=255)
                                                                                                                            gTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                             <  2:405068/100‑1 (MQ=255)
                                                                                                                               aGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAgg                        >  1:904750/1‑150 (MQ=255)
                                                                                                                                  tGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATc                     <  2:1380028/150‑1 (MQ=255)
                                                                                                                                   gCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                    >  2:1363515/1‑150 (MQ=255)
                                                                                                                                   gCCGCGGCGTTACCGGGAGCTGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                    <  2:64514/150‑1 (MQ=255)
                                                                                                                                      gcgGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTa                 >  1:1464969/1‑150 (MQ=255)
                                                                                                                                        ggCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt               <  2:398977/150‑1 (MQ=255)
                                                                                                                                        ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGATGTGAAAGGGCAGGATCGCTAtt               >  1:1091218/4‑150 (MQ=255)
                                                                                                                                        ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt               >  2:1537082/4‑150 (MQ=255)
                                                                                                                                        ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt               >  2:1340802/4‑150 (MQ=255)
                                                                                                                                         tgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAttt              >  2:1048422/3‑150 (MQ=255)
                                                                                                                                         gcgttACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAttt              <  2:60844/147‑1 (MQ=255)
                                                                                                                                         gCTTTAGCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAttt              <  1:440234/150‑1 (MQ=255)
                                                                                                                                          gtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATTTTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTTTTGTCGGAAGCATTTGTGGAAGTGAAAGGGCTGGATCGCTAtttt             >  1:371962/2‑150 (MQ=255)
                                                                                                                                          gtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt             >  1:87322/2‑150 (MQ=255)
                                                                                                                                          gtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGGGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt             >  1:1323252/2‑150 (MQ=255)
                                                                                                                                          cTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGAAATATTATGTGCCTGGATGTGATGCGCGTGCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt             <  1:86968/150‑1 (MQ=255)
                                                                                                                                           tttACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTg            <  1:327519/150‑1 (MQ=255)
                                                                                                                                            ttACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGACTGGATGTCTTGCGCGTTCTCAATAAGCAATCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGAGCGCTATTTTGa           >  2:121637/1‑150 (MQ=255)
                                                                                                                                              aCCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATc         <  2:408918/150‑1 (MQ=255)
                                                                                                                                              aCCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATc         <  1:1565215/150‑1 (MQ=255)
                                                                                                                                               ccGGGAGATGCCGGTAAACAGTATTTGGGCAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTTGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcg        <  2:499627/150‑1 (MQ=255)
                                                                                                                                               ccGGGAGATGCCGGTAAACAGTAATTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcg        >  2:899406/1‑150 (MQ=255)
                                                                                                                                                   gagaTGCCGGTAAACAGTATATGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcgcgcg    <  2:605968/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCGGT  >  NC_000913/4415376‑4415674

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