Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 605,022 C→G 5.9% D135H (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,0220CG 5.9% 230.2 / 5.5 68D135H (GAT→CAT) nfsBdihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base G (2/2);  total (37/31)
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.14e-01

TAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTG  >  NC_000913/604877‑605168
                                                                                                                                                 |                                                                                                                                                  
tAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcat                                                                                                                                                  >  2:369342/1‑148 (MQ=255)
  gCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatcat                                                                                                                                               >  2:391964/1‑149 (MQ=255)
  gCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatcat                                                                                                                                               <  6:168063/149‑1 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCGCGGATGGCGGCGTCAAACCCTTAGATCTTTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGc                                                                                                                                             <  2:93797/149‑1 (MQ=255)
      ttCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATg                                                                                                                                              >  8:371580/1‑146 (MQ=255)
        ctctTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGa                                                                                                                                          <  1:18138/148‑1 (MQ=255)
           tttCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGAt                                                                                                                                         <  6:163217/146‑1 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                         <  1:192779/144‑2 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                         >  2:192779/1‑143 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                    >  7:299146/1‑149 (MQ=255)
                 aCCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAACCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGt                                                                                                                                >  6:375784/1‑149 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTg                                                                                                                               >  2:107371/1‑148 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTg                                                                                                                               >  7:351860/1‑148 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                              <  5:531/149‑1 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                              <  8:351860/149‑1 (MQ=255)
                            cATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                       <  6:76573/147‑1 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAACCCTGTTTTGCCATCCACTCTGCATCATCATGCAGCTCTTTACGGTGCATATCAGCg                                                                                                                     >  2:19725/1‑148 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                    <  1:68094/149‑1 (MQ=255)
                              tCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaa                                                                                                                   <  8:150791/149‑1 (MQ=255)
                                 aGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaagaa                                                                                                                >  3:225176/1‑149 (MQ=255)
                                          gcttCAAAACCTTAGATGGGTACCGCGTCCAGAACCAGAGCCGACACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGAcc                                                                                                       <  7:370399/146‑1 (MQ=255)
                                           cgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCt                                                                                                      <  1:107371/149‑1 (MQ=255)
                                           cgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCt                                                                                                      <  8:299146/149‑1 (MQ=255)
                                             tCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCtt                                                                                                     >  8:335152/1‑148 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                  <  4:48924/149‑1 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATc                                                                                                 <  4:218202/149‑1 (MQ=255)
                                                 aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCg                                                                                                >  2:18596/1‑149 (MQ=255)
                                                  aCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGt                                                                                               <  4:200680/149‑1 (MQ=255)
                                                         aTGGGTACCGCGTCCAGACCCAGAGCCGCTACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGc                                                                                        >  8:164070/1‑149 (MQ=255)
                                                         aTGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGc                                                                                        >  2:191536/1‑149 (MQ=255)
                                                             gTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTc                                                                                    <  7:371580/149‑1 (MQ=255)
                                                                 cgcgTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCtt                                                                                >  2:344626/1‑149 (MQ=255)
                                                                              agagCCGCCACGCCGAGCATGACTTTACCGACGTTGAGATAAACCTGTGTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTCTCGGTCGCGGCTTTCGCTTCCGGCGGGGCaaa                                                                   <  6:15941/149‑1 (MQ=255)
                                                                                    gccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGAt                                                                                                                                         >  1:253730/1‑73 (MQ=255)
                                                                                    gccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGAt                                                                                                                                         <  2:253730/73‑1 (MQ=255)
                                                                                    gccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGcc                                                             <  6:229068/149‑1 (MQ=255)
                                                                                     ccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCGGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCa                                                            >  8:331070/1‑149 (MQ=255)
                                                                                        cGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCg                                                         <  1:344626/149‑1 (MQ=255)
                                                                                         gCCGAGCAGGAAGTTACCGAATTTGAGATAAACCTGTATTTCCATCCACTCTGCATCATCATGCAGATATTTACGGTGCATATCAGCGAAGAACTTGATACCTTTATCGTTCGCGGCTTTCGCTTCCGGAGTGGCAAAGCGGCCAGCgg                                                        >  8:233197/1‑149 (MQ=255)
                                                                                           cGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCCGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCa                                                      <  2:241550/149‑1 (MQ=255)
                                                                                           cGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCa                                                      >  5:294342/1‑149 (MQ=255)
                                                                                                aGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                                 <  7:335152/149‑1 (MQ=255)
                                                                                                      ttACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCTTCTTCCTGGt                                            >  6:350704/1‑148 (MQ=255)
                                                                                                        aCCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTcaa                                         >  4:211662/1‑149 (MQ=255)
                                                                                                            aCGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                     >  6:321415/1‑149 (MQ=255)
                                                                                                            aCGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                     >  2:109756/1‑149 (MQ=255)
                                                                                                                 gagaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                                 <  3:109647/148‑1 (MQ=255)
                                                                                                                    aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                               >  7:74581/1‑147 (MQ=255)
                                                                                                                    aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCa                             <  5:321415/149‑1 (MQ=255)
                                                                                                                    aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCa                             <  1:191536/149‑1 (MQ=255)
                                                                                                                    aTAAACATGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                               >  2:264166/1‑147 (MQ=255)
                                                                                                                        aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGc                           >  8:223513/1‑147 (MQ=255)
                                                                                                                          cTGTTTTGCAATCCACTCTGCATCATCATGCAGATCTATACGGTTCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCATGGTCAACCACCAGCTTCAGCCAGa                       >  6:343916/1‑149 (MQ=255)
                                                                                                                           tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                     >  3:144936/1‑86 (MQ=255)
                                                                                                                           tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                     <  4:144936/86‑1 (MQ=255)
                                                                                                                           tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGAc                      <  1:264166/149‑1 (MQ=255)
                                                                                                                            gTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACa                     >  1:113802/1‑149 (MQ=255)
                                                                                                                             ttttGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                    >  8:245224/1‑149 (MQ=255)
                                                                                                                                 gCCATCCACTCTGCATCATCATTCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCTCTTCCTGCGGGGCAAAGCGGCCATCGGCATCTTCCTTGTCAACAACCATCTTCAGCCAGACATCGt                 <  7:233197/148‑1 (MQ=255)
                                                                                                                                    aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                               >  3:199292/1‑131 (MQ=255)
                                                                                                                                    aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                               <  4:199292/131‑1 (MQ=255)
                                                                                                                                    aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt             >  3:64652/1‑149 (MQ=255)
                                                                                                                                     tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATc            >  3:266741/1‑149 (MQ=255)
                                                                                                                                     tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATc            >  4:358496/1‑149 (MQ=255)
                                                                                                                                        aCTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgcg         >  7:92972/1‑149 (MQ=255)
                                                                                                                                         atctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                         <  1:187735/83‑1 (MQ=255)
                                                                                                                                         atctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                         >  2:187735/2‑84 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATc                   <  7:275149/137‑1 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATc                   >  8:275149/1‑137 (MQ=255)
                                                                                                                                            tGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttt     >  3:211369/1‑149 (MQ=255)
                                                                                                                                             gCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt    >  5:369667/1‑149 (MQ=255)
                                                                                                                                              catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGTGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttttt   <  7:279943/149‑1 (MQ=255)
                                                                                                                                              catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt       <  4:225176/145‑1 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCTTCCATCGCGGTTTTTg  >  5:341411/1‑149 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGt                 >  7:384404/1‑134 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGt                 <  8:384404/134‑1 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTg  >  6:73986/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTG  >  NC_000913/604877‑605168

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