Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,601 T→G 6.7% intergenic (+195/+45) yfeW → / ← yfeX penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,549,6010TG 6.7% 350.4 / 11.8 120intergenic (+195/+45)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (57/55);  new base G (2/6);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 126 total aligned reads displayed.

ACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAA  >  NC_000913/2549452‑2549750
                                                                                                                                                     |                                                                                                                                                     
 cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtc                                                                                                                                                      >  2:1587994/1‑150 (MQ=255)
  aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTTTACAACGCGCACTTGGTCAACAATCTGAGGCcagtca                                                                                                                                                     <  2:356180/150‑1 (MQ=255)
   aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtcag                                                                                                                                                    >  2:247927/1‑150 (MQ=255)
     cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTa                                                                                                                                                  >  2:296549/1‑150 (MQ=255)
     cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTa                                                                                                                                                  >  1:1057107/1‑150 (MQ=255)
     cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTa                                                                                                                                                  >  2:1125074/1‑150 (MQ=255)
      gCAAAACTGCCTGATGCGCTACGCTTTTCAGGCCTACGTGGAGCCCACTTTATATTGAAATAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGAATCCGGCGTGTACACCGCGCACTGGGTCAACAATCTGAGGCCAGTCAGTaa                                                                                                                                                 >  2:117076/1‑150 (MQ=25)
        aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg                                                                                                                                               >  1:52313/1‑150 (MQ=255)
            cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTAc                                                                                                                                           >  2:1035444/1‑150 (MQ=255)
              gCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTg                                                                                                                                         >  2:100752/1‑150 (MQ=255)
               ccTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGc                                                                                                                                        >  2:644171/1‑150 (MQ=255)
                  gATGCGCTACGCGTATCAGGCCTACGTTGATACAAAAATATATTGAATTAGCACTATTTCTTAGGCCTTATAAGGCGTTCACGCCGCTTCCGGCGTGTACAACGCTCACTTTGTATACAATATTAGGCCAGTCAGTAATGTTACTGCCTg                                                                                                                                     <  2:937460/150‑1 (MQ=255)
                   aTGCGCTACGCTTATCCGGCCTACGTGGATCCCACAATATTTTGTATTAGCACGTTTTCGTAGGCCTGATAAGGCGTTCACTCCGCATCCGGCGTGTACAACGCGCACTTGGTCTAAAATCTGAGGCCAGTCAGTAATGTTACTGCCTgg                                                                                                                                    <  2:1605076/150‑1 (MQ=38)
                   aTGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTgg                                                                                                                                    <  2:1057107/150‑1 (MQ=255)
                    tGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGc                                                                                                                                   <  2:1055394/150‑1 (MQ=255)
                     gcgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGcc                                                                                                                                  <  1:1654132/150‑1 (MQ=255)
                      cgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCt                                                                                                                                 <  1:429652/150‑1 (MQ=255)
                         tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTGGGTCAACTAACAGAGCCCAGTCAGTAATTTTACTGCCTGTACTTTa                                                                                                                              >  2:534438/1‑150 (MQ=25)
                         gtagcgcaTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTa                                                                                                                              <  1:216083/143‑1 (MQ=255)
                           cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTATGCCGTATAAGGCGTTCACGCCGCATCCGGCGTTTACAACGCTAACTTTGTAAACAATCTTAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                                                                                                            <  2:943050/150‑1 (MQ=37)
                           cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                                                                                                            <  1:814180/150‑1 (MQ=255)
                           cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                                                                                                            <  1:63632/150‑1 (MQ=255)
                                 tCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGa                                                                                                                      >  1:46061/1‑150 (MQ=255)
                                   aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAAt                                                                                                                    >  2:10323/1‑150 (MQ=255)
                                    ggCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATc                                                                                                                   <  1:1587994/150‑1 (MQ=255)
                                      ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAg                                                                                                                 <  2:293557/150‑1 (MQ=255)
                                      ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAg                                                                                                                 >  1:1574284/1‑150 (MQ=255)
                                       cTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGCTACTGCCTGGCCTTTAATCAATGAATCAGa                                                                                                                <  1:1120230/150‑1 (MQ=255)
                                         aCGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaaa                                                                                                              <  2:708738/150‑1 (MQ=255)
                                           gTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACg                                                                                                            >  1:722408/1‑150 (MQ=255)
                                              gAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCtt                                                                                                         <  1:759897/150‑1 (MQ=255)
                                              gAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCtt                                                                                                         <  1:947502/150‑1 (MQ=255)
                                                  ccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAg                                                                                                     >  2:788910/1‑150 (MQ=255)
                                                   cacaATATATTGAATTAGCACGATTTCGAAGGCCGGATAAGGCGTTCACGCCGCATCCGGCTTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGATACTGCCTGGCATTTAATCAAAGAATCAGAACCGATTaaatc                                                                                                    >  1:159397/1‑146 (MQ=25)
                                                     caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcgc                                                                                                  <  2:907391/150‑1 (MQ=255)
                                                       atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCa                                                                                                >  2:1395273/1‑150 (MQ=255)
                                                              gAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                         >  1:1177225/1‑150 (MQ=255)
                                                              gAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                         <  2:1624730/150‑1 (MQ=255)
                                                               aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACTCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTg                                                                                        <  2:1713173/150‑1 (MQ=255)
                                                                  tAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTcc                                                                                     <  1:644171/150‑1 (MQ=255)
                                                                   aGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCTCTTGGTCAACAATCTGAGGCCAGTCACTAATGTTACTGCCTGGCCTTTAATCAAAGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCa                                                                                    >  2:284999/1‑150 (MQ=255)
                                                                   aGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTATAATCAATGAATCAGAAACGCTTACCGCGCCATCAACTTGTCCa                                                                                    >  2:152868/1‑150 (MQ=255)
                                                                    gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCTCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAg                                                                                   <  1:296549/150‑1 (MQ=255)
                                                                     cACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGc                                                                                  <  2:1378160/150‑1 (MQ=255)
                                                                     cACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCAGCAACTTGTCCAGc                                                                                  <  2:722408/150‑1 (MQ=255)
                                                                       cGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGa                                                                                >  2:176941/1‑150 (MQ=255)
                                                                         aTTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACg                                                                              <  2:43546/150‑1 (MQ=255)
                                                                           ttCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGCAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGt                                                                            >  2:247574/1‑150 (MQ=255)
                                                                            tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTg                                                                           >  1:1645931/1‑150 (MQ=255)
                                                                             cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGGGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGc                                                                          <  2:831339/150‑1 (MQ=255)
                                                                             cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGc                                                                          >  2:1584959/1‑150 (MQ=255)
                                                                                    cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaataa                                                                   >  1:991353/1‑150 (MQ=255)
                                                                                     ggATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaataat                                                                  >  2:1495232/1‑150 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                <  1:135662/150‑1 (MQ=255)
                                                                                       aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                <  1:69874/150‑1 (MQ=255)
                                                                                         aaGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgcc                                                              <  1:653124/150‑1 (MQ=255)
                                                                                         aaGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgcc                                                              <  1:1035444/150‑1 (MQ=255)
                                                                                          aGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccg                                                             >  1:997162/1‑150 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccgc                                                            >  2:387989/1‑150 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccgc                                                            <  1:451223/150‑1 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccgc                                                            >  2:74942/1‑150 (MQ=255)
                                                                                             cGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                          <  2:684444/150‑1 (MQ=255)
                                                                                              gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                         <  2:885504/150‑1 (MQ=255)
                                                                                               tacACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGt                                                        >  2:358259/3‑150 (MQ=255)
                                                                                                 cACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                      >  1:141133/1‑150 (MQ=255)
                                                                                                   cgccggATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTAcc                                                    >  2:480301/1‑150 (MQ=255)
                                                                                                    gccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCg                                                   >  2:1159498/1‑150 (MQ=255)
                                                                                                     ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCgg                                                  >  1:462811/1‑150 (MQ=255)
                                                                                                        cATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGttt                                               <  2:1087731/150‑1 (MQ=255)
                                                                                                          tCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCTAAATAATAGCCGCCGGTTACCGGTTTgg                                             <  1:1392285/150‑1 (MQ=255)
                                                                                                          tCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTgg                                             <  1:1039557/150‑1 (MQ=255)
                                                                                                            cGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGACACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTg                                           >  2:695145/1‑150 (MQ=255)
                                                                                                                gtgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                       <  1:1125074/150‑1 (MQ=255)
                                                                                                                gtgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                       <  1:187990/150‑1 (MQ=255)
                                                                                                                gtgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                       <  2:515481/150‑1 (MQ=255)
                                                                                                                gtgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                       <  2:175150/150‑1 (MQ=255)
                                                                                                                 tgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACg                                      <  2:1651543/150‑1 (MQ=255)
                                                                                                                 tgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACg                                      >  2:1387773/1‑150 (MQ=255)
                                                                                                                   tACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCa                                    <  1:788910/150‑1 (MQ=255)
                                                                                                                        cgcgCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAACCGCTTACAGCGCCATCAACTTGTCCAGCGCCGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATc                               >  2:154888/1‑150 (MQ=255)
                                                                                                                          cgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGc                             <  1:1007803/150‑1 (MQ=255)
                                                                                                                          cgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGc                             <  2:1079227/150‑1 (MQ=255)
                                                                                                                          cgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGc                             >  1:702464/1‑150 (MQ=255)
                                                                                                                              cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCa                         <  1:358259/150‑1 (MQ=255)
                                                                                                                              cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCa                         <  2:141133/150‑1 (MQ=255)
                                                                                                                               ttGGTCAACATTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCAc                        <  1:285910/150‑1 (MQ=255)
                                                                                                                                  gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGTATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCt                     >  2:1511349/1‑150 (MQ=255)
                                                                                                                                    aaacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                   <  1:1694686/149‑1 (MQ=255)
                                                                                                                                     aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTAc                  >  1:1650419/1‑150 (MQ=255)
                                                                                                                                      acaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTAcc                 >  2:650549/1‑150 (MQ=255)
                                                                                                                                      acaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTAcc                 >  1:613031/1‑150 (MQ=255)
                                                                                                                                       caaTCTGATGCCAGTCAGTAATGTTACTGCCTTGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAATTAATAGCCGCCGGTTACCGGTTTGGTGAAACGCACCATCGCATCACGCTTCCCa                >  2:865032/1‑150 (MQ=255)
                                                                                                                                          tCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCg             <  2:52313/150‑1 (MQ=255)
                                                                                                                                           cTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCgg            <  2:991353/150‑1 (MQ=255)
                                                                                                                                              aGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACAGGTTTGGTGAAACGCCACATCGCATCACGCTTACCATCGGTAt         >  2:31203/1‑150 (MQ=255)
                                                                                                                                               ggCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCATCTTGTCCATCGACGGTGCGAAATTATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATc        <  1:152868/150‑1 (MQ=255)
                                                                                                                                                gCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCg       <  1:247927/150‑1 (MQ=255)
                                                                                                                                                gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCg       >  1:1421550/1‑150 (MQ=255)
                                                                                                                                                 ccAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGc      >  1:1197084/1‑150 (MQ=255)
                                                                                                                                                  cagtcagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGcc     >  2:277386/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAA  >  NC_000913/2549452‑2549750

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