Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,549 C→A 15.9% intergenic (+143/+97) 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,5490CA15.9% 262.1 / 35.0 82intergenic (+143/+97)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base A (4/9);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

AATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATA  >  NC_000913/2549405‑2549687
                                                                                                                                                |                                                                                                                                          
aataaAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                        <  2:1678706/149‑1 (MQ=255)
     aaaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                   <  3:1469586/149‑5 (MQ=37)
      aaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                  <  4:1877683/149‑6 (MQ=25)
        ggCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                >  3:686889/1‑149 (MQ=255)
            agtcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                                                                                           <  4:1259019/134‑1 (MQ=255)
            agtcagtAATGTTACTGGCTGGTCTGAGAGTGCTAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTAGTAGGCCGGATAAGGCGTTCa                                                                                                                                           >  3:1259019/1‑134 (MQ=255)
             gtcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTACGCCGGATGCGGCGTGAACGCCTTATCCGGccta                                                                                                                           <  4:215901/149‑4 (MQ=25)
              tcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtgt                                                                                                                          >  3:1086059/1‑149 (MQ=255)
                 gtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGAAGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACa                                                                                                                       >  2:557395/1‑149 (MQ=255)
                   aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAc                                                                                                                     <  1:557395/149‑1 (MQ=255)
                   aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTACGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                            >  1:1005022/1‑140 (MQ=21)
                   aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTACGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                            <  2:964920/142‑3 (MQ=21)
                   aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTACGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                            <  2:1005022/142‑3 (MQ=21)
                   aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTACGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                            >  1:964920/1‑140 (MQ=25)
                               tggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                   <  4:937863/123‑5 (MQ=21)
                               tggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                   >  3:937863/1‑119 (MQ=21)
                                   cTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGt                                                                                                          >  3:566082/1‑144 (MQ=255)
                                   cTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGt                                                                                                          <  4:566082/144‑1 (MQ=255)
                                    tGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaacaa                                                                                                    <  4:686889/149‑1 (MQ=255)
                                                 aaaCGCAAAACTGCCTGATGCGTTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtc                                                                                       >  1:1551275/1‑149 (MQ=255)
                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                              <  4:1240617/110‑1 (MQ=18)
                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                              >  3:1240617/1‑110 (MQ=18)
                                                  aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAcgc                                                                                                                   <  3:1834979/120‑1 (MQ=255)
                                                  aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAcgc                                                                                                                   >  4:1834979/1‑120 (MQ=255)
                                                     gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                             <  4:744010/107‑2 (MQ=18)
                                                            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)
                                                                               aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                <  2:1486608/78‑1 (MQ=255)
                                                                               aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                >  1:1486608/1‑78 (MQ=255)
                                                                               aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                         <  3:108238/149‑1 (MQ=255)
                                                                                      cTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                                                                                           <  1:349471/60‑1 (MQ=255)
                                                                                      cTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                                                                                           >  2:349471/1‑60 (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)
                                                                                                   acaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTaacat                                                                               <  4:1151658/107‑4 (MQ=255)
                                                                                                   acaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTaacat                                                                               >  3:1151658/1‑104 (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)
                                                                                                            tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTaacatta                                                                             >  2:1212070/1‑95 (MQ=255)
                                                                                                            tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTaacatta                                                                             <  1:1212070/100‑6 (MQ=255)
                                                                                                              aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtt                          >  3:1492530/1‑149 (MQ=255)
                                                                                                               aTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                                               <  1:1105688/79‑1 (MQ=255)
                                                                                                               aTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                                               >  2:1105688/1‑79 (MQ=255)
                                                                                                                  aGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                            <  2:247173/79‑1 (MQ=255)
                                                                                                                  aGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                            >  1:247173/1‑79 (MQ=255)
                                                                                                                   ggaCGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTc                                                                                                         >  4:1362782/3‑65 (MQ=255)
                                                                                                                   gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTc                                                                                                         <  3:1362782/65‑1 (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)
                                                                                                                             gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                                               >  1:956925/1‑65 (MQ=255)
                                                                                                                             gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                                               <  2:956925/65‑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                                                           <  1:539653/97‑1 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCa                                                           >  2:539653/1‑97 (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                                                                  <  4:556892/78‑1 (MQ=255)
                                                                                                                                             gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCt                                                                  >  3:556892/1‑78 (MQ=255)
                                                                                                                                              tacACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCt                                                                       <  4:1540104/70‑1 (MQ=38)
                                                                                                                                              tacACGCCGGATGCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCt                                                                       >  3:1540104/3‑72 (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)
                                                                                                                                                |                                                                                                                                          
AATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATA  >  NC_000913/2549405‑2549687

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