Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 605,025 C→G 5.6% 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.6% 389.6 / 10.9 107D134H (GAT→CAT) nfsBdihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base C (52/49);  new base G (3/3);  total (55/52)
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.82e-01

Only 100 of 120 total aligned reads displayed.

GCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGCACA  >  NC_000913/604879‑605172
                                                                                                                                                  |                                                                                                                                                   
gCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatca                                                                                                                                                    <  1:55668/148‑1 (MQ=255)
               aCCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                             >  3:66594/1‑139 (MQ=255)
               aCCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                             <  4:66594/140‑2 (MQ=255)
                 cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGa                                                                                                                                              >  2:383871/1‑137 (MQ=255)
                 cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTg                                                                                                                                   >  4:408800/1‑148 (MQ=255)
                  aaaTTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCtt                                                                                                                                          >  3:52627/1‑140 (MQ=255)
                  aaaTTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCtt                                                                                                                                          <  4:52627/140‑1 (MQ=255)
                  aaaTTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGc                                                                                                                                  >  4:1113216/1‑148 (MQ=255)
                        tGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATc                                                                                                                                            <  3:428528/132‑1 (MQ=255)
                             cGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAg                                                                                                                                               >  3:227053/1‑124 (MQ=255)
                             cGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCa                                                                                                                                 >  2:724854/1‑138 (MQ=255)
                             cGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCa                                                                                                                                 <  1:724854/138‑1 (MQ=255)
                                    ggcggcGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTg                                                                                                                <  3:408800/148‑1 (MQ=255)
                                        gcgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        <  3:1097809/120‑1 (MQ=255)
                                        gcgTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTa                                                                                                                                        >  4:1097809/1‑120 (MQ=255)
                                              aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGAcc                                                                                                           >  2:274670/1‑143 (MQ=255)
                                              aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGAcc                                                                                                           <  1:274670/143‑1 (MQ=255)
                                              aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  4:393995/1‑148 (MQ=255)
                                              aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTCTTGCCATCCACTCTGCATCATCATGCAGATCTTTACgg                                                                                                                                     >  3:1077650/1‑117 (MQ=255)
                                               aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                      >  4:168922/1‑115 (MQ=255)
                                               aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                      <  3:168922/115‑1 (MQ=255)
                                                aCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                                  >  1:295943/1‑134 (MQ=255)
                                                   ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                      >  3:913761/1‑143 (MQ=255)
                                                    tCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       >  4:130285/1‑109 (MQ=255)
                                                    tCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                       <  3:130285/109‑1 (MQ=255)
                                                       aTGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATTCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCgg                                                                                             >  2:146533/1‑148 (MQ=255)
                                                             aCCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCg                                                                                       <  3:393995/148‑1 (MQ=255)
                                                                       gACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATc                                                                                                                             <  4:352258/100‑1 (MQ=255)
                                                                       gACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATc                                                                                                                             >  3:352258/1‑100 (MQ=255)
                                                                           cAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCa                                                                                                                                 <  4:50375/92‑1 (MQ=255)
                                                                                 cgccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCgg                                                                   >  3:1053085/1‑148 (MQ=255)
                                                                                   ccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                                  <  4:1068715/99‑1 (MQ=255)
                                                                                   ccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                                  >  3:1068715/1‑99 (MQ=255)
                                                                                      cGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATc                                                              >  3:587243/1‑148 (MQ=255)
                                                                                            gCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATc                                                        >  2:1087659/1‑148 (MQ=255)
                                                                                               ggAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                                     <  4:1053085/148‑1 (MQ=255)
                                                                                                 aaGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   <  3:1089614/148‑1 (MQ=255)
                                                                                                 aaGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   >  4:1089614/1‑148 (MQ=255)
                                                                                                    ttACCGACGTTGAGATAAACCTGTTTTGCCATTCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                <  1:146533/148‑1 (MQ=255)
                                                                                                         gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTcg                                                                                                <  4:402672/95‑1 (MQ=255)
                                                                                                         gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCaa                                                                        <  2:709564/119‑1 (MQ=255)
                                                                                                         gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCaa                                                                        >  1:709564/1‑119 (MQ=255)
                                                                                                             ttGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  2:472072/1‑148 (MQ=255)
                                                                                                             ttGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       <  4:343686/148‑1 (MQ=255)
                                                                                                              tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                        <  3:540890/66‑1 (MQ=255)
                                                                                                              tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAg                                      >  4:40012/1‑148 (MQ=255)
                                                                                                               gagaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                                     <  1:1087659/148‑1 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       >  3:62235/1‑79 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       <  1:316283/79‑1 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       <  4:62235/79‑1 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                       >  2:316283/1‑79 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                         >  2:1072366/1‑141 (MQ=255)
                                                                                                                  aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                         <  1:1072366/141‑1 (MQ=255)
                                                                                                                   tAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTcc                                                    >  2:842110/1‑129 (MQ=255)
                                                                                                                   tAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTcc                                                    <  1:842110/129‑1 (MQ=255)
                                                                                                                     aaCCTGTTTTTCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGc                               <  4:587243/148‑1 (MQ=255)
                                                                                                                      aCCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGcc                              >  3:481395/1‑148 (MQ=255)
                                                                                                                       ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                         >  1:1021343/1‑88 (MQ=255)
                                                                                                                       ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                         >  3:386260/1‑88 (MQ=255)
                                                                                                                       ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                         <  2:1021343/88‑1 (MQ=255)
                                                                                                                       ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                                   >  2:952468/1‑142 (MQ=255)
                                                                                                                       ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                                   <  1:952468/142‑1 (MQ=255)
                                                                                                                         tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCttt                                                                                                        >  3:292307/1‑71 (MQ=255)
                                                                                                                         tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCttt                                                                                                        <  4:292307/71‑1 (MQ=255)
                                                                                                                            tttGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACAt                        <  3:40012/148‑1 (MQ=255)
                                                                                                                                 cATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                <  4:1045967/119‑1 (MQ=255)
                                                                                                                                  aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                <  2:300640/118‑1 (MQ=255)
                                                                                                                                  aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                                >  1:300640/1‑118 (MQ=255)
                                                                                                                                  aTCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCa                  <  4:481395/148‑1 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           <  1:1037138/42‑1 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           >  2:1037138/1‑42 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCtt                                                      <  1:633489/111‑1 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCtt                                                      >  2:633489/1‑111 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                                   <  4:667648/114‑1 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt                 >  3:498068/1‑148 (MQ=255)
                                                                                                                                   tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt                 >  4:393948/1‑148 (MQ=255)
                                                                                                                                    ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           >  4:25700/1‑41 (MQ=255)
                                                                                                                                    ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                           <  3:25700/41‑1 (MQ=255)
                                                                                                                                    ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCg                                                                    <  2:15343/96‑1 (MQ=255)
                                                                                                                                    ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCg                      >  3:343754/1‑142 (MQ=255)
                                                                                                                                    ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCg                      <  4:343754/142‑1 (MQ=255)
                                                                                                                                       ctctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTgg                                                 <  1:528380/112‑1 (MQ=255)
                                                                                                                                        tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCa                                                                >  1:448033/1‑96 (MQ=255)
                                                                                                                                        tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCa                                                                <  2:448033/96‑1 (MQ=255)
                                                                                                                                         ctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgcg             <  1:552604/146‑1 (MQ=255)
                                                                                                                                         ctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgcg             >  2:552604/1‑146 (MQ=255)
                                                                                                                                         ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATc                                                              >  2:717191/1‑97 (MQ=255)
                                                                                                                                         ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATc                                                              <  1:717191/97‑1 (MQ=255)
                                                                                                                                         ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                                     >  2:838814/1‑106 (MQ=255)
                                                                                                                                         ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt           >  1:852268/1‑148 (MQ=255)
                                                                                                                                           gCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGttt         <  2:852268/148‑1 (MQ=255)
                                                                                                                                            catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt        >  1:37210/1‑148 (MQ=255)
                                                                                                                                            catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt        <  4:498068/148‑1 (MQ=255)
                                                                                                                                            catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt        >  4:805639/1‑148 (MQ=255)
                                                                                                                                                atcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTg                                                                            <  1:787528/76‑1 (MQ=255)
                                                                                                                                                atcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTg                                                                            >  2:787528/1‑76 (MQ=255)
                                                                                                                                                atcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       <  3:45720/113‑1 (MQ=255)
                                                                                                                                                atcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                       >  4:45720/1‑113 (MQ=255)
                                                                                                                                                atcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGca    >  3:763087/1‑148 (MQ=255)
                                                                                                                                                 tcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGc                               <  1:748712/120‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTTTGCACA  >  NC_000913/604879‑605172

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