Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 96,089 C→A 5.2% L30M (CTG→ATG)  mraY → phospho‑N‑acetylmuramoyl‑pentapeptide transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091396,0890CA 5.2% 327.6 / 5.7 105L30M (CTG→ATG) mraYphospho‑N‑acetylmuramoyl‑pentapeptide transferase
Reads supporting (aligned to +/- strand):  ref base C (42/57);  new base A (2/4);  total (44/61)
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.80e-01

Only 100 of 113 total aligned reads displayed.

GTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGATTATGAT  >  NC_000913/95940‑96237
                                                                                                                                                     |                                                                                                                                                    
gTTAATGGTTCACGTAGTGCCGACATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTTTTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGcc                                                                                                                                                      <  2:804990/150‑1 (MQ=255)
   aaGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc                                                                                                                                                   <  1:729020/150‑1 (MQ=255)
   aaGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgct                                                                                                                                                  >  2:404121/1‑150 (MQ=255)
   aaGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGGTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc                                                                                                                                                   <  1:62451/150‑1 (MQ=255)
     gggTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgctg                                                                                                                                                 <  2:354264/150‑1 (MQ=255)
     gggTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCaggctg                                                                                                                                                 <  2:761675/150‑1 (MQ=255)
      ggTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGa                                                                                                                                                >  1:1012188/1‑150 (MQ=255)
         tCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcg                                                                                                                                             <  1:1499622/150‑1 (MQ=255)
          cACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc                                                                                                                                            >  1:98576/1‑150 (MQ=255)
          cACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc                                                                                                                                            >  2:1337687/1‑150 (MQ=255)
          cACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc                                                                                                                                            <  2:1431989/150‑1 (MQ=255)
          cACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACAGTTGGTCAAATATTAGTCCGTCTTTAACGTCTTTTCCTAACTGACGTATCGCGCCATCGTCAGCCTGCTGACcgc                                                                                                                                            >  1:1577020/1‑150 (MQ=255)
            cGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCgcgc                                                                                                                                          <  1:593384/150‑1 (MQ=255)
            cGTAGAGCAGCAATGGAAGAGGTATTACGCGCTTTACAGGAGAATGGTACATGTTAGTTTGGCTGTCCGAACATTTTGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCTTCAGCAGGCTGAcgatgg                                                                                                                                          <  1:865346/150‑6 (MQ=255)
              tAGTGCCGCCATGGTAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCAGGCTGAcgatggcg                                                                                                                                        <  1:595782/150‑8 (MQ=255)
               aGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGTCGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGt                                                                                                                                       <  1:380058/150‑1 (MQ=255)
                gTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGtt                                                                                                                                      >  2:530950/1‑150 (MQ=255)
                  gccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCa                                                                                                                                    >  1:1245374/1‑150 (MQ=255)
                  gccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCAGGCTGAcgatg                                                                                                                                           <  2:1701939/143‑5 (MQ=255)
                  gccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCAGGCTGAcgatg                                                                                                                                           >  1:1701939/1‑139 (MQ=255)
                   ccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAt                                                                                                                                   <  2:711361/150‑1 (MQ=255)
                     gccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAtct                                                                                                                                 >  1:977623/1‑150 (MQ=255)
                     gccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAtct                                                                                                                                 <  2:686694/150‑1 (MQ=255)
                     gccATGGAAGAGGTAGTACGCGCTTTACACGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTACTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAtct                                                                                                                                 >  1:30012/1‑150 (MQ=255)
                     gccATGGAAGAGGTAGTACGAGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAtct                                                                                                                                 <  2:649133/150‑1 (MQ=255)
                      ccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCAtctc                                                                                                                                <  2:943560/150‑1 (MQ=255)
                        aTGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCAt                                                                                                                              <  1:274417/150‑1 (MQ=255)
                         tGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGCCCGCGCTGTTCATCTCAtt                                                                                                                             >  1:1172750/1‑150 (MQ=255)
                              gagGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGa                                                                                                                        <  1:652476/150‑1 (MQ=255)
                                    gtaCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGcc                                                                                                                  >  1:525092/1‑150 (MQ=255)
                                     taCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGccc                                                                                                                 <  1:1272249/150‑1 (MQ=255)
                                     taCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGGCAAATATTATTCCGGCTGTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCGGACCGCGCTGTTCATCTCATTGTGGATGGGccc                                                                                                                 <  1:343560/150‑1 (MQ=255)
                                      aCGCGCTTTACAGGAGAATGGTACATTTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTACCCTCATTACCTATCTGACCTATCGAGCCATCGTCAGCCTCCTTACCGCGCTGTTCAGCTGCTGGTCGATGGGCCcg                                                                                                                >  2:675755/1‑150 (MQ=255)
                                      aCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCcg                                                                                                                >  2:962618/1‑150 (MQ=255)
                                         cgcTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGt                                                                                                             >  2:1343753/1‑150 (MQ=255)
                                           cTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTAt                                                                                                           <  2:1118336/150‑1 (MQ=255)
                                                cAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTg                                                                                                      <  2:1147727/150‑1 (MQ=255)
                                                cAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTg                                                                                                      <  2:300585/150‑1 (MQ=255)
                                                 aGGAGACGGGGACAAGTTAGTTTGGCTGGCCGAACATTTGGTCAAATAGTCTTCCGGCTGTACCGGCTTTTCCTCTCTGACGTTTCGCGCCATCGTCAGCCTGCCGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGc                                                                                                     <  1:214085/150‑1 (MQ=255)
                                                  ggAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCATCCTCCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGAGTATGATTGCt                                                                                                    >  1:70398/1‑150 (MQ=255)
                                                      aaTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAtt                                                                                                <  1:356440/150‑1 (MQ=255)
                                                       aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAttt                                                                                               <  1:408285/150‑1 (MQ=255)
                                                          ggACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCACATTGTGCATGGGCCCGCGTATGATTGCTCATTTGCa                                                                                            >  2:581569/1‑150 (MQ=255)
                                                               tGTTAGTTTGGCTGGCCGCACATTTGCTCAAAAATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAAc                                                                                       <  2:1172750/150‑1 (MQ=255)
                                                               tGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAAc                                                                                       <  1:519770/150‑1 (MQ=255)
                                                                   aGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTc                                                                                   <  2:30012/150‑1 (MQ=255)
                                                                        ggctggcCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTATCCTTTg                                                                              >  2:206300/1‑150 (MQ=255)
                                                                         gctggcCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTgg                                                                             >  1:747679/1‑150 (MQ=255)
                                                                          ctggcCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGt                                                                            <  1:800484/150‑1 (MQ=255)
                                                                           tggcCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTc                                                                           >  2:1254400/1‑150 (MQ=255)
                                                                           tggcCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTATCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTc                                                                           <  2:1403334/150‑1 (MQ=255)
                                                                            ggcCGAACATTTGGGCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCa                                                                          <  1:1569829/150‑1 (MQ=255)
                                                                               cGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGAGCTGGTCATCTCATAGTGGAAGGGCCCGCGTATGATTGATCATTTGCAAAAACTTTCCTTTGGTCAggt                                                                       >  1:1196294/1‑150 (MQ=255)
                                                                                gAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAggtg                                                                      >  1:1730650/1‑150 (MQ=255)
                                                                                   cATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCGGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGgtggtg                                                                   <  2:1187655/150‑1 (MQ=255)
                                                                                    aTTTGGGCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCTCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGc                                                                  <  1:1327482/150‑1 (MQ=255)
                                                                                       tGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTa                                                               <  2:207422/150‑1 (MQ=255)
                                                                                        ggTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTaa                                                              >  1:1262939/1‑150 (MQ=255)
                                                                                         gTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGTCGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAac                                                             <  1:157129/150‑1 (MQ=255)
                                                                                         gTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAac                                                             >  1:222305/1‑150 (MQ=255)
                                                                                         gTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGACTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCAGAGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAac                                                             >  1:1646919/1‑150 (MQ=255)
                                                                                          tCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAacg                                                            <  2:626785/150‑1 (MQ=255)
                                                                                           cAAATAATATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATCGCCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGAGTAacga                                                           >  2:29147/1‑150 (MQ=255)
                                                                                            aaaTATTATTCCTGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGGCAGCCTGCTGCCCGCGCTGTTCATCTCATGGGGGATGGGCCCGCGTATGATGGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAacgac                                                          <  1:1570863/150‑1 (MQ=255)
                                                                                            aaaTATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGGCAGGTGGTGCGTAacgac                                                          <  1:1065520/150‑1 (MQ=255)
                                                                                             aaTATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTACTTTGGTCAGGTGGTGCGTAacgacg                                                         >  2:787315/1‑150 (MQ=255)
                                                                                              atctaattCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACgg                                                        <  2:958834/145‑1 (MQ=255)
                                                                                               tattatTCCGCCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCCGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGt                                                       <  2:239542/150‑1 (MQ=255)
                                                                                                attattCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTc                                                      <  2:1150044/150‑1 (MQ=255)
                                                                                                 ttattCCGGCTTTAACGTCTTTTCCTATCTGAGGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTAATTGTGTATGGGCCCGCGTATGATTGCTCATTTTCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTcc                                                     <  1:437665/150‑1 (MQ=255)
                                                                                                 ttattCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTcc                                                     <  1:1430923/150‑1 (MQ=255)
                                                                                                  tattCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCt                                                    <  1:588828/150‑1 (MQ=255)
                                                                                                   tttCCGGCGTTACCGGCTTTTCCTATCTGACGTGTCGCGCCATCGGCAGCCGGCTGACCGCGCTGTTCATCTCATTGTGGAGGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTg                                                   <  2:263390/149‑1 (MQ=255)
                                                                                                    ttCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGa                                                  >  1:271397/1‑150 (MQ=255)
                                                                                                      ccGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAAt                                                >  1:1035769/1‑150 (MQ=255)
                                                                                                      ccGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAAt                                                >  1:206017/1‑150 (MQ=255)
                                                                                                      ccGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAAt                                                >  2:1437835/1‑150 (MQ=255)
                                                                                                      ccGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGGTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAAt                                                <  1:1431973/150‑1 (MQ=255)
                                                                                                              aaCGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTc                                        <  2:486880/150‑1 (MQ=255)
                                                                                                                cGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAg                                      <  2:1035769/150‑1 (MQ=255)
                                                                                                                cGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAg                                      <  2:203067/150‑1 (MQ=255)
                                                                                                                cGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAg                                      >  1:173963/1‑150 (MQ=255)
                                                                                                                        ccTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGcgcg                              >  1:1483011/1‑150 (MQ=255)
                                                                                                                        ccTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGAAGGTCCTGAATCACACTTCAGCAAGcgcg                              >  1:805806/1‑150 (MQ=255)
                                                                                                                         cTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCGGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCgg                             <  1:1195427/150‑1 (MQ=255)
                                                                                                                          tctCTGACTGTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTTTGATTGCTCATTTGCAAAAACTTTCCTATGGTATGGTGGTGCGGAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGt                            <  2:1500053/148‑1 (MQ=255)
                                                                                                                              tGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATGTCATTGTGGAAGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGc                        <  1:1081011/150‑1 (MQ=255)
                                                                                                                              tGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGc                        <  1:845406/150‑1 (MQ=255)
                                                                                                                                 cGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGa                     <  1:1141064/150‑1 (MQ=255)
                                                                                                                                       gccccATCGTCACCCAGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCCTATGATTGATCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGAAACGACGGTCCTGAATAACACTGCAGCAAGCGCGGAACCCCGCCCATgg               >  2:686854/4‑150 (MQ=255)
                                                                                                                                        cgcCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACTGTCCTGAATCACACATCAGCAAGCGCGGTACGCCGACCATggg              >  2:919305/1‑150 (MQ=255)
                                                                                                                                         gcCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGc             >  1:1383967/1‑150 (MQ=255)
                                                                                                                                         gcCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGc             <  1:1107581/150‑1 (MQ=255)
                                                                                                                                           cATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCgg           <  2:1245374/150‑1 (MQ=255)
                                                                                                                                            aTCGTCAGCCTGCTGACCGCGCTGTTCATCACATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCggg          <  2:337122/150‑1 (MQ=255)
                                                                                                                                              cGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGAt        <  1:749803/150‑1 (MQ=255)
                                                                                                                                               gTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGAtt       <  1:1343753/150‑1 (MQ=255)
                                                                                                                                                  aGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGATTATg    >  1:117726/1‑150 (MQ=255)
                                                                                                                                                   gCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGATTATGa   >  2:451830/1‑150 (MQ=255)
                                                                                                                                                    ccTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGATTATGAt  <  2:222305/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCATTTGCAAAAACTTTCCTTTGGTCAGGTGGTGCGTAACGACGGTCCTGAATCACACTTCAGCAAGCGCGGTACGCCGACCATGGGCGGGATTATGAT  >  NC_000913/95940‑96237

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