Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,616 C→A 8.0% intergenic (+210/+30) 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,6160CA 8.0% 391.9 / 18.6 112intergenic (+210/+30)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base C (56/47);  new base A (6/3);  total (62/50)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

Only 100 of 112 total aligned reads displayed.

TACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAG  >  NC_000913/2549477‑2549762
                                                                                                                                           |                                                                                                                                                  
tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                                                                                               <  2:1662753/145‑1 (MQ=255)
tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                                                                                               >  1:1662753/1‑145 (MQ=255)
  cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTaa                                                                                                                                         <  4:1086059/149‑1 (MQ=255)
   gCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                                                                                                                        >  4:108238/1‑149 (MQ=255)
      tATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                     <  3:1464003/149‑1 (MQ=255)
       aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                                                                                    <  3:108238/149‑1 (MQ=255)
              cTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAAt                                                                                                                                >  3:1614508/1‑146 (MQ=255)
              cTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAAt                                                                                                                                <  4:1614508/146‑1 (MQ=255)
                     gAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCt                                                                                                                      >  2:427333/1‑149 (MQ=255)
                        cccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                                                                                                                            <  1:1817793/124‑1 (MQ=255)
                        cccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                                                                                                                            >  2:1817793/1‑124 (MQ=255)
                              atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcgcg                                                                                                             >  2:1739573/1‑148 (MQ=255)
                              atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGcc                                                                                                             >  4:7306/1‑149 (MQ=255)
                                    tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGcc                                                                                                                                              >  2:335562/1‑110 (MQ=255)
                                    tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGcc                                                                                                                                              <  1:335562/110‑1 (MQ=255)
                                      aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                                     >  3:1492530/1‑149 (MQ=255)
                                           gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                     >  1:1741858/1‑112 (MQ=255)
                                           gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                     <  2:1741858/112‑1 (MQ=255)
                                               gATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGa                                                                                            <  1:427333/149‑1 (MQ=255)
                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCGATCAACTTGTCCAGCGACGGTGCGa                                                                                    <  1:1739573/149‑1 (MQ=255)
                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGAc                                                                                           <  1:199169/142‑1 (MQ=255)
                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGAc                                                                                           >  2:199169/1‑142 (MQ=255)
                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGa                                                                                    <  3:7306/149‑1 (MQ=255)
                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGt                                                                                        <  1:577393/144‑1 (MQ=255)
                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGt                                                                                        >  2:577393/1‑144 (MQ=255)
                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCa                                                                                                                                      >  2:539653/1‑97 (MQ=255)
                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCa                                                                                                                                      <  1:539653/97‑1 (MQ=255)
                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaat                                                                                 >  3:1546975/1‑149 (MQ=255)
                                                           cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaata                                                                                <  4:1546975/149‑1 (MQ=255)
                                                            ggATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATc                                                                                                                                       >  3:81336/1‑93 (MQ=255)
                                                            ggATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATc                                                                                                                                       <  4:81336/93‑1 (MQ=255)
                                                             gatgcggCGTGAACGCCTTATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaaa                                                                                                                          <  1:633987/100‑1 (MQ=255)
                                                              aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                     <  2:329713/93‑1 (MQ=255)
                                                              aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                     >  1:329713/1‑93 (MQ=255)
                                                                     gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCt                                                                                                                                             >  3:556892/1‑78 (MQ=255)
                                                                     gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCt                                                                                                                                             <  4:556892/78‑1 (MQ=255)
                                                                      tacACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCt                                                                                                                                                  <  4:1540104/70‑1 (MQ=38)
                                                                       tCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAAc                                                                                                       <  2:1397203/114‑1 (MQ=255)
                                                                       tCACGCCGCATACGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAAc                                                                                                       >  1:1397203/1‑114 (MQ=255)
                                                                       acACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                             >  2:1520762/2‑140 (MQ=255)
                                                                         aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                                  >  2:60555/1‑149 (MQ=255)
                                                                         aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                                  >  1:868569/1‑149 (MQ=255)
                                                                            ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCg                                                               >  2:1379256/1‑149 (MQ=255)
                                                                            ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCg                                                               <  4:1492530/149‑1 (MQ=255)
                                                                              gcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGt                                                             >  2:167325/1‑149 (MQ=255)
                                                                                 tCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTg                                                          <  2:1852173/149‑1 (MQ=255)
                                                                                            cAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCaa                                               <  1:1379256/149‑1 (MQ=255)
                                                                                             aaCGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                              <  2:868569/149‑1 (MQ=255)
                                                                                             aaCGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                              >  1:1823243/1‑149 (MQ=255)
                                                                                                gcgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATc                                           >  4:35701/1‑149 (MQ=255)
                                                                                                     cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGa                                                                                                                                  <  3:1703954/57‑1 (MQ=255)
                                                                                                     cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGa                                                                                                                                  >  4:1703954/1‑57 (MQ=255)
                                                                                                         gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                               <  2:625091/72‑1 (MQ=255)
                                                                                                         gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                               >  1:625091/1‑72 (MQ=255)
                                                                                                         gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGc                                  <  1:167325/149‑1 (MQ=255)
                                                                                                          tCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTGGTGAAACGCAACATCGCATCACGCtt                                >  1:590668/1‑149 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCATCTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                                     <  2:593924/111‑1 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                                     >  1:1149141/1‑79 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                             <  1:21133/103‑1 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                        <  1:1670166/124‑1 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                                   >  3:772742/1‑129 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                                   <  4:772742/129‑1 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                               >  3:786099/1‑149 (MQ=255)
                                                                                                            aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                               >  3:1663097/1‑149 (MQ=255)
                                                                                                                 tCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                             <  3:1771941/98‑1 (MQ=255)
                                                                                                                  cTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCg                         <  1:60555/149‑1 (MQ=255)
                                                                                                                     aGGCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTa                      >  1:1201250/1‑149 (MQ=255)
                                                                                                                       gCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATc                    <  4:786099/149‑1 (MQ=255)
                                                                                                                       gCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATc                    <  2:1823243/149‑1 (MQ=255)
                                                                                                                       gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                        <  3:1176962/113‑1 (MQ=39)
                                                                                                                       gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                        >  4:1176962/1‑113 (MQ=39)
                                                                                                                        ccAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTg                                                       <  1:1071989/113‑1 (MQ=255)
                                                                                                                        ccAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTg                                                       >  2:1071989/1‑113 (MQ=255)
                                                                                                                         caggcagtAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATc                    >  1:1208619/1‑147 (MQ=255)
                                                                                                                         caggcagtAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGc                  >  2:1335661/1‑149 (MQ=255)
                                                                                                                         caggcagtAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGc                  >  4:248999/1‑149 (MQ=255)
                                                                                                                           gtcagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGcc                 <  3:1183922/148‑1 (MQ=255)
                                                                                                                           gtcagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGcc                 >  4:1183922/1‑148 (MQ=255)
                                                                                                                           gtcagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCa                >  1:1448464/1‑149 (MQ=255)
                                                                                                                              agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                             <  1:1293439/85‑1 (MQ=255)
                                                                                                                              agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                             >  2:1293439/1‑85 (MQ=255)
                                                                                                                              agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                            >  4:712663/1‑86 (MQ=255)
                                                                                                                              agtAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAAc             <  2:1201250/149‑1 (MQ=255)
                                                                                                                               gtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                                      <  2:948472/91‑1 (MQ=255)
                                                                                                                               gtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                                      >  1:948472/1‑91 (MQ=255)
                                                                                                                                tcaTGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTTCAGCGACGCTGCGAAATAATACCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCACACAt           >  1:728563/3‑149 (MQ=255)
                                                                                                                                tAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACa                                             <  1:642115/115‑1 (MQ=255)
                                                                                                                                 aaTGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATg          <  3:35701/149‑1 (MQ=255)
                                                                                                                                  acattACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGc         <  3:286369/146‑1 (MQ=255)
                                                                                                                                  aTGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGc         >  1:1027773/1‑149 (MQ=255)
                                                                                                                                    gTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTc       >  2:1270988/1‑149 (MQ=255)
                                                                                                                                     ttACTGCGTGGGCTTTAATCAATAAGTCAGAAACGCTTACAGCGCCAACAACTTGTCCAGCGATCGTGCGAAAGAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATGGCATGACGCTTACCATCGGTATCGCCAAACATGCTca      >  4:1321330/1‑149 (MQ=255)
                                                                                                                                     ttACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTca      >  4:1577006/1‑149 (MQ=255)
                                                                                                                                      tACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcag     >  1:902833/1‑149 (MQ=255)
                                                                                                                                       aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                        <  4:1535000/97‑1 (MQ=255)
                                                                                                                                       aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                        >  3:1535000/1‑97 (MQ=255)
                                                                                                                                       aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagc    >  1:161889/1‑149 (MQ=255)
                                                                                                                                       aCTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagc    >  2:283247/1‑149 (MQ=255)
                                                                                                                                        cTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagca   <  4:1663097/149‑1 (MQ=255)
                                                                                                                                         tGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTcagcag  >  2:1682831/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
TACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAG  >  NC_000913/2549477‑2549762

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