Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 605,025 C→G 6.2% 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 6.2% 338.4 / 12.1 96D134H (GAT→CAT) nfsBdihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base C (47/43);  new base G (3/3);  total (50/46)
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.51e-01

Only 100 of 102 total aligned reads displayed.

TAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGCACA  >  NC_000913/604877‑605172
                                                                                                                                                    |                                                                                                                                                   
tAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatc                                                                                                                                                     >  4:773445/1‑149 (MQ=255)
 aGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatca                                                                                                                                                    >  3:527693/1‑149 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGc                                                                                                                                                 >  4:290855/1‑149 (MQ=255)
    cTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGc                                                                                                                                                 >  2:452182/1‑149 (MQ=255)
     tttCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCa                                                                                                                                                >  4:726714/1‑149 (MQ=255)
     tttCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCa                                                                                                                                                >  3:346243/1‑149 (MQ=255)
      ttCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGATGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAg                                                                                                                                               <  4:736944/149‑1 (MQ=255)
       tctcTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGa                                                                                                                                              <  4:663086/149‑1 (MQ=255)
        ctctTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGAt                                                                                                                                             <  4:401015/149‑1 (MQ=255)
           tttCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCtt                                                                                                                                          >  2:475114/1‑149 (MQ=255)
             tCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        <  1:144758/149‑1 (MQ=255)
             tCAGACCAAACTCTGCATCGAGGATGGCGGCGACAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGagtggat                                                                                                                                        <  2:322120/149‑7 (MQ=255)
              cAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       <  3:344499/149‑1 (MQ=255)
               aGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                      >  2:275711/1‑149 (MQ=255)
               aGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                      >  4:211663/1‑149 (MQ=255)
                  ccAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTg                                                                                                                                   >  4:66697/1‑149 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  <  1:389378/149‑1 (MQ=255)
                   cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  4:658828/1‑149 (MQ=255)
                        tCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATc                                                                                                                             <  4:783117/149‑1 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGTGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                        >  1:85753/1‑149 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGagtg                                                                                                                                           <  2:758875/130‑4 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGagtg                                                                                                                                           >  1:758875/1‑127 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCat                                                                                                                                <  1:538592/141‑1 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCat                                                                                                                                >  2:538592/1‑141 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                        >  2:495210/1‑149 (MQ=255)
                             aTCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                        >  3:526077/1‑149 (MQ=255)
                                 aGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCgaagaa                                                                                                                    <  1:587834/149‑1 (MQ=255)
                                    aTGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACtt                                                                                                                 <  3:290855/149‑1 (MQ=255)
                                     tGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTg                                                                                                                <  3:776761/149‑1 (MQ=255)
                                      ggcggcGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGc                                                                                                               <  3:73726/149‑1 (MQ=255)
                                      ggcggcGTCAAAACCTTCGATGGGGACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTGGAGGTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGc                                                                                                               >  4:385682/1‑149 (MQ=255)
                                           cgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCt                                                                                                          >  4:261503/1‑149 (MQ=255)
                                             tCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCttt                                                                                                        >  3:224612/1‑149 (MQ=255)
                                             tCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCtct                                                                                                        >  3:704617/1‑147 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                             >  3:323589/1‑109 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                             <  4:323589/110‑2 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  4:287978/1‑149 (MQ=255)
                                               aaaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      <  4:198866/149‑1 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATc                                                                                                     <  4:36519/149‑1 (MQ=255)
                                                aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCACCGAAGAACTTGCGACCCTTATc                                                                                                     <  2:476936/149‑1 (MQ=255)
                                                     ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       >  4:341525/1‑142 (MQ=255)
                                                     ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       <  3:341525/142‑1 (MQ=255)
                                                     ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTcg                                                                                                >  4:154363/1‑149 (MQ=255)
                                                          tGGGTACCGCGTCCAGACCCAGAGCCACCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCt                                                                                           <  2:29191/149‑1 (MQ=255)
                                                             gTACCGCGTCCAGACGCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTc                                                                                        >  4:308057/1‑149 (MQ=255)
                                                             gTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCt                                                                                                                                           >  2:25374/1‑98 (MQ=255)
                                                             gTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCt                                                                                                                                           <  1:25374/98‑1 (MQ=255)
                                                               aCCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGc                                                                                      <  1:795196/149‑1 (MQ=255)
                                                                  gcgTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTc                                                                                   <  1:452182/149‑1 (MQ=255)
                                                                        aGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGCCGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                             <  3:577805/149‑1 (MQ=255)
                                                                           cccAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGc                                                                          <  3:726714/149‑1 (MQ=255)
                                                                            ccAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCa                                                                         >  4:298307/1‑149 (MQ=255)
                                                                            ccAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCa                                                                         >  4:710134/1‑149 (MQ=255)
                                                                              agagCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCaaa                                                                       <  4:384734/149‑1 (MQ=255)
                                                                                agCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGc                                                                     <  3:211663/149‑1 (MQ=255)
                                                                                    gccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                             <  1:122468/137‑1 (MQ=255)
                                                                                    gccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGt                                                                             >  2:122468/1‑137 (MQ=255)
                                                                                       aCGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATc                                                              >  2:815575/1‑149 (MQ=255)
                                                                                        cGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCg                                                             <  3:864027/149‑1 (MQ=255)
                                                                                           cGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCa                                                          <  3:658828/149‑1 (MQ=255)
                                                                                            gAGCAGGAAGTTACCGACGTTGTGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCAt                                                         <  2:85753/149‑1 (MQ=255)
                                                                                            gAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCAt                                                         <  4:346243/149‑1 (MQ=255)
                                                                                             aGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATc                                                        <  4:208054/149‑1 (MQ=255)
                                                                                             aGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATc                                                        <  4:429833/149‑1 (MQ=255)
                                                                                                aGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                                     >  1:775341/1‑149 (MQ=255)
                                                                                                 ggAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTcc                                                    <  1:581542/149‑1 (MQ=255)
                                                                                                     gTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                <  1:320240/149‑1 (MQ=255)
                                                                                                        aCCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTcaa                                             >  1:695160/1‑149 (MQ=255)
                                                                                                            aCGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCCTCTTCCTGGTCaacaac                                         <  4:338908/149‑1 (MQ=255)
                                                                                                            aCGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                         <  1:678365/149‑1 (MQ=255)
                                                                                                              gTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  2:808285/1‑149 (MQ=255)
                                                                                                              gTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  1:90431/1‑149 (MQ=255)
                                                                                                              gTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  4:2381/1‑149 (MQ=255)
                                                                                                                   gaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGt                                                                                                                                    <  3:128129/51‑1 (MQ=255)
                                                                                                                   gaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGt                                                                                                                                    >  4:128129/1‑51 (MQ=255)
                                                                                                                    aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCa                                 <  4:406941/149‑1 (MQ=255)
                                                                                                                       aaCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGcc                              <  3:737772/149‑1 (MQ=255)
                                                                                                                        aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCa                             <  1:138467/149‑1 (MQ=255)
                                                                                                                         ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAg                            >  4:746115/1‑149 (MQ=255)
                                                                                                                         ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAg                            >  1:542570/1‑149 (MQ=255)
                                                                                                                            gTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACa                         >  1:343916/1‑149 (MQ=255)
                                                                                                                             ttttGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                        >  4:318361/1‑149 (MQ=255)
                                                                                                                             ttttGCCATCCACTCTGCATCAGCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                        <  3:122450/149‑1 (MQ=255)
                                                                                                                                 gCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTc                    <  2:695160/149‑1 (MQ=255)
                                                                                                                                  ccATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTcc                   >  4:760572/1‑149 (MQ=255)
                                                                                                                                    aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt                 >  3:98654/1‑149 (MQ=255)
                                                                                                                                    aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt                 >  2:256562/1‑149 (MQ=255)
                                                                                                                                     tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATc                >  1:533454/1‑149 (MQ=255)
                                                                                                                                       cACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgc              >  1:288941/1‑149 (MQ=255)
                                                                                                                                       cACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgc              >  2:214396/1‑149 (MQ=255)
                                                                                                                                         ctctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCgg            >  4:513971/1‑149 (MQ=255)
                                                                                                                                          tctGCATGATGATGCAGATCTTTACGGTCCATATCAGCGAAGAACTTGGGACCTTTATCGTTTGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCCTCGGCATCTCCCTCGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCggc           >  2:688708/1‑148 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt           >  2:673816/1‑149 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt           >  1:462819/1‑149 (MQ=255)
                                                                                                                                          tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt           <  3:513971/149‑1 (MQ=255)
                                                                                                                                           ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtt          <  4:526077/149‑1 (MQ=255)
                                                                                                                                               atcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTg      >  3:539206/1‑149 (MQ=255)
                                                                                                                                                 catcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                        >  4:614493/1‑129 (MQ=255)
                                                                                                                                                 catcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                        <  3:614493/129‑1 (MQ=255)
                                                                                                                                                   tcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGcaca  >  3:810688/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 < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: