Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,814,906 T→C 5.7% F347S (TTC→TCC)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,9060TC 5.7% 223.2 / 3.3 70F347S (TTC→TCC) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base T (39/27);  new base C (2/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATC  >  NC_000913/1814765‑1815054
                                                                                                                                             |                                                                                                                                                    
ccTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTc                                                                                                                                               >  4:84077/1‑149 (MQ=255)
ccTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTc                                                                                                                                               >  3:562903/1‑149 (MQ=255)
  tCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGATTTCAAGTTCGa                                                                                                                                             >  3:507751/1‑149 (MQ=255)
   cACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGAc                                                                                                                                            <  2:290331/149‑1 (MQ=255)
   cACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGAc                                                                                                                                            <  3:458481/149‑1 (MQ=255)
    aCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACt                                                                                                                                           >  3:121893/1‑149 (MQ=255)
     cTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACtt                                                                                                                                          >  1:463885/1‑149 (MQ=255)
      tCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGGAGATGAATTCAAGTTCGACTTc                                                                                                                                         >  1:237845/1‑149 (MQ=255)
       cGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCg                                                                                                                                        <  4:170427/149‑1 (MQ=255)
          ttGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtc                                                                                                                                     <  4:11692/149‑1 (MQ=255)
          ttGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtc                                                                                                                                     <  1:509370/149‑1 (MQ=255)
              gATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtcttct                                                                                                                                 >  2:130737/1‑149 (MQ=255)
                  aaGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGAt                                                                                                                             <  1:133256/149‑1 (MQ=255)
                      acaAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATccaa                                                                                                                         >  3:55065/1‑149 (MQ=255)
                        aaaaaCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATccaacc                                                                                                                       >  3:502350/1‑149 (MQ=255)
                        aaaaaCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATccaacc                                                                                                                       >  2:313884/1‑149 (MQ=255)
                                aCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTAt                                                                                                               >  2:369252/1‑149 (MQ=255)
                                    gACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCg                                                                                                           >  1:107344/1‑149 (MQ=255)
                                         gACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATCCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGgaaga                                                                                                      >  2:33578/1‑149 (MQ=255)
                                         gACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATCCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGgaaga                                                                                                      >  2:69637/1‑149 (MQ=255)
                                           cccGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAAc                                                                                                    <  1:130737/149‑1 (MQ=255)
                                            ccGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACt                                                                                                   >  1:179273/1‑149 (MQ=255)
                                            ccGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACt                                                                                                   >  1:183326/1‑149 (MQ=255)
                                            ccGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACt                                                                                                   >  1:418424/1‑149 (MQ=255)
                                              ggACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTgg                                                                                                 <  3:199733/149‑1 (MQ=255)
                                               gACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGt                                                                                                <  2:204561/149‑1 (MQ=255)
                                                  ttCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGcc                                                                                             >  4:168165/1‑149 (MQ=255)
                                                  ttCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGcc                                                                                             >  2:224940/1‑149 (MQ=255)
                                                     cACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGt                                                                                          <  2:463885/149‑1 (MQ=255)
                                                       ccgccgCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTc                                                                                        <  2:94773/149‑1 (MQ=255)
                                                       ccgccgCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTc                                                                                        >  3:385099/1‑149 (MQ=255)
                                                          ccgcGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGc                                                                                     <  2:290776/149‑1 (MQ=255)
                                                                    tgGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaa                                                                           >  4:463374/1‑149 (MQ=255)
                                                                     gggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaa                                                                          >  3:328390/1‑149 (MQ=255)
                                                                     gggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaa                                                                          >  2:385520/1‑149 (MQ=255)
                                                                     gggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaa                                                                          >  3:581334/1‑149 (MQ=255)
                                                                      ggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTTGGCaaaa                                                                         >  3:426861/1‑149 (MQ=255)
                                                                      ggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaaa                                                                         <  3:188846/149‑1 (MQ=255)
                                                                       gAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAAt                                                                        <  2:335725/149‑1 (MQ=255)
                                                                        aaGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATg                                                                       <  2:179273/149‑1 (MQ=255)
                                                                         aGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATgg                                                                      >  3:26334/1‑149 (MQ=255)
                                                                          gCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGt                                                                     >  4:295572/1‑149 (MQ=255)
                                                                             aTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCt                                                                  >  1:595616/1‑149 (MQ=255)
                                                                             aTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCt                                                                  >  2:256864/1‑149 (MQ=255)
                                                                             aTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCt                                                                  >  1:289859/1‑149 (MQ=255)
                                                                                       gCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAAc                                                        <  4:219808/149‑1 (MQ=255)
                                                                                         gATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAAccc                                                      <  1:466540/149‑1 (MQ=255)
                                                                                             ttCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGAt                                                  <  1:442659/149‑1 (MQ=255)
                                                                                             ttCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGAt                                                  <  4:581334/149‑1 (MQ=255)
                                                                                             ttCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGAt                                                  <  2:183326/149‑1 (MQ=255)
                                                                                                 ggAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAAct                                              >  3:226493/1‑149 (MQ=255)
                                                                                                              ggCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGaaaa                                 <  3:295572/149‑1 (MQ=255)
                                                                                                                 ttCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGa                              >  2:399241/1‑149 (MQ=255)
                                                                                                                  tCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGaa                             <  4:138634/149‑1 (MQ=255)
                                                                                                                   ccAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAAc                            >  2:391814/1‑149 (MQ=255)
                                                                                                                    cAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACa                           >  1:295626/1‑149 (MQ=255)
                                                                                                                      gTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGCCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAgg                         >  3:97463/1‑149 (MQ=255)
                                                                                                                       ttGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATaa                                                >  3:419609/1‑125 (MQ=255)
                                                                                                                       ttGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATaa                                                <  4:419609/125‑1 (MQ=255)
                                                                                                                           ttCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCt                    >  1:409832/1‑149 (MQ=255)
                                                                                                                             ccTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCttt                  <  4:55065/149‑1 (MQ=255)
                                                                                                                               tGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTcc                >  1:569287/1‑149 (MQ=255)
                                                                                                                               tGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTcc                <  4:502350/149‑1 (MQ=255)
                                                                                                                               tGAAGAAGATGAATCCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTcc                <  1:33578/149‑1 (MQ=255)
                                                                                                                               tGAAGAAGATGAATCCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTcc                <  1:69637/149‑1 (MQ=255)
                                                                                                                                gaagaaGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCa               >  1:62518/1‑149 (MQ=255)
                                                                                                                                 aagaagATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCAt              >  2:217324/1‑149 (MQ=255)
                                                                                                                                     agaTGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTg          <  1:137604/149‑1 (MQ=255)
                                                                                                                                           aTTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCata    >  3:331279/1‑149 (MQ=255)
                                                                                                                                           aTTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCata    <  1:369252/149‑1 (MQ=255)
                                                                                                                                             tCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATc  <  4:121893/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATC  >  NC_000913/1814765‑1815054

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