Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,616 C→A 7.4% intergenic (+210/+30) yfeW → / ← yfeX penicillin binding protein 4B/porphyrinogen peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,549,6160CA 7.4% 312.9 / 10.4 95intergenic (+210/+30)yfeW/yfeXpenicillin binding protein 4B/porphyrinogen peroxidase
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base A (3/4);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

GATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTG  >  NC_000913/2549470‑2549765
                                                                                                                                                  |                                                                                                                                                     
gATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCTAGTCAGTAATGTTACTGCCTg                                                                                                                                                    <  1:898315/150‑1 (MQ=255)
gATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTg                                                                                                                                                    >  1:1217727/1‑150 (MQ=255)
  tGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                                                                                                  >  2:259267/1‑150 (MQ=255)
   gcgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGcc                                                                                                                                                 <  1:1581646/150‑1 (MQ=255)
   gcgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGcc                                                                                                                                                 <  2:639936/150‑1 (MQ=255)
     gcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                                                                                                                               <  1:1190648/150‑1 (MQ=255)
      cTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCttt                                                                                                                                              <  2:1217727/150‑1 (MQ=255)
       tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTa                                                                                                                                             <  2:973192/150‑1 (MQ=255)
        aCGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTaa                                                                                                                                            <  2:999275/150‑1 (MQ=255)
          gCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATc                                                                                                                                          <  2:1051578/150‑1 (MQ=255)
          gCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGCCtcagatt                                                                                                                                          <  2:1579387/150‑7 (MQ=25)
          gCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGCCtcagatt                                                                                                                                          <  2:1573812/150‑7 (MQ=25)
            ttATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                        <  2:310882/150‑1 (MQ=255)
             tATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                                                                                       <  2:847586/150‑1 (MQ=255)
             tATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                                                                                       <  2:84161/150‑1 (MQ=255)
              aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATg                                                                                                                                      <  1:423642/150‑1 (MQ=255)
                cAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGaa                                                                                                                                    <  1:347369/150‑1 (MQ=255)
                  ggCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATc                                                                                                                                  >  1:1212300/1‑150 (MQ=255)
                    ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAg                                                                                                                                >  1:538112/1‑150 (MQ=255)
                    ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAg                                                                                                                                >  2:1517210/1‑150 (MQ=255)
                           ggAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCt                                                                                                                         >  2:1008075/1‑150 (MQ=255)
                           ggAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCt                                                                                                                         <  2:1518297/150‑1 (MQ=255)
                             aGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTa                                                                                                                       >  1:643393/1‑150 (MQ=255)
                                ccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAg                                                                                                                    >  2:567323/1‑150 (MQ=255)
                                 cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgc                                                                                                                   >  1:605149/1‑150 (MQ=255)
                                 cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgc                                                                                                                   >  1:351423/1‑150 (MQ=255)
                                  acaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                                  >  2:269652/1‑150 (MQ=255)
                                  acaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGAACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                                  >  2:934930/1‑150 (MQ=255)
                                    aaTATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGcc                                                                                                                <  1:1099632/150‑1 (MQ=255)
                                        tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCa                                                                                                            >  1:556405/1‑150 (MQ=255)
                                              aTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGt                                                                                                      >  1:538737/1‑150 (MQ=255)
                                                 aGCACGATTTAGTAGGCCGGATAAGGCGTTCACGCATCATCCGGCGTGTACAAAGCGCAATAGGTCAACAATCTGAGGACAGTCAGTAATGTTACTGCCTTGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCa                                                                                                   >  2:1526462/1‑150 (MQ=38)
                                                   cACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTAAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGc                                                                                                 >  2:878121/1‑150 (MQ=255)
                                                     cGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGa                                                                                               >  2:1387818/1‑150 (MQ=255)
                                                        tttCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACgg                                                                                            >  1:1198704/1‑150 (MQ=255)
                                                         ttCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGt                                                                                           >  2:1622133/1‑150 (MQ=255)
                                                          tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTg                                                                                          <  2:308804/150‑1 (MQ=255)
                                                          tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTg                                                                                          >  2:334206/1‑150 (MQ=255)
                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGaaa                                                                                     <  2:975714/150‑1 (MQ=255)
                                                                gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCTCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaat                                                                                    <  2:1265537/150‑1 (MQ=255)
                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaataa                                                                                  <  2:1212300/150‑1 (MQ=255)
                                                                     aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                               <  2:181878/150‑1 (MQ=255)
                                                                     aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                               <  1:1538394/150‑1 (MQ=255)
                                                                     aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                               <  2:1484086/150‑1 (MQ=255)
                                                                       aaGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgcc                                                                             <  1:1210292/150‑1 (MQ=255)
                                                                         ggCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccgc                                                                           <  1:197234/150‑1 (MQ=255)
                                                                           tgtacACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                                         >  2:141871/5‑150 (MQ=255)
                                                                            gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                                        <  1:1519462/150‑1 (MQ=255)
                                                                              tCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGtt                                                                      >  2:1043533/1‑150 (MQ=255)
                                                                               cACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                                     >  2:855655/1‑150 (MQ=255)
                                                                                 cgccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTAcc                                                                   <  2:1201268/150‑1 (MQ=255)
                                                                                 cgccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTAcc                                                                   <  2:615094/150‑1 (MQ=255)
                                                                                   ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCgg                                                                 <  2:78007/150‑1 (MQ=255)
                                                                                   ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCgg                                                                 >  1:1188712/1‑150 (MQ=255)
                                                                                     gcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGtt                                                               <  1:319777/150‑1 (MQ=255)
                                                                                          cGGCGTTTACAACGCGCACTTGGGCAACAATCTGAGGCCAGGCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTg                                                          <  2:688463/150‑1 (MQ=255)
                                                                                          cGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTg                                                          >  1:1085969/1‑150 (MQ=255)
                                                                                              gtgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                                      <  2:950099/150‑1 (MQ=255)
                                                                                               tgtACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACg                                                     <  2:1505420/150‑1 (MQ=255)
                                                                                                 tACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCa                                                   <  1:259267/150‑1 (MQ=255)
                                                                                                 tACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCa                                                   <  2:1093543/150‑1 (MQ=255)
                                                                                                 tACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCa                                                   >  1:1458128/1‑150 (MQ=255)
                                                                                                 tACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCa                                                   <  1:1373871/150‑1 (MQ=255)
                                                                                                   cAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                                 <  2:556405/150‑1 (MQ=255)
                                                                                                   cAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCAACAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                                 >  2:1260842/1‑150 (MQ=255)
                                                                                                   cAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAACGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGATGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                                 >  2:7029/1‑150 (MQ=255)
                                                                                                     aCGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACAt                                               >  1:1099246/1‑150 (MQ=255)
                                                                                                       gcgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCg                                             <  1:1421596/150‑1 (MQ=255)
                                                                                                         gcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCa                                           <  2:1188712/150‑1 (MQ=255)
                                                                                                            cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCa                                        >  2:848590/1‑150 (MQ=255)
                                                                                                            cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCa                                        >  1:529224/1‑150 (MQ=255)
                                                                                                                  caacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                                  >  1:461580/1‑150 (MQ=255)
                                                                                                                  caacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                                  >  1:130863/1‑150 (MQ=255)
                                                                                                                  caacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                                  <  2:1040623/150‑1 (MQ=255)
                                                                                                                    acaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTAcc                                >  2:1151915/1‑150 (MQ=255)
                                                                                                                     caaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCa                               <  2:351423/150‑1 (MQ=255)
                                                                                                                         cTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTTCTAAATAATAGCCGCCGGTTACCGTTTTGGTGAAACGCAACATCGCATCACGCTTACCATCgg                           <  2:605149/150‑1 (MQ=255)
                                                                                                                         cTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCgg                           <  1:269652/150‑1 (MQ=255)
                                                                                                                          aaaGGCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGt                          >  2:1597989/3‑150 (MQ=255)
                                                                                                                           gAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTa                         >  2:1285996/1‑150 (MQ=255)
                                                                                                                              gCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGCTGAAACGCAACATCGCATCACGCTTACCATCGGTATCg                      >  2:530114/1‑150 (MQ=255)
                                                                                                                               ccAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                                      <  2:825595/117‑1 (MQ=255)
                                                                                                                               ccAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                                      >  1:825595/1‑117 (MQ=255)
                                                                                                                               ccAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGc                     >  1:364918/1‑150 (MQ=255)
                                                                                                                               ccAGGCAGTAACATTACTGACTGGCCTTTAATCAAAGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGc                     >  1:451379/1‑150 (MQ=255)
                                                                                                                                    cagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAAc                >  2:1135809/1‑150 (MQ=255)
                                                                                                                                     agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACa               >  2:521344/1‑150 (MQ=255)
                                                                                                                                          tGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTc          <  2:1494259/150‑1 (MQ=255)
                                                                                                                                           gTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTca         <  1:1387818/150‑1 (MQ=255)
                                                                                                                                             tACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagc       <  2:1576852/150‑1 (MQ=255)
                                                                                                                                              aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAAc                >  2:1125451/1‑140 (MQ=255)
                                                                                                                                              aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAAc                <  1:1125451/140‑1 (MQ=255)
                                                                                                                                              aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagca      >  1:916296/1‑150 (MQ=255)
                                                                                                                                               cTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagcag     <  1:1526462/150‑1 (MQ=255)
                                                                                                                                                 gCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGtt   >  1:1507019/1‑150 (MQ=255)
                                                                                                                                                  ccTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTtg  <  2:806979/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTG  >  NC_000913/2549470‑2549765

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