Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 605,025 C→G 5.3% D134H (GAT→CAT)  nfsB ← dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913605,0250CG 5.3% 329.5 / 8.4 95D134H (GAT→CAT) nfsBdihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base C (38/51);  new base G (3/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 108 total aligned reads displayed.

TAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGCACA  >  NC_000913/604877‑605172
                                                                                                                                                    |                                                                                                                                                   
tAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatc                                                                                                                                                     >  2:506985/1‑149 (MQ=255)
tAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatc                                                                                                                                                     >  3:138379/1‑149 (MQ=255)
   ccTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATg                                                                                                                                                  >  1:873939/1‑149 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatcatcc                                                                                                                                                 >  4:773696/1‑147 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGc                                                                                                                                                 >  4:682337/1‑149 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCCAGCAGGGCGGCGTTAAAACTTTGGATGGGTACCGAGTCCAGACCCAGAGCCGCCACGCCGAACAGCAAGTTACCGACGTTTAGATAAACCTGTTATGCCATCCACACGGcatcaacatac                                                                                                                                                 >  4:1188/1‑147 (MQ=255)
     tttCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCa                                                                                                                                                <  1:879536/149‑1 (MQ=255)
          ctTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGCACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACTGACGTTGAGATAACCCTGTTTTGCCATCCACTGTGCATCATCATCCAGATCt                                                                                                                                           <  4:747447/149‑1 (MQ=255)
          ctTTCAGACCAAATTCTGCATCGAGGATGGCAGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCt                                                                                                                                           <  4:338993/149‑1 (MQ=255)
           tttCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCtt                                                                                                                                          >  2:149386/1‑149 (MQ=255)
           tttCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCtt                                                                                                                                          <  3:166151/149‑1 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        <  2:81920/149‑1 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        >  3:144930/1‑149 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATAGGTACCGCGTCCAGACCTAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        <  4:525554/149‑1 (MQ=255)
              cAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       <  4:124947/149‑1 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  2:585534/1‑149 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  3:397600/1‑149 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  3:700229/1‑149 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  2:302214/1‑149 (MQ=255)
                      aTTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCata                                                                                                                               >  3:510320/1‑149 (MQ=255)
                         cTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCa                                                                                                                            <  4:620252/149‑1 (MQ=255)
                          tGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           >  4:354115/1‑149 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGagt                                                                                                                                            >  1:341972/1‑127 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGagt                                                                                                                                            <  2:341972/129‑3 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                        >  4:477823/1‑149 (MQ=255)
                              tCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaa                                                                                                                       <  1:331187/149‑1 (MQ=255)
                              tCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaa                                                                                                                       >  1:290074/1‑149 (MQ=255)
                                gAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaaga                                                                                                                     >  3:40439/1‑149 (MQ=255)
                                           cgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCt                                                                                                          <  3:1817/149‑1 (MQ=255)
                                             tCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       >  4:819331/1‑118 (MQ=255)
                                             tCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       <  3:819331/118‑1 (MQ=255)
                                              cAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       <  3:135566/149‑1 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  4:412774/1‑149 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  4:545235/1‑149 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  1:534210/1‑149 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  1:212451/1‑149 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGc                                                                                                                          <  4:344695/128‑1 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGc                                                                                                                          >  3:344695/1‑128 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCATCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATc                                                                                                     <  3:725057/149‑1 (MQ=255)
                                                 aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCg                                                                                                    >  4:115101/1‑149 (MQ=255)
                                                     ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTcg                                                                                                >  3:538066/1‑149 (MQ=255)
                                                           gggTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCtt                                                                                          <  2:85368/149‑1 (MQ=255)
                                                            ggTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                         <  4:40439/149‑1 (MQ=255)
                                                             gTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTc                                                                                        <  3:687681/149‑1 (MQ=255)
                                                                ccGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCt                                                                                     <  2:845205/149‑1 (MQ=255)
                                                                      ccAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGc                                                                               >  3:840008/1‑149 (MQ=255)
                                                                        aGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           <  1:81900/103‑1 (MQ=255)
                                                                        aGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           >  2:81900/1‑103 (MQ=255)
                                                                        aGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                             <  2:828238/149‑1 (MQ=255)
                                                                         gACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTg                                                                            <  3:545235/149‑1 (MQ=255)
                                                                          aCCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTgg                                                                           <  2:365293/149‑1 (MQ=255)
                                                                               gagCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAg                                                                      <  3:354115/149‑1 (MQ=255)
                                                                                 gccgccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCg                                                                    <  2:534210/149‑1 (MQ=255)
                                                                                 gccgccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCg                                                                    >  4:57561/1‑149 (MQ=255)
                                                                                  ccgccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCgg                                                                   <  1:808525/149‑1 (MQ=255)
                                                                                      cACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCAt                                                               >  4:451968/1‑149 (MQ=255)
                                                                                       aCGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGc                                                                               >  3:840228/1‑132 (MQ=255)
                                                                                         gCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCgg                                                            <  1:149386/149‑1 (MQ=255)
                                                                                          ccGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGc                                                           <  1:846172/149‑1 (MQ=255)
                                                                                                aGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                                     <  1:480089/149‑1 (MQ=255)
                                                                                                  gAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   <  4:679905/149‑1 (MQ=255)
                                                                                                  gAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   >  4:510648/1‑149 (MQ=255)
                                                                                                  gAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCGTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   <  1:585534/149‑1 (MQ=255)
                                                                                                    aGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTgg                                                 <  3:458018/149‑1 (MQ=255)
                                                                                                     gTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                <  4:397600/149‑1 (MQ=255)
                                                                                                      ttACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTc                                               <  3:451968/149‑1 (MQ=255)
                                                                                                       tACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTca                                              <  4:78755/149‑1 (MQ=255)
                                                                                                         ccGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTcaac                                            <  3:372392/149‑1 (MQ=255)
                                                                                                              gTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  1:847052/1‑149 (MQ=255)
                                                                                                               ttGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAg                                      >  2:89350/1‑149 (MQ=255)
                                                                                                                tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                                     <  2:503331/149‑1 (MQ=255)
                                                                                                                   gaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTc                                  >  4:363957/1‑149 (MQ=255)
                                                                                                                    aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCa                                 >  3:665742/1‑149 (MQ=255)
                                                                                                                       aaCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGcc                              >  3:844833/1‑149 (MQ=255)
                                                                                                                        aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCTTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCa                             >  2:594627/1‑149 (MQ=255)
                                                                                                                        aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCa                             <  2:212451/149‑1 (MQ=255)
                                                                                                                        aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCa                             <  1:168301/149‑1 (MQ=255)
                                                                                                                          cTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGa                           <  3:205124/149‑1 (MQ=255)
                                                                                                                          cTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGa                           >  2:196152/1‑149 (MQ=255)
                                                                                                                          cTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGa                           <  2:521320/149‑1 (MQ=255)
                                                                                                                                 gCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTc                    <  2:69864/149‑1 (MQ=255)
                                                                                                                                       cACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgc              >  1:170478/1‑149 (MQ=255)
                                                                                                                                        aCTCTGCATCATTATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTc                                  >  3:66230/1‑128 (MQ=255)
                                                                                                                                        aCTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTc                                  <  4:66230/128‑1 (MQ=255)
                                                                                                                                          tctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                                     >  4:409662/1‑123 (MQ=255)
                                                                                                                                          tctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                                     <  3:409662/123‑1 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAg                            <  4:227860/132‑1 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAg                            >  3:227860/1‑132 (MQ=255)
                                                                                                                                           ctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtt          >  2:129568/1‑149 (MQ=255)
                                                                                                                                            tGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATAGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttt         <  4:700229/149‑1 (MQ=255)
                                                                                                                                             gCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt        <  1:89350/149‑1 (MQ=255)
                                                                                                                                              catgatgatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttttt       >  1:597077/1‑149 (MQ=255)
                                                                                                                                              catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttttt       <  3:477823/149‑1 (MQ=255)
                                                                                                                                              catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttttt       <  4:76926/149‑1 (MQ=255)
                                                                                                                                              catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttttt       <  2:871793/149‑1 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATGAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTg      >  1:712511/1‑149 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTg      <  3:510648/149‑1 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTg      <  1:302214/149‑1 (MQ=255)
                                                                                                                                                tcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGc     >  3:253113/1‑149 (MQ=255)
                                                                                                                                                   tcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGcaca  >  1:533100/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGCACA  >  NC_000913/604877‑605172

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: