Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,814,762 A→C K299T (AAA→ACA)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,7620AC99.0% 323.6 / NA 96K299T (AAA→ACA) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (41/54);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATT  >  NC_000913/1814614‑1814906
                                                                                                                                                    |                                                                                                                                                
ctGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTaa                                                                                                                                                  <  1:624304/149‑1 (MQ=255)
       aCGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCAc                                                                                                                                          >  2:218805/1‑150 (MQ=255)
       aCGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACt                                                                                                                                         <  2:100717/151‑1 (MQ=255)
         gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCg                                                                                                                                       >  2:489293/1‑151 (MQ=255)
           gATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGtt                                                                                                                                     >  1:1059539/1‑151 (MQ=255)
           gATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGtt                                                                                                                                     >  2:340461/1‑151 (MQ=255)
           gATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGtt                                                                                                                                     >  2:170975/1‑151 (MQ=255)
           gATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGtt                                                                                                                                     <  1:120316/151‑1 (MQ=255)
           gATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGt                                                                                                                                      <  1:649516/150‑1 (MQ=255)
              gtgGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTgg                                                                                                                                  <  2:92334/151‑1 (MQ=255)
              gtgGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTgg                                                                                                                                  <  2:731305/151‑1 (MQ=255)
                   cGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACTCGTTTGGGATg                                                                                                                              >  1:1176907/1‑150 (MQ=255)
                    gATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACTCGTTTGGGATg                                                                                                                              >  1:86190/1‑149 (MQ=255)
                    gATGTCGGATCGCGGCAACCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAATCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGa                                                                                                                             <  1:1273568/150‑1 (MQ=255)
                     aTGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGGTTGGGATg                                                                                                                              >  2:392201/1‑148 (MQ=255)
                     aTGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACTCGTTTGGGATg                                                                                                                              >  1:965292/1‑148 (MQ=255)
                     aTGTCGGATCGCGGCACCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAg                                                                                                                           <  1:954379/151‑1 (MQ=255)
                       gTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGca                                                                                                                         <  2:1219805/151‑1 (MQ=255)
                       gTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGc                                                                                                                          >  2:71884/1‑150 (MQ=255)
                       gTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCCTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGc                                                                                                                          <  1:480334/150‑1 (MQ=255)
                        tCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGcac                                                                                                                        >  1:33074/1‑151 (MQ=255)
                         cGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGcac                                                                                                                        <  1:1198689/150‑1 (MQ=255)
                          ggATCGCGGCATCCTCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCCCACCCGTTTGGGATGAAGcaca                                                                                                                       >  2:676454/1‑150 (MQ=255)
                            aTCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTTCAGCCCCACACGTTTGGGATGAAACACaaa                                                                                                                     >  2:1034743/1‑150 (MQ=255)
                                cgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTc                                                                                                                <  2:924766/151‑1 (MQ=255)
                                    aTCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCAcc                                                                                                             >  2:984435/1‑150 (MQ=255)
                                     tCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGGTGAAGCACAAAAACTCACCgg                                                                                                           >  1:1130492/1‑151 (MQ=255)
                                       ccccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGAc                                                                                                         <  1:738907/151‑1 (MQ=255)
                                          cGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTg                                                                                                      <  2:237730/151‑1 (MQ=255)
                                           gCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTg                                                                                                      >  2:660767/1‑150 (MQ=255)
                                                tACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCgg                                                                                                >  2:313077/1‑151 (MQ=255)
                                                tACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCg                                                                                                 <  2:140828/150‑1 (MQ=255)
                                                    gCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACtt                                                                                            >  2:431006/1‑151 (MQ=255)
                                                    gCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACtt                                                                                            >  2:353570/1‑151 (MQ=255)
                                                     cACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACtt                                                                                            >  1:1280600/1‑150 (MQ=255)
                                                      aCCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTcc                                                                                          <  1:170973/151‑1 (MQ=255)
                                                      aCCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGGGACCCGGACTTcc                                                                                          >  1:188514/1‑151 (MQ=255)
                                                       ccATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCa                                                                                         <  1:816194/151‑1 (MQ=255)
                                                        cATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGTAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTTACTCGCTTGGGATGAAGCACCAAAACTCACCGGACGTTACCCGGCCTTCCAc                                                                                        >  2:955151/1‑151 (MQ=255)
                                                         aTGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAcc                                                                                       <  1:574832/151‑1 (MQ=255)
                                                          tGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAcc                                                                                       <  1:1196515/150‑1 (MQ=255)
                                                            gAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAccgcc                                                                                    <  1:290144/151‑1 (MQ=255)
                                                             aaGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAACACAAAAACACAACGGCCGCGACCCGGACTTCCAccgcc                                                                                    >  2:974295/1‑150 (MQ=255)
                                                              aGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGGTGGAGCACAAAAACTCACCGGGCGCGACCCGGACTTCCAccgccg                                                                                   >  1:400181/1‑150 (MQ=255)
                                                               ggCTTCGGTATTCACCCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCcgcg                                                                                 <  1:518901/151‑1 (MQ=255)
                                                                   tCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAg                                                                               <  2:439812/149‑1 (MQ=255)
                                                                    cGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGt                                                                              <  2:811268/149‑1 (MQ=255)
                                                                       tATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACCCGTTTGGGATGAAGCACAAAAACTCCCCGCACGTGACCCGGACTTCCACCGCCGCGAGTtgt                                                                           <  2:945438/149‑1 (MQ=255)
                                                                         ttCACTCCTTCCGCCTGATTAATGCCGAAGGGAGGGCAACGTTTGTACTTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTggg                                                                        <  1:50677/150‑1 (MQ=255)
                                                                           cacaCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTggg                                                                        >  1:735309/1‑148 (MQ=255)
                                                                              acCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGcc                                                                   >  1:711174/1‑150 (MQ=255)
                                                                                 ttCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGGCTCACTCGTTTGGGCTGAAGCACAAAAACTCACCGGACGTGACCCCGGCTTCCACCGCCGCGAGGTGTGGGAAGCCATTg                                                               >  2:1248088/1‑151 (MQ=255)
                                                                                   ccGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCAt                                                                 >  1:805200/1‑147 (MQ=255)
                                                                                   ccGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGaa                                                             <  2:320700/151‑1 (MQ=255)
                                                                                   ccGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGa                                                              <  1:182322/150‑1 (MQ=255)
                                                                                     gCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACTCGTTTGGGGCGAAGCACAAAAACTCACCGGACGTGAACCGGACTTCCACCGCCGCGAGATGAGGGAAGCCATTGaa                                                             >  1:110177/1‑149 (MQ=255)
                                                                                      ccTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGc                                                           >  1:7052/1‑150 (MQ=255)
                                                                                         gATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAgg                                                        <  1:821980/150‑1 (MQ=255)
                                                                                         gATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAgg                                                        <  1:1078654/150‑1 (MQ=255)
                                                                                           ttAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCg                                                      <  1:1246949/150‑1 (MQ=255)
                                                                                           ttAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCg                                                      >  2:15581/1‑150 (MQ=255)
                                                                                              aTGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTTACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGGttt                                                  >  2:879943/1‑151 (MQ=255)
                                                                                               tGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAtttt                                                 >  2:648886/1‑151 (MQ=255)
                                                                                                 ccGAAGGGAAGGCAACGTTTGTCCGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTcc                                               <  1:184151/151‑1 (MQ=255)
                                                                                                     agggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGa                                            >  2:151584/1‑150 (MQ=255)
                                                                                                      gggaaggCAACGTTTGTACGTTTCCACTGGAACCCACGGGCAGGTACAGCCTCACTCGTTGGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAAt                                          <  1:1290758/151‑1 (MQ=255)
                                                                                                       ggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCATTTGGGCTGAAGCACAAAAACCCACCGGACGGGACCCGGACGTCCACCGCCGGGCGGTGGGGGAAGACATTGAAGCAGGCGATTTTTCGGaa                                           >  2:212602/1‑149 (MQ=255)
                                                                                                           ggCAACGTTTGTACTTTCCCACTGGAAACCACTGGCTGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGa                                      <  1:676431/150‑1 (MQ=255)
                                                                                                           ggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGa                                      <  2:1076446/150‑1 (MQ=255)
                                                                                                            gCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGCGACCCGGACTTCCACCGCCGCGAGGTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAAc                                    >  1:336815/1‑151 (MQ=255)
                                                                                                             cAACGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACt                                   <  1:485246/151‑1 (MQ=255)
                                                                                                               aCGTTTGTACGTTTCCACTGGAAACCACTGGCAGGCACAGCCTCACTCGGTTGGGGTGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGGGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTg                                  >  2:665409/1‑150 (MQ=255)
                                                                                                                cGTTTGTACGTTTCCACTGGAACCCACGGGCAGGTACAGCCTCCCTCTTTTGGGATGAAGCACAAAACCTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                <  1:139771/151‑1 (MQ=255)
                                                                                                                cGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                <  2:320101/151‑1 (MQ=255)
                                                                                                                  tttGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCAAAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                <  1:1252074/149‑1 (MQ=255)
                                                                                                                     gTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTcc                           >  1:432625/1‑151 (MQ=255)
                                                                                                                       aCGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCa                          <  2:965350/150‑1 (MQ=255)
                                                                                                                        cGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCa                          <  1:1026452/149‑1 (MQ=255)
                                                                                                                         gTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACCCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGtt                       <  1:489280/151‑1 (MQ=255)
                                                                                                                           ttCCACTGGAAACCGCTGGCAGTTCCAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTg                      <  2:922450/150‑1 (MQ=255)
                                                                                                                             ccACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAt                    <  2:926027/150‑1 (MQ=255)
                                                                                                                             ccACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAtt                   <  2:33074/151‑1 (MQ=255)
                                                                                                                              cACTGGAAACCACTGGCAGGGACAGCCTCACTCGTTTGGGAGGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAtt                   >  2:408968/1‑150 (MQ=255)
                                                                                                                               aCTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTcc                 >  2:143545/1‑151 (MQ=255)
                                                                                                                                cTGGAAACCACTGGCAGGTACCGCCTCACTCGTTTGGGATGAAGCACAAAACCTCACCGGACGTGCCCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTcc                 <  2:295749/150‑1 (MQ=255)
                                                                                                                                cTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGAAGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCCTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTcc                 >  2:1214110/1‑150 (MQ=255)
                                                                                                                                 tGGAAGCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCt                <  2:1005036/150‑1 (MQ=255)
                                                                                                                                  ggAACCCACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCCCAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCt                <  2:1152856/149‑1 (MQ=255)
                                                                                                                                       cctcTGGCAGGTACAGCCTCACTCTTTTGGGATGAAGCACAAAAACTCACCGGACGTTACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTgaag            <  1:1187618/145‑1 (MQ=255)
                                                                                                                                       ccACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGaagaag         <  2:156576/151‑1 (MQ=255)
                                                                                                                                        cACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAagaaga        <  2:211338/151‑1 (MQ=255)
                                                                                                                                        cACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAagaaga        <  2:475549/151‑1 (MQ=255)
                                                                                                                                        cACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAagaaga        <  1:519203/151‑1 (MQ=255)
                                                                                                                                         aCTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGAt       <  2:286926/151‑1 (MQ=255)
                                                                                                                                           tGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATg      >  2:290666/1‑150 (MQ=255)
                                                                                                                                               aGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAAt   <  1:948317/149‑1 (MQ=255)
                                                                                                                                               aGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGCACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAAtt  <  1:208565/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATT  >  NC_000913/1814614‑1814906

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