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,4510TG96.7% 197.5 / NA 61F529V (TTT→GTT) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (34/25);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAT  >  NC_000913/1815224‑1815680
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
tCGACACTAACCCGGCAAATTACGAACCGACCTCGATTACCGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                <  2:61394/235‑1 (MQ=255)
   acacTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGc                                                                                                                                                                                                                            >  1:66346/1‑236 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACATCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGGCAGAAGCCCTTTTAGCATCGCCCTATTGTCGATGGCTTTCGTGTTTAGATAAGCaacg                                                                                                                                                                                                                        >  2:696571/1‑234 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGGTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGGTTCAGTGTTGCGGTAAGCAAAg                                                                                                                                                                                                                        >  1:626724/1‑236 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCCTGGAAGGCAATAAAGGTCGCGAGCGCCGCCCATCGTTTGGCGAATATTATTCCCATCCGCGCCTGTTCTGGCTAAGTCAGACGCCATTTTAGGAGCGCCCTCTTGTCGAATGTTTCAGTGTTTAGTTAAggtaag                                                                                                                                                                                                                        >  2:267589/1‑231 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTTAATCATACCAGGAGCGCGGGGAAGGCAATAAAGTTCGCGGGGGCAGCCCATCCTTTGGCGAATATTACTCCCATCCGCGGCTGTTCTGGCTAAGTCAGGAGCCATTTGAGCAACGCCCTATTTGCGCTGGGTTACGGGTTTCGTTAAGCAAAGtcgttcg                                                                                                                                                                                                                 >  2:350040/1‑236 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCCTTTGAGCAGCGCCAAATTGTCCATGGTTTCAGGGGTGAGTTAAGCAAAGtcgttcg                                                                                                                                                                                                                 >  1:483170/1‑236 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCCCGTGGAAGGCAATAAAGGTCGCGAGCTCAGCCCCTCGTTTGGCGAAAATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCCTTTTAGCAGCGGCATATTTTCGATAGGTTTCGTGGTTAGTTAAGTAAAAtcatttg                                                                                                                                                                                                                 >  2:314086/1‑234 (MQ=255)
                   ttACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCCTACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGTGCGCAGCCCATCGTTTGGCGGATATTATTTCCATCCGCGTCTGTTCTGGCTAAGTCAGCCGCCATTTTGTTAAGCCCCTTATTTCGCTGGGTTCAGCGTTATGTTAAGCAAATTCGTTCGTCCgg                                                                                                                                                                                                            >  2:320312/1‑235 (MQ=255)
                     aCGAACCGAACTCGATTAACGAGAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGta                                                                                                                                                                                                           <  1:320312/235‑1 (MQ=255)
                                 cGATTAACGATAACGGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGa                                                                                                                                                                                                  <  1:267589/232‑1 (MQ=255)
                                   agtAACGATAACTGGCCGCGCGACACCCCGCCGGGGCCGAAACGCGGCTGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCgt                                                                                                                                                                                            <  2:383522/234‑1 (MQ=255)
                                            aaCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTTAGCAGCGCCATATTGTCGATGGTTTCAGTGTTTAGTTAAGCAAAGACGGTCGTCCGTATATTCGTGAGCGCGTTGTTGa                                                                                                                                                                                      >  2:184235/1‑233 (MQ=255)
                                              cTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAg                                                                                                                                                                                  >  1:740506/1‑235 (MQ=255)
                                                 gCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTgg                                                                                                                                                                              >  1:151138/1‑236 (MQ=255)
                                                  ccGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgc                                                                                                                                                                             <  2:66346/236‑1 (MQ=255)
                                                     ggcgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgcgc                                                                                                                                                                           <  1:184235/234‑1 (MQ=255)
                                                       cgAAACACCGCCGGGGCCGAAACGCGGCAGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCata                                                                                                                                                                        <  2:504876/236‑1 (MQ=255)
                                                                  cgGGGCCGAAACGCGGCTGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                             <  2:151138/236‑1 (MQ=255)
                                                                    gggCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCActct                                                                                                                                                           <  2:740506/236‑1 (MQ=255)
                                                                                            tCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                                                   <  2:113524/236‑1 (MQ=255)
                                                                                               tACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCgg                                                                                                                                >  2:480864/1‑236 (MQ=255)
                                                                                                          cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         >  1:587059/1‑232 (MQ=255)
                                                                                                          cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGAATCGaa                                                                                                                         >  1:567792/1‑232 (MQ=255)
                                                                                                                        aaGTTCGGGGTCGCAGCCCAGCGTATGGCGAATATTATTCCCATCCGCGTCCGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCAATGGTTTCAGTGTTGAGTTAAGCAAAGGCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTg                                                                                                       <  2:447115/236‑1 (MQ=255)
                                                                                                                             cgcgAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGAATCGAACTGACTGACGACCAGCTGAata                                                                                                   >  1:150655/1‑235 (MQ=255)
                                                                                                                                  gcgcAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                                             <  2:57159/236‑1 (MQ=255)
                                                                                                                                                  ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGTACGTCaa                                                                             >  1:335260/1‑236 (MQ=255)
                                                                                                                                                  ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACg                                                                                 >  2:139681/1‑232 (MQ=255)
                                                                                                                                                   gCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCaa                                                                             <  1:139681/235‑1 (MQ=255)
                                                                                                                                                             ttCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaaa                                                                  <  1:350040/236‑1 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCg                                                                                     >  2:121547/1‑212 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCg                                                                                     <  1:121547/212‑1 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                                             >  1:193322/1‑236 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                                             <  2:587059/236‑1 (MQ=255)
                                                                                                                                                                    ccGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                           <  1:500847/236‑1 (MQ=255)
                                                                                                                                                                     cgcgTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         <  1:653106/173‑1 (MQ=255)
                                                                                                                                                                     cgcgTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         >  2:653106/1‑173 (MQ=255)
                                                                                                                                                                           tGTTCTGGCTAAGTCAGGCGCCATTTGAGCAGCGCCATATTGTCGATGNGTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTa                                                    <  2:144824/236‑1 (MQ=255)
                                                                                                                                                                                            aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                <  2:221963/47‑1 (MQ=255)
                                                                                                                                                                                            aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                >  1:221963/1‑47 (MQ=255)
                                                                                                                                                                                                 aTTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                                                   >  1:363812/1‑199 (MQ=255)
                                                                                                                                                                                                 aTTTGAGCAGCGCCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                                                   <  2:363812/199‑1 (MQ=255)
                                                                                                                                                                                                           gcgcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGa                     >  1:1814/1‑236 (MQ=255)
                                                                                                                                                                                                           gcgcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGa                     >  2:280192/1‑236 (MQ=255)
                                                                                                                                                                                                            cgcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                  >  2:297287/1‑222 (MQ=255)
                                                                                                                                                                                                            cgcCATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                  <  1:297287/222‑1 (MQ=255)
                                                                                                                                                                                                              ccATATTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAg                  >  1:350773/1‑236 (MQ=255)
                                                                                                                                                                                                                atatTGTCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGaa                    <  2:350773/232‑1 (MQ=255)
                                                                                                                                                                                                                      tCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGTAGGTTATGGTAAAGGTCGCGTg          >  2:663690/1‑236 (MQ=255)
                                                                                                                                                                                                                      tCGATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTg          >  1:349825/1‑236 (MQ=255)
                                                                                                                                                                                                                        gATGGTTTCGGTTGTGAGTTAAGCAATGTCGTTCGTACGTATATTCGTGAGCGCGTTGTTGCCCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTGTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        <  2:150655/236‑1 (MQ=255)
                                                                                                                                                                                                                        gATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        <  1:280192/236‑1 (MQ=255)
                                                                                                                                                                                                                        gATGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        >  1:472113/1‑236 (MQ=255)
                                                                                                                                                                                                                         aTGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTa       <  1:480864/236‑1 (MQ=255)
                                                                                                                                                                                                                         aTGGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTa       >  2:89158/1‑236 (MQ=255)
                                                                                                                                                                                                                            ggTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGc     >  2:59252/1‑236 (MQ=255)
                                                                                                                                                                                                                              tttCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAACAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGCTGTGTAAGGTCGCCTGGTAGCGa   >  2:743481/1‑236 (MQ=255)
                                                                                                                                                                                                                                tCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaaa                                                                  >  1:736304/1‑171 (MQ=255)
                                                                                                                                                                                                                                tCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaaa                                                                  <  2:736304/171‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCCGGCGGGGGTGTGCAAGGTCCCGGTGTAGCgct  >  2:657344/1‑232 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
TCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAT  >  NC_000913/1815224‑1815680

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