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.7% 276.4 / NA 75E209D (GAA→GACyrbLMg(2+)‑starvation‑stimulated protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (37/37);  total (38/37)
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

TGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAA  >  NC_000913/3348936‑3349224
                                                                                                                                              |                                                                                                                                                  
tGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATc                                                                                                                                            >  1:785392/1‑151 (MQ=255)
 ggCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATcc                                                                                                                                           <  1:101696/151‑1 (MQ=255)
    gACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGc                                                                                                                                        <  1:457207/151‑1 (MQ=255)
             gAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGtgttgt                                                                                                                               <  1:879438/151‑1 (MQ=255)
                    tgttgtCTGCGTAAGCAAGGGCGGTTATGGAAACGATTTATCGCCCAAACGGCGCTGGCTATTGCCCTGCAGAAAGATATGCAACCCCGTGAGAGCAACACGCTCGCCCTCACTTCGCGCGACGCTTAATACAGCTGGTGTTGTTGCAt                                                                                                                          >  2:899675/1‑149 (MQ=255)
                          cTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGc                                                                                                                   >  1:925692/1‑150 (MQ=255)
                            gCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCt                                                                                                                  <  2:527260/149‑1 (MQ=255)
                            gCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTc                                                                                                                 <  1:104340/150‑1 (MQ=255)
                              gTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCATCCTCGTGAGAGCAAAACGCTGGCCCGCACTTCGCGCGCCGCTTAATCCAGCTGGTGTTGTTGCATAAACg                                                                                                                     <  2:505921/144‑1 (MQ=255)
                                aaGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCg                                                                                                             <  2:448754/150‑1 (MQ=255)
                                aaGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGc                                                                                                            >  1:811878/1‑151 (MQ=255)
                                aaGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGc                                                                                                            >  2:564634/1‑151 (MQ=255)
                                  gCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCt                                                                                                          >  2:556654/1‑151 (MQ=255)
                                    aaGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTa                                                                                                         >  1:833976/1‑150 (MQ=255)
                                    aaGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTa                                                                                                         >  2:165193/1‑150 (MQ=255)
                                       agCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTg                                                                                                     <  2:844520/151‑1 (MQ=255)
                                             tATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACAt                                                                                               >  2:489287/1‑151 (MQ=255)
                                                ggAAACGATTTATCTCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCt                                                                                             <  2:1170504/150‑1 (MQ=255)
                                                ggAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCt                                                                                             <  2:208432/150‑1 (MQ=255)
                                                ggAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTg                                                                                            >  1:839106/1‑151 (MQ=255)
                                                ggAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTg                                                                                            <  2:161661/151‑1 (MQ=255)
                                                ggAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTg                                                                                            >  1:411370/1‑151 (MQ=255)
                                                  aaaCGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGcc                                                                                          >  2:938287/1‑151 (MQ=255)
                                                   aaCGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCg                                                                                         >  2:629427/1‑151 (MQ=255)
                                                    aCGATTTATCGCCCAACCGTCGCTCGATATTGCCCTGCAGAAAGATCTGCACCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGCCGCGGAAGCCAGCGGGGGGGGGGGCAGAAACGGCGCAACGACGAACGGATACATCTGc                                                                                           <  2:172392/148‑1 (MQ=255)
                                                        tttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTaaa                                                                                     >  1:92808/1‑150 (MQ=255)
                                                         ttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAAt                                                                                   >  2:646914/1‑151 (MQ=255)
                                                         ttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAAt                                                                                   >  2:1034856/1‑151 (MQ=255)
                                                         ttATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAAt                                                                                   >  1:495216/1‑151 (MQ=255)
                                                            tCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATcc                                                                                 >  2:440841/1‑150 (MQ=255)
                                                             cGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCAc                                                                               >  1:804800/1‑151 (MQ=255)
                                                              gCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACg                                                                              >  1:61837/1‑151 (MQ=255)
                                                               cccAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGc                                                                             >  1:269717/1‑151 (MQ=255)
                                                                   aCCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTg                                                                          >  1:395202/1‑150 (MQ=255)
                                                                   aCCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGa                                                                         <  1:165200/151‑1 (MQ=255)
                                                                      ggCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCGCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACgg                                                                      <  2:807073/151‑1 (MQ=255)
                                                                      ggCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACgg                                                                      <  2:653704/151‑1 (MQ=255)
                                                                         gcTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTa                                                                   <  2:1041971/151‑1 (MQ=255)
                                                                          cTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTAc                                                                  >  1:768946/1‑151 (MQ=255)
                                                                            ggCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                <  1:1127145/151‑1 (MQ=255)
                                                                             gCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGc                                                                >  1:595498/1‑150 (MQ=255)
                                                                                aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACg                                                             <  2:777512/150‑1 (MQ=255)
                                                                                aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGGGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGt                                                            >  1:907153/1‑151 (MQ=255)
                                                                                 ttGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGt                                                            >  1:1024066/1‑150 (MQ=255)
                                                                                  tGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTa                                                           <  2:849355/150‑1 (MQ=255)
                                                                                          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATgg                                                  >  1:555587/1‑151 (MQ=255)
                                                                                           gAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGc                                                 <  2:73784/151‑1 (MQ=255)
                                                                                             aaGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGCGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCga                                               <  1:1135301/151‑1 (MQ=255)
                                                                                                  cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagg                                             >  2:956587/1‑148 (MQ=255)
                                                                                                    gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAACCGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                            <  2:614764/147‑1 (MQ=255)
                                                                                                    gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                            <  1:1087180/147‑1 (MQ=255)
                                                                                                    gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagg                                             <  1:165219/146‑1 (MQ=255)
                                                                                                    gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGAc                                         <  1:997216/150‑1 (MQ=255)
                                                                                                     cAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                       <  1:874600/151‑1 (MQ=255)
                                                                                                     cAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                       <  1:554206/151‑1 (MQ=255)
                                                                                                      aaCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                      >  1:102117/1‑151 (MQ=255)
                                                                                                      aaCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                      >  1:143818/1‑151 (MQ=255)
                                                                                                       aCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                      >  1:1080882/1‑150 (MQ=255)
                                                                                                            gTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAAc                                 <  2:405899/150‑1 (MQ=255)
                                                                                                            gTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACg                                <  1:790150/151‑1 (MQ=255)
                                                                                                               agagCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCCAGGAGACTTTGAAACg                                >  2:1036041/1‑148 (MQ=255)
                                                                                                                   cAAAACGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                         <  1:578165/151‑1 (MQ=255)
                                                                                                                      aaCGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTa                      <  1:543065/151‑1 (MQ=255)
                                                                                                                       aCGCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGAtt                       >  1:584764/1‑149 (MQ=255)
                                                                                                                         gCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAgg                    >  1:325705/1‑150 (MQ=255)
                                                                                                                         gCTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGt                   <  2:52759/151‑1 (MQ=255)
                                                                                                                          cTGGCCCTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGt                   <  2:848093/150‑1 (MQ=255)
                                                                                                                              cccTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc              >  2:468555/1‑151 (MQ=255)
                                                                                                                               ccTCGCTTTGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCGTAACGGATACATCTGCCGTAAAATCCACGCCTGACGGCTAAGCACGTAGCCCGATGTCGAGGAGACTTTGAAACGAAGTGGATTAGGAAATAc              >  2:525521/1‑150 (MQ=255)
                                                                                                                               ccTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc              >  1:808678/1‑150 (MQ=255)
                                                                                                                                cTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTTTTGCAAAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAGGGCGAGGAGACTTGGAAACGAAGTGGATTAGGTaa                 <  2:930586/146‑1 (MQ=255)
                                                                                                                                cTCACTTCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATAc              >  2:192566/1‑149 (MQ=255)
                                                                                                                                   aCTGCGCGCGACGCTTAATCCAGCTGGTGTTGTTGCATAAACTGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTg         <  2:719956/151‑1 (MQ=255)
                                                                                                                                         cgcgACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTa   <  1:277913/151‑1 (MQ=255)
                                                                                                                                          gcgACGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTaa  <  1:209761/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TGGCGACCTGGTCGAAGTGGTGTTGTCTGCGTAAGCAAGAGCGGTTATGGAAACGATTTATCGCCCAACCGGCGCTGGCTATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAA  >  NC_000913/3348936‑3349224

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