Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 605,022 C→G 5.6% 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.6% 386.2 / 10.6 107D135H (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.47e-01

Only 100 of 112 total aligned reads displayed.

TGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTT  >  NC_000913/604875‑605166
                                                                                                                                                   |                                                                                                                                                
tgtAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatc                                                                                                                                                  >  2:55668/1‑148 (MQ=255)
    gCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGcatcatca                                                                                                                                              <  1:55668/148‑1 (MQ=255)
                   aCCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                       <  4:66594/140‑2 (MQ=255)
                   aCCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAGag                                                                                                                                       >  3:66594/1‑139 (MQ=255)
                     cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGa                                                                                                                                        >  2:383871/1‑137 (MQ=255)
                     cAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGa                                                                                                                                        <  1:383871/137‑1 (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)
                           cTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATc                                                                                                                       <  3:1113216/148‑1 (MQ=255)
                            tGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATc                                                                                                                                      >  4:428528/1‑132 (MQ=255)
                            tGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATc                                                                                                                                      <  3:428528/132‑1 (MQ=255)
                                 cGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATGATGATGCAg                                                                                                                                         <  4:227053/124‑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)
                                                  aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACgg                                                                                                                               <  4:1077650/117‑1 (MQ=255)
                                                  aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGAcc                                                                                                     <  1:274670/143‑1 (MQ=255)
                                                  aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGAcc                                                                                                     >  2:274670/1‑143 (MQ=255)
                                                  aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                >  4:393995/1‑148 (MQ=255)
                                                  aaaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTCTTGCCATCCACTCTGCATCATCATGCAGATCTTTACgg                                                                                                                               >  3:1077650/1‑117 (MQ=255)
                                                   aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                <  3:168922/115‑1 (MQ=255)
                                                   aaCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACg                                                                                                                                >  4:168922/1‑115 (MQ=255)
                                                    aCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                            >  1:295943/1‑134 (MQ=255)
                                                    aCCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                            <  2:295943/134‑1 (MQ=255)
                                                       ttCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTAt                                                                                                >  3:913761/1‑143 (MQ=255)
                                                        tCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTAc                                                                                                                                 >  4:130285/1‑109 (MQ=255)
                                                           aTGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATTCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCgg                                                                                       >  2:146533/1‑148 (MQ=255)
                                                                 aCCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCg                                                                                 <  3:393995/148‑1 (MQ=255)
                                                                           gACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATc                                                                                                                       >  3:352258/1‑100 (MQ=255)
                                                                               cAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCa                                                                                                                           >  3:50375/1‑92 (MQ=255)
                                                                               cAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCa                                                                                                                           <  4:50375/92‑1 (MQ=255)
                                                                                     cgccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCgg                                                             >  3:1053085/1‑148 (MQ=255)
                                                                                       ccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                            >  3:1068715/1‑99 (MQ=255)
                                                                                       ccACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACt                                                                                                            <  4:1068715/99‑1 (MQ=255)
                                                                                               aGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCAt                                                   >  3:343686/1‑148 (MQ=255)
                                                                                                   ggAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                               <  4:1053085/148‑1 (MQ=255)
                                                                                                     aaGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                             <  3:1089614/148‑1 (MQ=255)
                                                                                                        ttACCGACGTTGAGATAAACCTGTTTTGCCATTCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                          <  1:146533/148‑1 (MQ=255)
                                                                                                             gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTcg                                                                                          <  4:402672/95‑1 (MQ=255)
                                                                                                             gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTcg                                                                                          >  3:402672/1‑95 (MQ=255)
                                                                                                             gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCaa                                                                  >  1:709564/1‑119 (MQ=255)
                                                                                                             gACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCaa                                                                  <  2:709564/119‑1 (MQ=255)
                                                                                                                 ttGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                 <  4:343686/148‑1 (MQ=255)
                                                                                                                 ttGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCa                                 >  2:472072/1‑148 (MQ=255)
                                                                                                                  tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                  >  4:540890/1‑66 (MQ=255)
                                                                                                                  tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCga                                                                                                                  <  3:540890/66‑1 (MQ=255)
                                                                                                                  tGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAg                                >  4:40012/1‑148 (MQ=255)
                                                                                                                   gagaTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGc                               <  1:1087659/148‑1 (MQ=255)
                                                                                                                      aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                 <  4:62235/79‑1 (MQ=255)
                                                                                                                      aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                 >  3:62235/1‑79 (MQ=255)
                                                                                                                      aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTa                                                                                                 >  2:316283/1‑79 (MQ=255)
                                                                                                                      aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                   <  1:1072366/141‑1 (MQ=255)
                                                                                                                      aTAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCaacaac                                   >  2:1072366/1‑141 (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                                                                                   <  2:1021343/88‑1 (MQ=255)
                                                                                                                           ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                   >  3:386260/1‑88 (MQ=255)
                                                                                                                           ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCttt                                                                                   <  4:386260/88‑1 (MQ=255)
                                                                                                                           ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                             <  1:952468/142‑1 (MQ=255)
                                                                                                                           ccTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCtt                             >  2:952468/1‑142 (MQ=255)
                                                                                                                             tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCttt                                                                                                  >  3:292307/1‑71 (MQ=255)
                                                                                                                             tGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCttt                                                                                                  <  4:292307/71‑1 (MQ=255)
                                                                                                                                     cATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                          <  4:1045967/119‑1 (MQ=255)
                                                                                                                                     cATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGt                                          >  3:1045967/1‑119 (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)
                                                                                                                                       tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCt                                             >  3:667648/1‑114 (MQ=255)
                                                                                                                                       tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt           >  3:498068/1‑148 (MQ=255)
                                                                                                                                       tCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCAt           >  4:393948/1‑148 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                     <  3:25700/41‑1 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAg                                                                                                                     >  4:25700/1‑41 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCg                                                              <  2:15343/96‑1 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCg                                                              >  1:15343/1‑96 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCg                <  4:343754/142‑1 (MQ=255)
                                                                                                                                        ccACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCg                >  3:343754/1‑142 (MQ=255)
                                                                                                                                           ctctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTgg                                           <  1:528380/112‑1 (MQ=255)
                                                                                                                                            tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCa                                                          <  2:448033/96‑1 (MQ=255)
                                                                                                                                            tctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCa                                                          >  1:448033/1‑96 (MQ=255)
                                                                                                                                             ctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgcg       >  2:552604/1‑146 (MQ=255)
                                                                                                                                             ctGCATGATGATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATcgcg       <  1:552604/146‑1 (MQ=255)
                                                                                                                                             ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATc                                                        >  2:717191/1‑97 (MQ=255)
                                                                                                                                             ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                               <  1:838814/106‑1 (MQ=255)
                                                                                                                                             ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTc                                               >  2:838814/1‑106 (MQ=255)
                                                                                                                                             ctGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGt     >  1:852268/1‑148 (MQ=255)
                                                                                                                                                catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt  >  1:37210/1‑148 (MQ=255)
                                                                                                                                                catcatcatGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGtttt  >  4:805639/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTGTTTTGCCATCCACTCTGCATCATCATGCAGATCTTTACGGTGCATATCAGCGAAGAACTTGCGACCTTTATCGTTCGCGGCTTTCGCTTCCGGCGTGGCAAAGCGGCCATCGGCATCTTCCTGGTCAACAACCAGCTTCAGCCAGACATCGTCCATCGCGGTTTT  >  NC_000913/604875‑605166

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: