Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 100,158 C→T 5.3% P172L (CCG→CTG)  murG → N‑acetylglucosaminyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913100,1580CT 5.3% 429.8 / 9.4 132P172L (CCG→CTG) murGN‑acetylglucosaminyl transferase
Reads supporting (aligned to +/- strand):  ref base C (72/52);  new base T (4/3);  total (76/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 132 total aligned reads displayed.

TTCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCAAAG  >  NC_000913/100011‑100304
                                                                                                                                                    |                                                                                                                                                  
ttCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCg                                                                                                                                                   >  4:400876/1‑149 (MQ=255)
ttCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCg                                                                                                                                                   >  1:664381/1‑149 (MQ=255)
ttCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCg                                                                                                                                                   >  3:814743/1‑149 (MQ=255)
 tCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCgc                                                                                                                                                  >  3:350050/1‑149 (MQ=255)
   aTGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCgcag                                                                                                                                                >  1:713221/1‑149 (MQ=255)
     gAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCgcagca                                                                                                                                              >  2:65922/1‑149 (MQ=255)
      aaCAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCaa                                                                                                                                             >  1:118045/1‑149 (MQ=255)
      aaCAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCaa                                                                                                                                             >  4:333734/1‑149 (MQ=255)
      aaCAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCaa                                                                                                                                             >  2:404716/1‑149 (MQ=255)
         aaaaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGt                                                                                                                                          >  3:849937/1‑149 (MQ=255)
          aaaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGtt                                                                                                                                         >  4:519107/1‑149 (MQ=255)
          aaaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGtt                                                                                                                                         <  4:814743/149‑1 (MQ=255)
          aaaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGtt                                                                                                                                         <  1:308079/149‑1 (MQ=255)
           aaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCGGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCGGCAACtttt                                                                                                                                        >  3:821013/1‑149 (MQ=255)
           aaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCTGCGGCAACggca                                                                                                                                        <  1:387221/149‑4 (MQ=255)
           aaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCTGCGGCAACggca                                                                                                                                        <  4:214593/149‑4 (MQ=255)
              ggTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGc                                                                                                                                     <  1:64518/149‑1 (MQ=255)
              ggTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGc                                                                                                                                     <  4:1015957/149‑1 (MQ=255)
               gTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    <  3:400876/149‑1 (MQ=255)
               gTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    >  2:272764/1‑149 (MQ=255)
                        gCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGaa                                                                                                                           <  4:849937/149‑1 (MQ=255)
                          ttAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAgg                                                                                                                         >  2:257535/1‑149 (MQ=255)
                          tcaaCCAATAAAAGGCTGGCGAAGATTGCCACCAACGTGATGCAGGCGATTCCAGGTGCTTGCCCGACTGAGGCAGTAGTGGATAACGCGGGGCGTAGCGATGTGCCGGCGCTGCAGTGGCCGCAGCAACGTTTGGCTGGACGTGAAgg                                                                                                                         >  4:485104/3‑149 (MQ=255)
                            aaCCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTc                                                                                                                       >  2:969579/1‑149 (MQ=255)
                              ccAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCg                                                                                                                     >  1:914343/1‑149 (MQ=255)
                                aaTAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGt                                                                                                                   <  4:224233/149‑1 (MQ=255)
                                aaTAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGt                                                                                                                   >  1:444952/1‑149 (MQ=255)
                                aaTAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTATCCCTAATGCGGAAGAAGTGGGTAACCCGGTGCGTACCCATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGt                                                                                                                   >  3:154727/1‑149 (MQ=255)
                                  tAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTc                                                                                                                 <  4:257163/149‑1 (MQ=255)
                                       ggctggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGtgtg                                                                                                            >  3:924651/1‑149 (MQ=255)
                                        gctggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGc                                                                                                           <  1:501058/149‑1 (MQ=255)
                                          tggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTg                                                                                                         >  3:627420/1‑149 (MQ=255)
                                           ggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTgg                                                                                                        <  2:118045/149‑1 (MQ=255)
                                             cGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTa                                                                                                      <  4:537550/149‑1 (MQ=255)
                                              gAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAg                                                                                                     >  2:469672/1‑149 (MQ=255)
                                              gAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAg                                                                                                     >  2:629303/1‑149 (MQ=255)
                                                 gATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgg                                                                                                  >  4:423113/1‑149 (MQ=255)
                                                    tGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgggtg                                                                                               <  1:404716/149‑1 (MQ=255)
                                                       caccaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGtt                                                                                            >  4:317603/1‑149 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           >  3:828122/1‑149 (MQ=255)
                                                         ccaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGCGGTtct                                                                                          >  2:406210/1‑149 (MQ=255)
                                                          caAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtctc                                                                                         <  2:917134/149‑1 (MQ=255)
                                                           aaaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCa                                                                                        <  1:226114/149‑1 (MQ=255)
                                                           aaaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGAGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCa                                                                                        <  2:819233/149‑1 (MQ=255)
                                                            aaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    >  3:363571/1‑104 (MQ=255)
                                                            aaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    <  4:363571/104‑1 (MQ=255)
                                                             aGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAgg                                                                                      <  1:257535/149‑1 (MQ=255)
                                                               tGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGgc                                                                                    >  3:117373/1‑149 (MQ=255)
                                                                  tGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGcgca                                                                                 >  2:284715/1‑149 (MQ=255)
                                                                     aGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGcgcacgc                                                                              >  1:727613/1‑149 (MQ=255)
                                                                     aGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGcgcacgc                                                                              >  4:897880/1‑149 (MQ=255)
                                                                     aGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACAGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGcgcacgc                                                                              >  2:914525/1‑149 (MQ=255)
                                                                          tttCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCt                                                                         <  1:1020490/149‑1 (MQ=255)
                                                                                ggTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTACTCAGGGCGCACGCATTCTTAACCa                                                                   <  3:317603/149‑1 (MQ=255)
                                                                                  tGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAAc                                                                     >  4:214202/1‑145 (MQ=255)
                                                                                  tGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGa                                                                 <  4:41600/149‑1 (MQ=255)
                                                                                  tGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGa                                                                 >  2:199664/1‑149 (MQ=255)
                                                                                     tttCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACaa                                                              >  2:826196/1‑149 (MQ=255)
                                                                                       tCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATg                                                            <  3:519107/149‑1 (MQ=255)
                                                                                         ccTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgggtgg                                                                                              <  1:471197/113‑1 (MQ=255)
                                                                                          cTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTGGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCg                                                         >  1:25842/1‑149 (MQ=255)
                                                                                           tAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGc                                                        >  3:929169/1‑149 (MQ=255)
                                                                                            aaTGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCa                                                       >  4:798141/1‑149 (MQ=255)
                                                                                              tGCGGACGTAGTGGTTACCCGGGGCGGTCACGAGGTGTTGGAGCTGCCGTTGCGGATGCCATGTTTGGGTGGACGTGCAGGTCCGGTTCGTGTGCGGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTACCCAGACAATGCCGCAgg                                                     <  2:818379/149‑1 (MQ=255)
                                                                                                 ggAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtg                                                  >  3:229768/1‑149 (MQ=255)
                                                                                                 ggAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtg                                                  >  4:783151/1‑149 (MQ=255)
                                                                                                 ggAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtg                                                  >  3:569959/1‑149 (MQ=255)
                                                                                                  gAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgc                                                 >  1:237910/1‑149 (MQ=255)
                                                                                                   aaGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgct                                                <  4:117373/149‑1 (MQ=255)
                                                                                                         tGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGaaa                                          >  4:617098/1‑149 (MQ=255)
                                                                                                            gTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTg                                       >  1:604620/1‑149 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      >  3:401160/1‑149 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  2:598028/149‑1 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  2:161974/149‑1 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGAGGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  4:821013/149‑1 (MQ=255)
                                                                                                              aaCCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTggg                                     <  1:629303/149‑1 (MQ=255)
                                                                                                               aCCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGt                                    <  4:431139/149‑1 (MQ=255)
                                                                                                               aCCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGt                                    <  3:897880/149‑1 (MQ=255)
                                                                                                                 ccGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGa                                  >  4:722096/1‑149 (MQ=255)
                                                                                                                  cGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGAt                                 >  1:830018/1‑149 (MQ=255)
                                                                                                                   ggTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGAtt                                <  3:798141/149‑1 (MQ=255)
                                                                                                                   agTGCGTACCGATGTGTTGGCGCTGCCGGTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAAATGGGTGAtt                                <  4:902324/148‑1 (MQ=255)
                                                                                                                    gTGCGTACCGATGTGTTGGAGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTc                               >  2:914827/1‑149 (MQ=255)
                                                                                                                         gTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTc                                                                                                                       >  4:586267/1‑57 (MQ=255)
                                                                                                                         gTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTc                           >  3:334316/1‑149 (MQ=255)
                                                                                                                               aTGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGt                                    >  4:354992/1‑134 (MQ=255)
                                                                                                                               aTGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGt                                    <  3:354992/134‑1 (MQ=255)
                                                                                                                               aTGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATc                     >  4:663689/1‑149 (MQ=255)
                                                                                                                                tgtgTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCt                    <  4:591310/149‑1 (MQ=255)
                                                                                                                                    ttGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCa                >  3:38930/1‑149 (MQ=255)
                                                                                                                                     tGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCAt               <  3:671398/149‑1 (MQ=255)
                                                                                                                                         gcTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCaga           >  3:432048/1‑149 (MQ=255)
                                                                                                                                           aaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  4:262551/117‑1 (MQ=255)
                                                                                                                                           aaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      >  3:262551/4‑120 (MQ=255)
                                                                                                                                            gCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCg        <  3:663689/149‑1 (MQ=255)
                                                                                                                                             acGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCgg       >  4:102914/2‑149 (MQ=255)
                                                                                                                                             acGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCgg       >  1:1053400/2‑149 (MQ=255)
                                                                                                                                               gTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCa     >  4:12411/1‑149 (MQ=255)
                                                                                                                                               gTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCa     >  2:278748/1‑149 (MQ=255)
                                                                                                                                               gTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCa     >  1:782992/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCAAAG  >  NC_000913/100011‑100304

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