Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,614,830 G→A 5.4% Q388Q (CAG→CAA nikA → nickel/heme ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,614,8300GA 5.4% 376.4 / 5.7 110Q388Q (CAG→CAAnikAnickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (52/52);  new base A (3/3);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCG  >  NC_000913/3614684‑3614974
                                                                                                                                                  |                                                                                                                                                
ttGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGAt                                                                                                                                                <  3:1280056/149‑1 (MQ=255)
  gCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTg                                                                                                                                              <  2:1354290/149‑1 (MQ=255)
        aaaaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCat                                                                                                                                        <  3:536298/149‑2 (MQ=255)
          aaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGAt                                                                                                                                      <  1:1305409/149‑1 (MQ=255)
              cGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGtct                                                                                                                                  <  3:1213274/149‑1 (MQ=255)
                    gACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGa                                                                                                                            >  3:30645/1‑149 (MQ=255)
                         tgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       <  1:1157258/149‑1 (MQ=255)
                         tgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       >  1:627083/1‑149 (MQ=255)
                                 ggCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               <  4:30645/149‑1 (MQ=255)
                                  gCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCat                                                                                                                                        >  3:338580/1‑122 (MQ=255)
                                    aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 >  2:468011/1‑144 (MQ=255)
                                    aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 <  1:468011/144‑1 (MQ=255)
                                    aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagagag                                                                                                            <  2:627083/149‑1 (MQ=255)
                                    aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagagag                                                                                                            >  2:488694/1‑149 (MQ=255)
                                      aGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCa                                                                                                          >  4:438620/1‑149 (MQ=255)
                                             cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTgg                                                                                                                        >  4:485241/1‑128 (MQ=255)
                                             cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTgg                                                                                                                        <  3:485241/128‑1 (MQ=255)
                                             cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTa                                                                                                   >  2:1360717/1‑149 (MQ=255)
                                             cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGgcgc                                                                                                                                          >  2:744176/1‑110 (MQ=255)
                                                gagaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGc                                                                                                <  3:438620/149‑1 (MQ=255)
                                                  gaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGtctc                                                                                                                                 >  1:674788/1‑114 (MQ=255)
                                                  gaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGtctc                                                                                                                                 <  2:674788/114‑1 (MQ=255)
                                                     aaaTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTc                                                                                           >  3:801954/1‑149 (MQ=255)
                                                         ggTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTAt                                                                                                  >  1:412555/1‑138 (MQ=255)
                                                         ggTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTAt                                                                                                  <  2:412555/138‑1 (MQ=255)
                                                           tCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGa                                                                                                                                       <  3:1179831/99‑1 (MQ=255)
                                                           tCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGa                                                                                                                                       >  4:1179831/1‑99 (MQ=255)
                                                           tCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGt                                                                                                        >  2:247195/1‑130 (MQ=255)
                                                           tCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGt                                                                                                        <  1:247195/130‑1 (MQ=255)
                                                              gCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 <  1:1241549/118‑1 (MQ=255)
                                                              gCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 >  2:1241549/1‑118 (MQ=255)
                                                               ccGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCt                                                                                                    <  2:513055/130‑1 (MQ=255)
                                                                        aTTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGc                                                                                                                               >  1:35160/1‑94 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGa                                                                                                                                       >  3:894730/1‑85 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGc                                                                                                <  1:1338841/124‑1 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGc                                                                                                >  2:1338841/1‑124 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAgcg                                                                                       >  3:1036381/1‑133 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAgcg                                                                                       <  4:1036381/133‑1 (MQ=255)
                                                                         ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGc                                                                       <  4:801954/149‑1 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       <  2:1465273/95‑1 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       >  1:1465273/1‑95 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               <  3:883878/103‑1 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               >  4:883878/1‑103 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCa                                                                                                          >  4:157515/1‑108 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCa                                                                                                          <  3:157515/108‑1 (MQ=255)
                                                                               tttCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCat                                                                                                                                        <  3:594277/78‑2 (MQ=255)
                                                                                 tCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCa                                                                                                                                         <  2:1436130/75‑1 (MQ=255)
                                                                                 tCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCa                                                                                                                                         >  1:1436130/1‑75 (MQ=255)
                                                                                 tCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGAtttt                                                               >  4:42544/1‑149 (MQ=255)
                                                                                   gTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGa                                                                                    <  1:1262032/126‑1 (MQ=255)
                                                                                    ttCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCa                                                            <  1:488694/149‑1 (MQ=255)
                                                                                    ttCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCa                                                            >  3:1435028/1‑149 (MQ=255)
                                                                                     tCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGa                                                                                    >  2:1443635/1‑124 (MQ=255)
                                                                                     tCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                           >  3:1441812/1‑149 (MQ=255)
                                                                                     tCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                           <  1:1360717/149‑1 (MQ=255)
                                                                                      cATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               >  4:547637/1‑96 (MQ=255)
                                                                                      cATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTc                                                              <  3:780773/145‑1 (MQ=255)
                                                                                      cATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTc                                                              >  4:780773/1‑145 (MQ=255)
                                                                                      cATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCAAagaaga                                                                                                               <  3:547637/96‑1 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 <  1:1449055/89‑1 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaagaa                                                                                                                 >  2:1449055/1‑89 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTa                                                                                                   >  3:1008587/1‑103 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGttt                                                                           >  3:131664/1‑127 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGttt                                                                           <  4:131664/127‑1 (MQ=255)
                                                                                               cGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTggg                                                 <  3:42544/149‑1 (MQ=255)
                                                                                                     gTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       <  2:848924/73‑1 (MQ=255)
                                                                                                     gTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                       >  1:848924/1‑73 (MQ=255)
                                                                                                      ttAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGc                                                                                                >  4:753114/1‑95 (MQ=255)
                                                                                                      ttAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGc                                                                                                <  3:753114/95‑1 (MQ=255)
                                                                                                      ttAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGcc                                          >  3:350934/1‑149 (MQ=255)
                                                                                                          gCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCat                                        <  4:733340/147‑1 (MQ=255)
                                                                                                          gCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCat                                        >  3:733340/1‑147 (MQ=255)
                                                                                                            aaaTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGa                                    >  1:48109/1‑149 (MQ=255)
                                                                                                             aaTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGtttt                                                                          >  3:517153/1‑110 (MQ=255)
                                                                                                             aaTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGtttt                                                                          <  4:517153/110‑1 (MQ=255)
                                                                                                                      cGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCt                                                    <  3:1009298/123‑1 (MQ=255)
                                                                                                                      cGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCt                                                    >  4:1009298/1‑123 (MQ=255)
                                                                                                                       ggAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               >  3:326818/1‑63 (MQ=255)
                                                                                                                       ggAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAagaaga                                                                                                               <  4:326818/63‑1 (MQ=255)
                                                                                                                        gAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCa                                                                                          >  1:783667/1‑83 (MQ=255)
                                                                                                                        gAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCa                                                                                          <  2:783667/83‑1 (MQ=255)
                                                                                                                          aaTCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATc                                  <  4:369484/137‑1 (MQ=255)
                                                                                                                              aTTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCcaca                              >  3:543874/1‑137 (MQ=255)
                                                                                                                              aTTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCcaca                              <  4:543874/137‑1 (MQ=255)
                                                                                                                              aTTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGc                           <  4:676774/140‑1 (MQ=255)
                                                                                                                               ttCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGt                 >  3:739863/1‑149 (MQ=255)
                                                                                                                               ttCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGt                 <  2:48109/149‑1 (MQ=255)
                                                                                                                                 cAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTc               >  4:537516/1‑149 (MQ=255)
                                                                                                                                  aGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATc                                                                                                     >  1:1062112/1‑62 (MQ=255)
                                                                                                                                  aGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATc                                                                                                     <  2:1062112/62‑1 (MQ=255)
                                                                                                                                   ggCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCaa             <  4:1441812/149‑1 (MQ=255)
                                                                                                                                        atatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaa                                                                                                                    >  2:1294007/1‑41 (MQ=255)
                                                                                                                                        atatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaa                                                                                                                    <  1:1294007/41‑1 (MQ=255)
                                                                                                                                        atatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcgc        <  1:707410/149‑1 (MQ=255)
                                                                                                                                              gcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCg  >  4:119914/1‑149 (MQ=255)
                                                                                                                                              gcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCg  >  4:6917/1‑149 (MQ=255)
                                                                                                                                               ccAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATcc                                 <  2:704027/117‑1 (MQ=255)
                                                                                                                                               ccAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATcc                                 >  1:704027/1‑117 (MQ=255)
                                                                                                                                                  gATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCAt                                                                     >  3:833871/1‑78 (MQ=255)
                                                                                                                                                  gATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCAt                                                                     <  4:833871/78‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCG  >  NC_000913/3614684‑3614974

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