Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,815,451 T→G F529V (TTT→GTT)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,815,4510TG100.0% 285.0 / NA 74F529V (TTT→GTT) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (36/38);  total (36/38)

GACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAG  >  NC_000913/1815226‑1815676
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
gACACTACCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCaa                                                                                                                                                                                                                      <  1:677988/239‑1 (MQ=255)
          cGGCGAATTACGAACCGAACACGATTACCGATAACTGGCCGCCCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGtcgttcgt                                                                                                                                                                                                            <  2:834934/239‑1 (MQ=255)
            gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGGTTGGCGAAAATTATTCCCATCCGCGTCTGTTCTGGCCAAGTCAGACGCCATTTTAGCAGCGCCATATTGTCGATGGTTTCAGTGTTTAGGTAAGCAAAGACGTTCGTcc                                                                                                                                                                                                          >  2:928578/1‑239 (MQ=255)
                            acgcGATAAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGa                                                                                                                                                                                              <  1:280402/232‑1 (MQ=255)
                                         tAACTGGCCGCGCGAAGCACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGATCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAg                                                                                                                                                                              <  1:355799/238‑1 (MQ=255)
                                                ccGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgcgc                                                                                                                                                                       <  1:222393/238‑1 (MQ=255)
                                                  gcgcgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCata                                                                                                                                                                    <  1:719332/239‑1 (MQ=255)
                                                        aaCACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCCGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGAtc                                                                                                                                                              <  2:641132/239‑1 (MQ=255)
                                                             cgccgGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGa                                                                                                                                                                >  2:269284/1‑232 (MQ=255)
                                                                              ggcggTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCCTCGTTTTGCGAATATTATTCCCATCCCCGTCTGTTCCGGCTAAGTCAGACGCCATTTGAGCCGCGCCCTTTTGTCGATGGTTTTAGTGTTGAGTTAAACAACGTCGTTCGTCCGTATCTTCGTGACCGCGTTGTTGACCAGCTGGCGCACTTTGATCCTCCTCTGGCCCAGGCGggggt                                                                                                                                        >  2:261241/1‑238 (MQ=255)
                                                                                   ttttGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGaaa                                                                                                                                    <  1:928578/238‑1 (MQ=255)
                                                                                             tACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtc                                                                                                                                >  2:176886/1‑232 (MQ=255)
                                                                                             tACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtc                                                                                                                                >  1:709620/1‑232 (MQ=255)
                                                                                                    aGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCCGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                     <  2:362197/236‑1 (MQ=255)
                                                                                                     gcgcGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGttgttg                                                                                                                                                                                   >  1:900936/1‑173 (MQ=255)
                                                                                                     gcgcGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGttgttg                                                                                                                                                                                   <  2:900936/173‑1 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAgcg                                                                                                                                                                                           <  1:872090/162‑1 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAgcg                                                                                                                                                                                           >  2:872090/1‑162 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgc                                                                                                                                                                         >  2:608941/1‑180 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                         >  1:1058124/1‑196 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                         <  2:1058124/196‑1 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTAt                                                                                                                         <  2:790806/228‑1 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTAt                                                                                                                         >  1:790806/1‑228 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCAATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgc                                                                                                                                                                         <  1:608941/180‑1 (MQ=255)
                                                                                                        cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGGCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAactg                                                                                                                  >  1:709275/1‑235 (MQ=255)
                                                                                                        cgTGGAAGACAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAactgactg                                                                                                              >  2:235005/1‑239 (MQ=255)
                                                                                                         gTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAactgactg                                                                                                              <  2:709620/238‑1 (MQ=255)
                                                                                                          tGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                         >  1:883978/1‑194 (MQ=255)
                                                                                                          tGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                         <  2:883978/194‑1 (MQ=255)
                                                                                                              aGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTgacgac                                                                                                         <  1:1045851/238‑1 (MQ=255)
                                                                                                               ggCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCa                                                                                                       >  1:695120/1‑239 (MQ=255)
                                                                                                                 cAATAAAGTTCGCGAGCTCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAg                                                                                                      <  1:793673/238‑1 (MQ=255)
                                                                                                                  aaTAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCa                                                                                                       >  2:793673/1‑236 (MQ=255)
                                                                                                                          tCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATc                                                                                             >  2:335829/1‑238 (MQ=255)
                                                                                                                                      cccATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCg                                                                                 >  2:1020448/1‑238 (MQ=255)
                                                                                                                                                ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTcc                                                                                                                                                                                                          <  2:988710/107‑1 (MQ=255)
                                                                                                                                                ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTcc                                                                                                                                                                                                          >  1:988710/1‑107 (MQ=255)
                                                                                                                                                ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCatat                                                                                                                                                                   >  2:247000/1‑146 (MQ=255)
                                                                                                                                                ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCatat                                                                                                                                                                   <  1:247000/146‑1 (MQ=255)
                                                                                                                                                  cGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGtcgttcgt                                                                                                                                                                                                            <  1:384960/103‑1 (MQ=255)
                                                                                                                                                  cGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGtcgttcgt                                                                                                                                                                                                            >  2:384960/1‑103 (MQ=255)
                                                                                                                                                           ttCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAgg                                                           >  1:1025758/1‑239 (MQ=255)
                                                                                                                                                             cccATCCGCGTCTCTTCTGGCGAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGGTTCAGTGTTGAGTTAAGCAAAGGCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                                         <  2:27027/239‑1 (MQ=255)
                                                                                                                                                               cATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                       <  1:363178/239‑1 (MQ=255)
                                                                                                                                                                aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                                               >  2:244797/1‑166 (MQ=255)
                                                                                                                                                                aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                                               <  1:244797/166‑1 (MQ=255)
                                                                                                                                                                aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTAAACGGTCTGAAAAAGGatcc                                                       >  2:1051289/1‑238 (MQ=255)
                                                                                                                                                                 tCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCt                                                                                                                                                                            >  2:168842/1‑120 (MQ=255)
                                                                                                                                                                 tCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCt                                                                                                                                                                            <  1:168842/120‑1 (MQ=255)
                                                                                                                                                                     cgTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCata                                                                                                                                                                    <  2:343410/124‑1 (MQ=255)
                                                                                                                                                                     cgTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCata                                                                                                                                                                    >  1:343410/1‑124 (MQ=255)
                                                                                                                                                                              tGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGt                                         >  2:1013166/1‑238 (MQ=255)
                                                                                                                                                                                   aaGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCa                                   <  2:93729/239‑1 (MQ=255)
                                                                                                                                                                                     gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGt                                                                                                                                                                                                >  1:639961/1‑80 (MQ=255)
                                                                                                                                                                                     gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGt                                                                                                                                                                                                <  2:639961/80‑1 (MQ=255)
                                                                                                                                                                                     gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAgcg                                                                                                                                                                                           <  1:765450/85‑1 (MQ=255)
                                                                                                                                                                                     gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAgcg                                                                                                                                                                                           >  2:765450/1‑85 (MQ=255)
                                                                                                                                                                                     gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCAt                                  >  2:476679/1‑238 (MQ=255)
                                                                                                                                                                                          aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGGCCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCAtt                                 >  2:798644/1‑234 (MQ=255)
                                                                                                                                                                                               aTTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCt                                                                                                    >  2:756825/1‑162 (MQ=255)
                                                                                                                                                                                               aTTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCt                                                                                                    <  1:756825/162‑1 (MQ=255)
                                                                                                                                                                                                        gcgcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATATGAAAg              >  1:1021185/1‑239 (MQ=255)
                                                                                                                                                                                                          gcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGaa                <  1:335829/235‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt    <  1:261241/239‑1 (MQ=255)
                                                                                                                                                                                                                   tCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                       <  1:721938/187‑1 (MQ=255)
                                                                                                                                                                                                                   tCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                       >  2:721938/1‑187 (MQ=255)
                                                                                                                                                                                                                    cGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAg  <  1:738333/239‑1 (MQ=255)
                                                                                                                                                                                                                    cGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAg  <  1:1051289/239‑1 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGt                                                                     <  1:176038/166‑1 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGt                                                                     >  2:176038/1‑166 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCa                                   >  1:922018/1‑200 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCa                                   <  2:922018/200‑1 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTg                       <  1:825930/212‑1 (MQ=255)
                                                                                                                                                                                                                          tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTg                       >  2:825930/1‑212 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
GACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAG  >  NC_000913/1815226‑1815676

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