Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 837,635 G→A 5.4% intergenic (+199/+30) hcxB → / ← ybiJ hydroxycarboxylate dehydrogenase B/DUF1471 domain‑containing protein YbiJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913837,6350GA 5.4% 295.9 / 4.1 92intergenic (+199/+30)hcxB/ybiJhydroxycarboxylate dehydrogenase B/DUF1471 domain‑containing protein YbiJ
Reads supporting (aligned to +/- strand):  ref base G (37/50);  new base A (2/3);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AATGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGA  >  NC_000913/837486‑837783
                                                                                                                                                     |                                                                                                                                                    
aaTGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATg                                                                                                                                                      >  2:505192/1‑150 (MQ=255)
    aGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCCGGCAACAGATGCCTGATgcgac                                                                                                                                                  <  2:59590/150‑3 (MQ=38)
      ggTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGt                                                                                                                                                <  2:97240/150‑1 (MQ=255)
      ggTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATAAGGCGt                                                                                                                                                >  2:271173/1‑150 (MQ=255)
        ttAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTaa                                                                                                                                              >  1:995966/1‑150 (MQ=255)
         tAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTaaa                                                                                                                                             >  2:952368/1‑150 (MQ=255)
         tAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTaaa                                                                                                                                             >  2:460928/1‑150 (MQ=255)
           gagCATATGCGTGTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACg                                                                                                                                           >  1:1608144/1‑150 (MQ=255)
           gagCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACg                                                                                                                                           <  1:385674/150‑1 (MQ=255)
            agCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                          <  1:1492140/150‑1 (MQ=255)
             gCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                         <  2:961606/150‑1 (MQ=255)
             gCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                         >  1:449195/1‑150 (MQ=255)
             gCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGGAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                         <  2:1437407/150‑1 (MQ=255)
               atatGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                       <  1:505192/150‑1 (MQ=255)
                 atGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                     <  1:1124058/150‑1 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    >  2:783091/1‑150 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    >  2:808244/1‑150 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    <  2:1204698/150‑1 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    >  1:1596146/1‑150 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    <  1:79481/150‑1 (MQ=255)
                  tGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                    <  1:386161/150‑1 (MQ=255)
                    cGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                  <  1:1425903/150‑1 (MQ=255)
                     gTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGt                                                                                                                                 >  1:132332/1‑150 (MQ=255)
                     gTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                 <  1:1041577/150‑2 (MQ=38)
                        tctCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCt                                                                                                                              >  1:819016/3‑150 (MQ=255)
                        tGTCGGCAAACATACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                              <  1:1417096/150‑1 (MQ=255)
                        tGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATAAGGCGTTTACGCCGCATCCggcaa                                                                                                                              <  1:1178610/150‑5 (MQ=37)
                          tCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATAAGGCGTTTACGCCGCATCCggcaaca                                                                                                                            <  1:271173/150‑7 (MQ=37)
                               aaaCAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAtt                                                                                                                       <  1:1345777/150‑1 (MQ=255)
                                  cagacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACt                                                                                                                    <  2:687746/150‑1 (MQ=255)
                                    gacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTta                                                                                                                  >  1:118059/1‑150 (MQ=255)
                                    gacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTta                                                                                                                  <  1:433818/150‑1 (MQ=255)
                                      cagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTtata                                                                                                                <  2:1332674/150‑1 (MQ=255)
                                        gggAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATaaa                                                                                                              >  1:809719/1‑150 (MQ=255)
                                               aCTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGc                                                                                                       <  1:865401/150‑1 (MQ=255)
                                               aCTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGc                                                                                                       <  1:575637/150‑1 (MQ=255)
                                                cTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                      >  1:689182/1‑150 (MQ=255)
                                                    tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTAc                                                                                                  <  1:1019070/150‑1 (MQ=255)
                                                    tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTAc                                                                                                  <  2:112739/150‑1 (MQ=255)
                                                     gCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTAcc                                                                                                 >  1:19679/1‑150 (MQ=255)
                                                         gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCtt                                                                                             <  2:19679/150‑1 (MQ=255)
                                                           cGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTaa                                                                                           >  1:1359650/1‑150 (MQ=255)
                                                             tAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCtt                                                                                             <  2:870411/146‑1 (MQ=255)
                                                             tAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCtt                                                                                             >  1:870411/1‑146 (MQ=255)
                                                                 gTAACCGCCTGATCCGTCCATCGGTTCGGATTTTGTCGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                     <  1:710722/150‑1 (MQ=255)
                                                                  tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttg                                                                                    <  1:321578/150‑1 (MQ=255)
                                                                  tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttg                                                                                    <  2:1186543/150‑1 (MQ=255)
                                                                  tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttg                                                                                    <  2:718763/150‑1 (MQ=255)
                                                                    aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgtt                                                                                  <  2:436665/150‑1 (MQ=255)
                                                                      cGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgt                                                                                >  1:975839/1‑150 (MQ=255)
                                                                        ccTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgttg                                                                              >  2:1093476/1‑150 (MQ=255)
                                                                          tGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggt                                                                            <  1:808244/150‑1 (MQ=255)
                                                                             tCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGc                                                                         >  2:1462549/1‑150 (MQ=255)
                                                                               cGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCgg                                                                       >  1:379436/1‑150 (MQ=255)
                                                                               cGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCgg                                                                       >  2:69123/1‑150 (MQ=255)
                                                                                gTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGa                                                                      <  1:673206/150‑1 (MQ=255)
                                                                                gTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGa                                                                      <  2:742781/150‑1 (MQ=255)
                                                                                  ccAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                    <  2:1608144/150‑1 (MQ=255)
                                                                                   cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGt                                                                   >  2:527491/1‑150 (MQ=255)
                                                                                   cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGt                                                                   >  2:802323/1‑150 (MQ=255)
                                                                                    aaCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTg                                                                  >  1:1439317/1‑150 (MQ=255)
                                                                                     aCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGa                                                                 <  1:1336375/150‑1 (MQ=255)
                                                                                      cGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGAt                                                                <  2:305186/150‑1 (MQ=255)
                                                                                       ggTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATg                                                               <  1:1542833/150‑1 (MQ=255)
                                                                                       ggTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCAGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATg                                                               >  1:1491246/1‑150 (MQ=255)
                                                                                         ttCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCt                                                             <  1:460928/150‑1 (MQ=255)
                                                                                          tCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTg                                                            <  2:132332/150‑1 (MQ=255)
                                                                                          tCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTg                                                            <  1:1615168/150‑1 (MQ=255)
                                                                                           cGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGATTTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGGGGAAGTGATGCTGt                                                           <  2:1577563/150‑1 (MQ=255)
                                                                                           cGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGt                                                           <  2:889422/150‑1 (MQ=255)
                                                                                            ggATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTa                                                          >  1:644754/1‑150 (MQ=255)
                                                                                            ggATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTa                                                          <  2:360486/150‑1 (MQ=255)
                                                                                            ggATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTa                                                          >  2:272692/1‑150 (MQ=255)
                                                                                            ggATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTa                                                          <  1:977263/150‑1 (MQ=255)
                                                                                               ttttGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCAACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATcc                                                       <  1:1168042/150‑1 (MQ=255)
                                                                                                   gTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGTATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTg                                                   >  1:192924/1‑150 (MQ=255)
                                                                                                        ccTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCAcc                                              >  1:741248/1‑150 (MQ=255)
                                                                                                         cTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCg                                             <  2:689462/150‑1 (MQ=255)
                                                                                                              aaGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGc                                        >  1:1389014/1‑150 (MQ=255)
                                                                                                              aaGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGc                                        >  1:61618/1‑150 (MQ=255)
                                                                                                                  cgcgGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCcgcg                                    <  2:572045/150‑1 (MQ=255)
                                                                                                                        aGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTtc                              <  1:502477/150‑1 (MQ=255)
                                                                                                                        aGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTtc                              <  1:783091/150‑1 (MQ=255)
                                                                                                                           gTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCa                           >  2:647731/1‑150 (MQ=255)
                                                                                                                            tCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAg                          <  2:1596146/150‑1 (MQ=255)
                                                                                                                               cATCAGGCATCTGTTGCCGGATACGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCa                       >  1:31658/1‑150 (MQ=255)
                                                                                                                                  cAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGtt                    >  1:1604685/1‑150 (MQ=255)
                                                                                                                                       aTCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCt               <  1:944699/150‑1 (MQ=255)
                                                                                                                                       aTCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCt               <  2:819016/150‑1 (MQ=255)
                                                                                                                                         cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTc             >  2:276409/1‑150 (MQ=255)
                                                                                                                                          tGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTcc            >  1:503702/1‑150 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAggg        <  1:527491/150‑1 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTAAAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAggg        >  2:812192/1‑150 (MQ=255)
                                                                                                                                                    tGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGa  >  2:879456/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATGAGGGTTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGA  >  NC_000913/837486‑837783

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