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

TCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAT  >  NC_000913/1815224‑1815680
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
tCGACACTAACCCGGCAAATTACGAACCGACCTCGATTACCGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                <  2:61394/235‑1 (MQ=255)
   acacTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGc                                                                                                                                                                                                                            >  1:66346/1‑236 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACATCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGGCAGAAGCCC‑TTTTAGCATCGCCCTATTGTCGAT‑GGCTTTCGTGTTTAGATAAGCaacg                                                                                                                                                                                                                        >  2:696571/1‑234 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGGTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGGTTCAGTGTTGCGGTAAGCAAAg                                                                                                                                                                                                                        >  1:626724/1‑236 (MQ=255)
       tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCCTGGAAGGCAATAAAGGTCGCGAGCGCCGCCCATCGTTTGGCGAATATTATTCCCATCCGCGCCTGTTCTGGCTAAGTCAGACGCCA‑TTTTAGGAGCGCCCTCTTGTCGAA‑TGTTTCAGTGTTTAGTTAAggtaag                                                                                                                                                                                                                        >  2:267589/1‑231 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTTAATCATACCAGGAGCGCGGGGAAGGCAATAAAGTTCGCGGGGGCAGCCCATCCTTTGGCGAATATTACTCCCATCCGCGGCTGTTCTGGCTAAGTCAGGAGCCA‑TTTGAGCAACGCCCTATTTGCGCT‑GGGTTACGGGTTTCGTTAAGCAAAGtcgttcg                                                                                                                                                                                                                 >  2:350040/1‑236 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCC‑TTTGAGCAGCGCCAAATTGTCCAT‑GGTTTCAGGGGTGAGTTAAGCAAAGtcgttcg                                                                                                                                                                                                                 >  1:483170/1‑236 (MQ=255)
              gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCCCGTGGAAGGCAATAAAGGTCGCGAGCTCAGCCCCTCGTTTGGCGAAAATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCC‑TTTTAGCAGCGGCATATTTTCGATAGGTTTC‑GTGGTTAGTTAAGTAAAAtcatttg                                                                                                                                                                                                                 >  2:314086/1‑234 (MQ=255)
                   ttACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCCTACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGTGCGCAGCCCATCGTTTGGCGGATATTATTTCCATCCGCGTCTGTTCTGGCTAAGTCAGCCGCCATTTTGTTAAGCCCCTTATT‑TCGCT‑GGGTTCAGCGTTATGTTAAGCAAATTCGTTCGTCCgg                                                                                                                                                                                                            >  2:320312/1‑235 (MQ=255)
                     aCGAACCGAACTCGATTAACGAGAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGta                                                                                                                                                                                                           <  1:320312/235‑1 (MQ=255)
                                 cGATTAACGATAACGGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGa                                                                                                                                                                                                  <  1:267589/232‑1 (MQ=255)
                                   agtAACGATAACTGGCCGCGCGACACCCCGCCGGGGCCGAAACGCGGCTGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCgt                                                                                                                                                                                            <  2:383522/234‑1 (MQ=255)
                                            aaCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTTAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTTAGTTAAGCAAAGACGGTCGTCCGTATATTCGTGAGCGCGTTGTTGa                                                                                                                                                                                      >  2:184235/1‑233 (MQ=255)
                                              cTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAg                                                                                                                                                                                  >  1:740506/1‑235 (MQ=255)
                                                 gCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTgg                                                                                                                                                                              >  1:151138/1‑236 (MQ=255)
                                                  ccGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgc                                                                                                                                                                             <  2:66346/236‑1 (MQ=255)
                                                     ggcgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGgcgc                                                                                                                                                                           <  1:184235/234‑1 (MQ=255)
                                                       cgAAACACCGCCGGGGCCGAAACGCGGCAGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCata                                                                                                                                                                        <  2:504876/236‑1 (MQ=255)
                                                                  cgGGGCCGAAACGCGGCTGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCAct                                                                                                                                                             <  2:151138/236‑1 (MQ=255)
                                                                    gggCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCActct                                                                                                                                                           <  2:740506/236‑1 (MQ=255)
                                                                                            tCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                                                   <  2:113524/236‑1 (MQ=255)
                                                                                               tACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCgg                                                                                                                                >  2:480864/1‑236 (MQ=255)
                                                                                                          cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         >  1:587059/1‑232 (MQ=255)
                                                                                                          cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGAATCGaa                                                                                                                         >  1:567792/1‑232 (MQ=255)
                                                                                                                        aaGTTCGGGGTCGCAGCCCAGCGTATGGCGAATATTATTCCCATCCGCGTCCGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCAAT‑GGTTTCAGTGTTGAGTTAAGCAAAGGCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTg                                                                                                       <  2:447115/236‑1 (MQ=255)
                                                                                                                             cgcgAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGAATCGAACTGACTGACGACCAGCTGAata                                                                                                   >  1:150655/1‑235 (MQ=255)
                                                                                                                                  gcgcAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                                             <  2:57159/236‑1 (MQ=255)
                                                                                                                                                  ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGTACGTCaa                                                                             >  1:335260/1‑236 (MQ=255)
                                                                                                                                                  ggCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACg                                                                                 >  2:139681/1‑232 (MQ=255)
                                                                                                                                                   gCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCaa                                                                             <  1:139681/235‑1 (MQ=255)
                                                                                                                                                             ttCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaaa                                                                  <  1:350040/236‑1 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCg                                                                                     >  2:121547/1‑212 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCg                                                                                     <  1:121547/212‑1 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                                             >  1:193322/1‑236 (MQ=255)
                                                                                                                                                                  aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                                             <  2:587059/236‑1 (MQ=255)
                                                                                                                                                                    ccGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                           <  1:500847/236‑1 (MQ=255)
                                                                                                                                                                     cgcgTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         <  1:653106/173‑1 (MQ=255)
                                                                                                                                                                     cgcgTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                                                         >  2:653106/1‑173 (MQ=255)
                                                                                                                                                                           tGTTCTGGCTAAGTCAGGCGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GNGTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTa                                                    <  2:144824/236‑1 (MQ=255)
                                                                                                                                                                                            aCGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                <  2:221963/47‑1 (MQ=255)
                                                                                                                                                                                            aCGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGtt                                                                                                                                                                                                                                >  1:221963/1‑47 (MQ=255)
                                                                                                                                                                                                 a‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                                                   >  1:363812/1‑199 (MQ=255)
                                                                                                                                                                                                 a‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                                                   <  2:363812/199‑1 (MQ=255)
                                                                                                                                                                                                           gcgcCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGa                     >  1:1814/1‑236 (MQ=255)
                                                                                                                                                                                                           gcgcCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGa                     >  2:280192/1‑236 (MQ=255)
                                                                                                                                                                                                            cgcCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                  >  2:297287/1‑222 (MQ=255)
                                                                                                                                                                                                            cgcCATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                  <  1:297287/222‑1 (MQ=255)
                                                                                                                                                                                                              ccATATTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAg                  >  1:350773/1‑236 (MQ=255)
                                                                                                                                                                                                                atatTGTCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGaa                    <  2:350773/232‑1 (MQ=255)
                                                                                                                                                                                                                      tCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGTAGGTTATGGTAAAGGTCGCGTg          >  2:663690/1‑236 (MQ=255)
                                                                                                                                                                                                                      tCGAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTg          >  1:349825/1‑236 (MQ=255)
                                                                                                                                                                                                                        gAT‑GGTTTCGGTTGTGAGTTAAGCAATGTCGTTCGTACGTATATTCGTGAGCGCGTTGTTGCCCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTGTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        <  2:150655/236‑1 (MQ=255)
                                                                                                                                                                                                                        gAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        <  1:280192/236‑1 (MQ=255)
                                                                                                                                                                                                                        gAT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt        >  1:472113/1‑236 (MQ=255)
                                                                                                                                                                                                                         aT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTa       <  1:480864/236‑1 (MQ=255)
                                                                                                                                                                                                                         aT‑GGTTTCAGTGTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTa       >  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)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
TCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCA‑TTTGAGCAGCGCCATATTGTCGAT‑GGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAT  >  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: ‑