Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,608 T→C 7.7% intergenic (+202/+38) 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,6080TC 7.7% 341.6 / 14.0 104intergenic (+202/+38)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (49/47);  new base C (6/2);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 104 total aligned reads displayed.

TGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGC  >  NC_000913/2549465‑2549755
                                                                                                                                               |                                                                                                                                                   
tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTAc                                                                                                                                                <  2:1551275/149‑1 (MQ=255)
            tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                                                                                        >  1:1662753/1‑145 (MQ=255)
            tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                                                                                        <  2:1662753/145‑1 (MQ=255)
              cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTaa                                                                                                                                  <  4:1086059/149‑1 (MQ=255)
               gCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                                                                                                                 >  4:108238/1‑149 (MQ=255)
                 ttATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCa                                                                                                                               >  4:1464003/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                                                                                                                                     >  2:1817793/1‑124 (MQ=255)
                                    cccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                                                                                                                     <  1:1817793/124‑1 (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)
                                                        cACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAg                                                                                        >  1:1852173/1‑149 (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                                                                                 >  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                                                                                                                                <  4:81336/93‑1 (MQ=255)
                                                                        ggATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATc                                                                                                                                >  3:81336/1‑93 (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                                                                                                                                           >  3:1540104/3‑72 (MQ=38)
                                                                                  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)
                                                                                   acACGCCGGATGCGGCATGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                      <  1:1520762/139‑1 (MQ=255)
                                                                                     aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGt                                                                                                                                                   >  1:1578538/1‑61 (MQ=255)
                                                                                     aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGt                                                                                                                                                   <  2:1578538/61‑1 (MQ=255)
                                                                                     aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                           >  1:868569/1‑149 (MQ=255)
                                                                                     aCGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTa                                                           >  2:60555/1‑149 (MQ=255)
                                                                                        ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCg                                                        <  4:1492530/149‑1 (MQ=255)
                                                                                        ccgcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCg                                                        >  2:1379256/1‑149 (MQ=255)
                                                                                          gcATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGt                                                      >  2:167325/1‑149 (MQ=255)
                                                                                             tCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTg                                                   <  2:1852173/149‑1 (MQ=255)
                                                                                                        cAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCaa                                        <  1:1379256/149‑1 (MQ=255)
                                                                                                         aaCGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                       >  1:1823243/1‑149 (MQ=255)
                                                                                                         aaCGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAAc                                       <  2:868569/149‑1 (MQ=255)
                                                                                                            gcgcACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATc                                    >  4:35701/1‑149 (MQ=255)
                                                                                                                 cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGa                                                                                                                           <  3:1703954/57‑1 (MQ=255)
                                                                                                                 cTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGa                                                                                                                           >  4:1703954/1‑57 (MQ=255)
                                                                                                                     gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                        >  1:625091/1‑72 (MQ=255)
                                                                                                                     gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcg                                                                                                        <  2:625091/72‑1 (MQ=255)
                                                                                                                     gTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGc                           <  1:167325/149‑1 (MQ=255)
                                                                                                                      tCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTGGTGAAACGCAACATCGCATCACGCtt                         >  1:590668/1‑149 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCATCTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                              <  2:593924/111‑1 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                              <  2:1149141/79‑1 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                                                                                              >  1:1149141/1‑79 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                      <  1:21133/103‑1 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                      >  2:21133/1‑103 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                              >  1:593924/1‑111 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                 <  1:1670166/124‑1 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                 >  2:1670166/1‑124 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                            <  4:772742/129‑1 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAAc                                            >  3:772742/1‑129 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                        >  3:786099/1‑149 (MQ=255)
                                                                                                                        aacaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTa                        >  3:1663097/1‑149 (MQ=255)
                                                                                                                             tCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                      >  4:1771941/1‑98 (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)
                                                                                                                                   gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                 >  4:1176962/1‑113 (MQ=39)
                                                                                                                                   gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGt                                                 <  3:1176962/113‑1 (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                                                                      >  2:1293439/1‑85 (MQ=255)
                                                                                                                                          agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAataata                                                                      <  1:1293439/85‑1 (MQ=255)
                                                                                                                                          agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                     >  4:712663/1‑86 (MQ=255)
                                                                                                                                          agtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg                                                                     <  3:712663/86‑1 (MQ=255)
                                                                                                                                          agtAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAAc      <  2:1201250/149‑1 (MQ=255)
                                                                                                                                           gtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                               >  1:948472/1‑91 (MQ=255)
                                                                                                                                           gtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg                                                               <  2:948472/91‑1 (MQ=255)
                                                                                                                                            tcaTGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTTCAGCGACGCTGCGAAATAATACCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCACACAt    >  1:728563/3‑149 (MQ=255)
                                                                                                                                            tAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACa                                      <  1:642115/115‑1 (MQ=255)
                                                                                                                                            tAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACa                                      >  2:642115/1‑115 (MQ=255)
                                                                                                                                              aTGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGc  >  1:1027773/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGC  >  NC_000913/2549465‑2549755

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