Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,349,078 A→C E209D (GAA→GAC yrbL → Mg(2+)‑starvation‑stimulated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,349,0780AC98.8% 301.8 / ‑4.3 82E209D (GAA→GACyrbLMg(2+)‑starvation‑stimulated protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (49/32);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

CCACTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCT  >  NC_000913/3348932‑3349216
                                                                                                                                                  |                                                                                                                                          
ccACTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCtt                                                                                                                                        <  2:889700/151‑1 (MQ=255)
ccACTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCtt                                                                                                                                        >  2:1280574/1‑151 (MQ=255)
 cACTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTa                                                                                                                                       <  2:242938/151‑1 (MQ=255)
  aCTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTa                                                                                                                                       >  1:939568/1‑150 (MQ=255)
     ggCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATc                                                                                                                                    <  1:1269550/150‑1 (MQ=255)
            tGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTgg                                                                                                                             >  1:778089/1‑150 (MQ=255)
            tGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGgt                                                                                                                            >  1:730484/1‑151 (MQ=255)
             ggTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGgt                                                                                                                            <  2:1115676/150‑1 (MQ=255)
              gTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGgtgt                                                                                                                          >  2:158924/1‑151 (MQ=255)
                cGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGtgttg                                                                                                                        >  1:350362/1‑151 (MQ=255)
                      ggtgTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCAt                                                                                                                  >  2:998978/1‑151 (MQ=255)
                           tgtCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACg                                                                                                             >  2:65439/1‑151 (MQ=255)
                            gtCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACgg                                                                                                            >  1:597794/1‑151 (MQ=255)
                               tGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTc                                                                                                         >  2:103686/1‑151 (MQ=255)
                                gCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTc                                                                                                         >  1:519676/1‑150 (MQ=255)
                                acggaatcacgaacGGTCATAGAAACGATGTATAGCCCAACCGGCGCTGGCGATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTGGCGCGGCGATTAATCCAGCTGGTGTTGCTGCGTAGACGGcgc                                                                                                         >  1:1011917/14‑148 (MQ=255)
                                     aGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGcc                                                                                                   >  1:1093318/1‑151 (MQ=255)
                                      gCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCt                                                                                                  <  1:209214/151‑1 (MQ=255)
                                       cAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCt                                                                                                  >  2:1048809/1‑150 (MQ=255)
                                          gagCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACt                                                                                              <  2:9108/151‑1 (MQ=255)
                                              ggTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGAt                                                                                           >  2:1226201/1‑150 (MQ=255)
                                              ggTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGAt                                                                                           >  2:1286031/1‑150 (MQ=255)
                                              ggTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGAt                                                                                           >  1:1017959/1‑150 (MQ=255)
                                                ttATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATAc                                                                                         >  1:55553/1‑150 (MQ=255)
                                                ttATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATAc                                                                                         >  1:1204519/1‑150 (MQ=255)
                                                 tATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACa                                                                                        <  2:189927/150‑1 (MQ=255)
                                                 tATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACAt                                                                                       >  1:1027750/1‑151 (MQ=255)
                                                 tATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACAt                                                                                       >  1:285659/1‑151 (MQ=255)
                                                    ggAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTg                                                                                    <  1:205880/151‑1 (MQ=255)
                                                      aaaCGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGcc                                                                                  <  1:140330/151‑1 (MQ=255)
                                                        aCGATTTATCCCCCAACCGGCACTGGCTATTGCCCTGCAGTAAGATCTGCAACCTCGTGAGAGCAAAACGCTGCCCCTCACTTCGCGCGACGCTTAATCCAGCGGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGt                                                                                <  2:417096/151‑1 (MQ=255)
                                                         cGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTCCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTa                                                                               <  1:464036/151‑1 (MQ=255)
                                                           aTTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTaa                                                                              >  1:721830/1‑150 (MQ=255)
                                                            tttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTaaa                                                                             >  2:254945/1‑150 (MQ=255)
                                                             ttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTaaaa                                                                            <  1:479614/150‑1 (MQ=255)
                                                               aTCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAACATCTGCACCCTCGTGAGAGCAAAAGGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCAGCGCCTACCTGATACATCTGCCGTAAAATcc                                                                         <  2:330951/151‑1 (MQ=255)
                                                                 cGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCAc                                                                       >  2:287338/1‑151 (MQ=255)
                                                                  gCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACg                                                                      <  2:1049395/151‑1 (MQ=255)
                                                                    ccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGcc                                                                    <  2:1227582/151‑1 (MQ=255)
                                                                      aaCCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTg                                                                  <  1:1032362/151‑1 (MQ=255)
                                                                      aaCCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTg                                                                  >  1:55373/1‑151 (MQ=255)
                                                                          ggCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTTCATAAACGGCTCGCCGCCTAACTGATACATCTGCCGTAAAATCCACGcc                                                                    <  1:862022/145‑1 (MQ=255)
                                                                            cgcTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCt                                                            >  2:605732/1‑151 (MQ=255)
                                                                             gcTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCt                                                            >  1:561622/1‑150 (MQ=255)
                                                                               tGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCt                                                                   <  2:1282676/141‑1 (MQ=255)
                                                                               tGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCt                                                                   >  1:1282597/1‑141 (MQ=255)
                                                                                ggCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACg                                                         <  2:960260/150‑1 (MQ=255)
                                                                                  cTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                      >  1:141410/1‑151 (MQ=255)
                                                                                   tATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                     >  2:304516/1‑151 (MQ=255)
                                                                                   tATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                     >  1:1225720/1‑151 (MQ=255)
                                                                                    aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCAc                                                      <  1:1290211/149‑1 (MQ=255)
                                                                                    aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGt                                                    <  1:1245868/151‑1 (MQ=255)
                                                                                      tGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGt                                                    >  2:955697/1‑149 (MQ=255)
                                                                                       gCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGc                                                 >  1:1138121/1‑151 (MQ=255)
                                                                                          cTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCg                                              >  2:1097337/1‑151 (MQ=255)
                                                                                            gCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCg                                              <  2:165527/149‑1 (MQ=255)
                                                                                             cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATg                                           >  2:212488/1‑151 (MQ=255)
                                                                                             cAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATg                                           >  2:763589/1‑151 (MQ=255)
                                                                                                  aGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgag                                      <  1:114173/151‑1 (MQ=255)
                                                                                                  aGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCg                                        >  2:687955/1‑149 (MQ=255)
                                                                                                    aTCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagg                                     >  2:353622/1‑150 (MQ=255)
                                                                                                       tGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGAc                                 >  2:704020/1‑151 (MQ=255)
                                                                                                        gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgaggag                                   <  1:1140826/148‑1 (MQ=255)
                                                                                                        gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGgaga                                  <  1:1220149/149‑1 (MQ=255)
                                                                                                            ccTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTg                             >  1:868755/1‑150 (MQ=255)
                                                                                                               cGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAAc                         >  1:1044927/1‑151 (MQ=255)
                                                                                                                gTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACg                        <  1:430861/151‑1 (MQ=255)
                                                                                                                 tGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACg                        >  1:194317/1‑150 (MQ=255)
                                                                                                                  gagaGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGaa                      <  1:978908/151‑1 (MQ=255)
                                                                                                                   agagCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGa                       >  1:1300840/1‑149 (MQ=255)
                                                                                                                   agagCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAg                     <  2:823669/151‑1 (MQ=255)
                                                                                                                    cagCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGt                    >  1:1228563/2‑151 (MQ=255)
                                                                                                                     agCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTg                   >  2:850044/1‑151 (MQ=255)
                                                                                                                      gCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTgg                  >  2:588851/1‑151 (MQ=255)
                                                                                                                       cAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                 <  1:892811/151‑1 (MQ=255)
                                                                                                                        aaaaCGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAt                >  2:647088/1‑151 (MQ=255)
                                                                                                                        aaaaCGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAt                >  1:772308/1‑151 (MQ=255)
                                                                                                                         aaaCGCTGGACCTCACCTCGCGCGTCGCTTAATCCAGCTGTTGTTGTTGCATAAACGGCTCACCGCCTAACTGATATATCTGCCGTAAAATCCACGCCTGTCGTATACTCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAtt               <  1:1020679/151‑1 (MQ=255)
                                                                                                                          aaCGCTGGCCCTCTCTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTa              <  2:1050065/151‑1 (MQ=255)
                                                                                                                                   ccTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc      >  2:497028/1‑150 (MQ=255)
                                                                                                                                   ccTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACa     <  1:1145154/151‑1 (MQ=255)
                                                                                                                                      cACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCt  <  2:82669/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
CCACTGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCT  >  NC_000913/3348932‑3349216

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