Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 100,158 C→T 5.8% 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.8% 320.5 / 5.2 98P172L (CCG→CTG) murGN‑acetylglucosaminyl transferase
Reads supporting (aligned to +/- strand):  ref base C (47/45);  new base T (4/2);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCAAAG  >  NC_000913/100011‑100304
                                                                                                                                                   |                                                                                                                                                  
ttCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCg                                                                                                                                                   >  3:449774/1‑149 (MQ=255)
ttCATGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCg                                                                                                                                                   >  4:681944/1‑149 (MQ=255)
   aTGAACAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCgcag                                                                                                                                                >  3:465205/1‑149 (MQ=255)
       aCAAAACGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAAc                                                                                                                                            <  4:532807/149‑1 (MQ=255)
          aaaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGTGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGtt                                                                                                                                         <  4:600057/149‑1 (MQ=255)
           aaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGttt                                                                                                                                        >  2:80755/1‑149 (MQ=255)
           aaCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGttt                                                                                                                                        <  1:614017/149‑1 (MQ=255)
            aCGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTATTGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTg                                                                                                                                       >  1:96931/1‑149 (MQ=255)
             cGGTATTGCGGGCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTgg                                                                                                                                      <  1:735040/149‑1 (MQ=255)
                        gCTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGAGATTTCCCTAATGAGGAAGAAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTAGTCGTGaa                                                                                                                           >  4:591482/1‑149 (MQ=255)
                         cTTAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAg                                                                                                                          <  1:403297/149‑1 (MQ=255)
                          ttAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAgg                                                                                                                         >  1:546642/1‑149 (MQ=255)
                           tAACCAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGt                                                                                                                        <  3:489328/149‑1 (MQ=255)
                               cAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCgg                                                                                                                    >  3:451945/1‑149 (MQ=255)
                               cAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCgg                                                                                                                    >  1:653939/1‑149 (MQ=255)
                               cAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCgg                                                                                                                    >  3:237606/1‑149 (MQ=255)
                               cAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCgg                                                                                                                    >  3:198981/1‑149 (MQ=255)
                               cAATAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCgg                                                                                                                    >  1:252779/1‑149 (MQ=255)
                                aaTAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGt                                                                                                                   >  1:641376/1‑149 (MQ=255)
                                  tAAATGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTc                                                                                                                 <  3:32972/149‑1 (MQ=255)
                                   aaaTGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCg                                                                                                                <  3:667361/149‑1 (MQ=255)
                                   aaaGGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCg                                                                                                                <  1:80755/149‑1 (MQ=255)
                                    aaTGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCgt                                                                                                               <  1:499256/149‑1 (MQ=255)
                                     aTGGCTGGCGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCgtg                                                                                                              >  2:302129/1‑149 (MQ=255)
                                        gctggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGc                                                                                                           <  4:48042/149‑1 (MQ=255)
                                        gctggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGc                                                                                                           >  3:513217/1‑149 (MQ=255)
                                        gctggcGAAGATTGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGc                                                                                                           <  4:744141/149‑1 (MQ=255)
                                                    tGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgggtg                                                                                               <  2:333981/149‑1 (MQ=255)
                                                    tGCCACCAAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgggtg                                                                                               <  1:112581/149‑1 (MQ=255)
                                                    tGCCACCAAAGTGATGCAGGCGTTTACAGGTGCTTTCCCTAATGCGGCAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtgggtg                                                                                               <  4:449774/149‑1 (MQ=255)
                                                      ccaccaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGggtggt                                                                                             >  2:55336/1‑149 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTTATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           >  4:551317/1‑149 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           <  2:653939/149‑1 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           >  1:671686/1‑149 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           >  3:97019/1‑149 (MQ=255)
                                                        accaAAGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTtc                                                                                           <  4:237606/149‑1 (MQ=255)
                                                           aaaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCa                                                                                        <  2:563848/149‑1 (MQ=255)
                                                            aaGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAg                                                                                       <  1:549682/149‑1 (MQ=255)
                                                             aGTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAgg                                                                                      <  1:625867/149‑1 (MQ=255)
                                                              gTGATGCAGGCGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCTGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAggg                                                                                     >  1:427068/1‑149 (MQ=255)
                                                                    cAGGCGTTTCCAGGTGCCTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGcgcacg                                                                               >  1:512604/1‑149 (MQ=255)
                                                                        cGTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCAtt                                                                           <  2:483340/149‑1 (MQ=255)
                                                                         gTTTCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGAGGGTGGTTCTCAGGGCGCACGCATTc                                                                          <  4:192797/149‑1 (MQ=255)
                                                                          tttCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCt                                                                         <  2:252779/149‑1 (MQ=255)
                                                                            tCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    <  3:397061/88‑1 (MQ=255)
                                                                            tCCAGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCt                                                                                                                                    >  4:397061/1‑88 (MQ=255)
                                                                              cttgTGCTGTCCCTCATGCGGAAGTCGTGGGTAACCCGGTGCGTACCGATGTGCTGGCGCTTCCGTTGACGGAACAACGTTTGGTTGGACGTGGAGGTCCGCTTCGTGTGCTGGTAGTGGGTGGTAATCCGGGCACACGCCTTCTTAAc                                                                     <  2:456660/146‑1 (MQ=255)
                                                                               aGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAAcc                                                                    >  1:384190/1‑149 (MQ=255)
                                                                               aGGTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAAcc                                                                    >  3:200470/1‑149 (MQ=255)
                                                                                ggTGCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCa                                                                   <  4:200470/149‑1 (MQ=255)
                                                                                   gCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCt                                                                         <  1:39098/140‑1 (MQ=255)
                                                                                   gCTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCt                                                                         >  2:39098/1‑140 (MQ=255)
                                                                                    cTTTCCCTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACa                                                               >  3:707697/1‑149 (MQ=255)
                                                                                         ccTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGcc                                                          >  1:493993/1‑149 (MQ=255)
                                                                                         ccTAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGcc                                                          >  3:45051/1‑149 (MQ=255)
                                                                                          cTAATGCAGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTgg                                                                                                        >  3:727477/1‑102 (MQ=255)
                                                                                          cTAATGCAGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTgg                                                                                                        <  4:727477/102‑1 (MQ=255)
                                                                                           tAATGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGc                                                        <  4:465205/149‑1 (MQ=255)
                                                                                            aaTGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCa                                                       <  4:513217/149‑1 (MQ=255)
                                                                                              tGCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAgg                                                     >  2:242253/1‑149 (MQ=255)
                                                                                               gCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTc                                                                                                                 <  1:411753/88‑1 (MQ=255)
                                                                                               gCGGAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTc                                                                                                                 >  2:411753/1‑88 (MQ=255)
                                                                                                 ggAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtg                                                  >  1:729766/1‑149 (MQ=255)
                                                                                                 ggAAGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtg                                                  >  3:562025/1‑149 (MQ=255)
                                                                                                   aaGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgct                                                >  1:372959/1‑149 (MQ=255)
                                                                                                   aaGTAGTGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgct                                                <  3:230628/149‑1 (MQ=255)
                                                                                                     gtagtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgctgc                                              <  4:244556/149‑1 (MQ=255)
                                                                                                     gtagtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTtgctgc                                              >  3:797/1‑149 (MQ=255)
                                                                                                        gtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGaa                                           >  2:190910/1‑149 (MQ=255)
                                                                                                        gtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGaa                                           >  3:315196/1‑149 (MQ=255)
                                                                                                        gtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGaa                                           >  1:645111/1‑149 (MQ=255)
                                                                                                        gtGGGTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGaa                                           >  2:125249/1‑149 (MQ=255)
                                                                                                          gggTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtggg                                                                                                 <  1:757891/93‑1 (MQ=255)
                                                                                                          gggTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAgtggg                                                                                                 >  2:757891/1‑93 (MQ=255)
                                                                                                            gTAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTg                                       >  3:189963/1‑149 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  3:478345/149‑1 (MQ=255)
                                                                                                             tAACCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGAATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTgg                                      <  2:512604/149‑1 (MQ=255)
                                                                                                              aaCCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTggg                                     <  3:749874/149‑1 (MQ=255)
                                                                                                               aCCCGGTGCGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGt                                    <  2:493993/149‑1 (MQ=255)
                                                                                                                 gggttGCGTTCTGATGTGTTAGCGCTGCGGGTGCTGCAGCACCGTTTGGCTGGGCGTGAAGGTCTGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCTTTCTGAGCCAGACAATGCCCCAGGTTGCTGCGAAAGTGGTTGa                                  <  4:461651/145‑1 (MQ=255)
                                                                                                                       cGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGt                            <  1:125249/149‑1 (MQ=255)
                                                                                                                       cGTACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGt                            <  1:400294/149‑1 (MQ=255)
                                                                                                                         tACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCa                          >  2:122744/1‑149 (MQ=255)
                                                                                                                         tACCGATGTGTTGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCa                          >  1:236436/1‑149 (MQ=255)
                                                                                                                                   ttGGCGCTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCACACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCa                >  1:278910/1‑149 (MQ=255)
                                                                                                                                      gcgcTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCa             <  4:451945/149‑1 (MQ=255)
                                                                                                                                      gcgcTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGAGGGTGGTTCTCAGGGCGCACGCATTCTTACCCAGACAATGCCGTAGGTTGCTGCGAGACTGGGTGATTCAGTCACTATCTGGCATCa             <  3:591482/149‑1 (MQ=255)
                                                                                                                                        gcTGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCaga           >  4:85951/1‑149 (MQ=255)
                                                                                                                                        ccaaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCaga           >  4:205554/6‑149 (MQ=255)
                                                                                                                                          tGCCGTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGc         <  2:641376/149‑1 (MQ=255)
                                                                                                                                          aaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAAt                                                             <  3:267772/94‑1 (MQ=255)
                                                                                                                                          aaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAAt                                                             >  4:267772/4‑97 (MQ=255)
                                                                                                                                          aaacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGc         >  4:93196/4‑149 (MQ=255)
                                                                                                                                           aacGTTGCTGCGGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTAGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCg        >  2:382733/3‑149 (MQ=255)
                                                                                                                                              gTTGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCa     >  1:574712/1‑149 (MQ=255)
                                                                                                                                                tGCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCaaa   <  4:97019/149‑1 (MQ=255)
                                                                                                                                                 gCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTTCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCAAAg  >  2:423102/1‑149 (MQ=255)
                                                                                                                                                 gCCGCAGCAACGTTTGGCTGGACGTGAAGGTCCGGTTCGTGTGCTGGTAGTGGGTGGTCCTCAGGGCGCACGCATTCTTAACCAGACAATGCCGCAGGTTGCTGCGAAACTGGGTGATTCAGTCACTATCTGGCATCAGAGCGGCAAAg  <  3:551317/149‑1 (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: