Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 99,898 T→C 6.9% R85R (CGT→CGC murG → N‑acetylglucosaminyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091399,8980TC 6.9% 257.2 / 7.7 89R85R (CGT→CGCmurGN‑acetylglucosaminyl transferase
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base C (3/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GTTCGCTGGCTGGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCA  >  NC_000913/99749‑100043
                                                                                                                                                     |                                                                                                                                                 
gTTCGCTGGCTGGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCGGCGTATCTTCAACGCCTGGCGt                                                                                                                                                   >  2:1023384/1‑150 (MQ=255)
    gctggctgGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAgg                                                                                                                                               <  1:32344/150‑1 (MQ=255)
        gctggGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGACGCCAGgcgtt                                                                                                                                           <  2:850225/150‑3 (MQ=255)
        gctggGGACTGCCGACCGTAAGGGAGCGGACTTAGTGCCAAAACATTGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGGTGCCCCGCTGCGTATCTTCAACTCCTGACGCCAGgcgtt                                                                                                                                           <  1:1239845/150‑3 (MQ=255)
         ctggGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCgg                                                                                                                                          >  1:1526230/1‑150 (MQ=255)
          tggGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCggg                                                                                                                                         <  1:69233/150‑1 (MQ=255)
           ggggACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGc                                                                                                                                        >  2:1245086/1‑150 (MQ=255)
            gggACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCg                                                                                                                                       <  1:236136/150‑1 (MQ=255)
             ggACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGa                                                                                                                                      >  1:1413379/1‑150 (MQ=255)
               aCTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGAtt                                                                                                                                    <  2:1209837/150‑1 (MQ=255)
                cTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTa                                                                                                                                   <  1:181950/150‑1 (MQ=255)
                 tGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTAt                                                                                                                                  <  2:1539090/150‑1 (MQ=255)
                  gCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTTTAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATg                                                                                                                                 >  2:755042/1‑150 (MQ=255)
                  gCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATg                                                                                                                                 >  2:1251488/1‑150 (MQ=255)
                  gCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATg                                                                                                                                 <  2:228009/150‑1 (MQ=255)
                   ccGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGa                                                                                                                                >  2:1281702/1‑150 (MQ=255)
                   ccGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGa                                                                                                                                >  2:83302/1‑150 (MQ=255)
                        cGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCg                                                                                                                           <  1:310419/150‑1 (MQ=255)
                        cGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCg                                                                                                                           <  2:51462/150‑1 (MQ=255)
                         gTATGGAAGAGGACTTAGTGCCAAAACATGGAATCGAAATTGATTTCATTCGTATCTATGTTCTGCGTGGAAAAGGTATAAAAGCACTGATAGATGACCAGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGt                                                                                                                          >  2:1546090/1‑150 (MQ=255)
                           aTGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTAc                                                                                                                        >  2:1518757/1‑150 (MQ=255)
                             ggAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACaa                                                                                                                      <  1:78862/150‑1 (MQ=255)
                               aaGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAAc                                                                                                                    >  1:944078/1‑150 (MQ=255)
                                  cGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTg                                                                                                                 <  1:83302/150‑1 (MQ=255)
                                  cGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTg                                                                                                                 >  1:369319/1‑150 (MQ=255)
                                        tAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtgg                                                                                                           >  1:414918/1‑150 (MQ=255)
                                        tAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACgtgg                                                                                                           <  1:1606360/150‑1 (MQ=255)
                                           tGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGc                                                                                                        >  2:1620680/1‑150 (MQ=255)
                                            gCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCt                                                                                                       >  2:212813/1‑150 (MQ=255)
                                            gCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCt                                                                                                       >  1:412462/1‑150 (MQ=255)
                                                aaaCATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGt                                                                                                   <  1:353157/150‑1 (MQ=255)
                                                    aTGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATgg                                                                                               <  1:1620680/150‑1 (MQ=255)
                                                     tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATggg                                                                                              <  1:1263714/150‑1 (MQ=255)
                                                     tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATggg                                                                                              <  1:1245086/150‑1 (MQ=255)
                                                     tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATggg                                                                                              <  1:738208/150‑1 (MQ=255)
                                                     tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGAGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATggg                                                                                              <  2:94056/150‑1 (MQ=255)
                                                     tGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGGATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGGCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATggg                                                                                              <  2:88301/150‑1 (MQ=255)
                                                       gCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAg                                                                                            <  2:1526230/150‑1 (MQ=255)
                                                        cATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAgg                                                                                           >  1:169312/1‑150 (MQ=255)
                                                         aTCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGc                                                                                          <  2:360461/150‑1 (MQ=255)
                                                             aaaTTGATTTCATTCGTATCTCTGGGCTGCGGGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACg                                                                                      <  2:911483/150‑1 (MQ=255)
                                                                   aTTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAg                                                                                <  1:268241/150‑1 (MQ=255)
                                                                     ttCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTCACGTGGTGCTCGGTATGGGAGGCTACGTGTCAgct                                                                              >  2:540368/1‑148 (MQ=255)
                                                                       cATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTcc                                                                            <  1:1139919/150‑1 (MQ=255)
                                                                             tATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAACGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAggtgg                                                                      <  1:1239110/150‑1 (MQ=255)
                                                                             tATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAggtgg                                                                      >  2:157053/1‑150 (MQ=255)
                                                                             tATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCGGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAggtgg                                                                      <  2:1110138/150‑1 (MQ=255)
                                                                               tctcTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTc                                                                    >  1:339484/1‑150 (MQ=255)
                                                                                    ggTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTTCTCGGTATGGTAGGCTACGTGTCAGGTCCAGGTGGTATGGcc                                                               >  2:770502/1‑150 (MQ=255)
                                                                                    ggTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGcc                                                               >  1:631062/1‑150 (MQ=255)
                                                                                    ggTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGGCCAGGTGGTCTGGcc                                                               <  2:555890/150‑1 (MQ=255)
                                                                                        tGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGt                                                           >  2:359908/1‑150 (MQ=255)
                                                                                        tGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGt                                                           >  1:1072022/1‑150 (MQ=255)
                                                                                         gCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTg                                                          >  2:1409814/1‑150 (MQ=255)
                                                                                           gTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGt                                                        >  2:1595265/1‑150 (MQ=255)
                                                                                             ggAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCg                                                      >  2:556505/1‑150 (MQ=255)
                                                                                                aaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTa                                                   <  1:648531/150‑1 (MQ=255)
                                                                                                aaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTa                                                   >  2:1141285/1‑150 (MQ=255)
                                                                                                aaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTa                                                   <  1:1518757/150‑1 (MQ=255)
                                                                                                aaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTa                                                   >  2:1135636/1‑150 (MQ=255)
                                                                                                aaaGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTa                                                   <  1:212813/150‑1 (MQ=255)
                                                                                                  aGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAgg                                                 <  1:1251488/150‑1 (MQ=255)
                                                                                                         aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCg                                          >  2:643909/1‑150 (MQ=255)
                                                                                                         aaaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCg                                          <  2:1148018/150‑1 (MQ=255)
                                                                                                          aaGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCgg                                         <  2:788380/150‑1 (MQ=255)
                                                                                                           aGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGt                                        <  1:556505/150‑1 (MQ=255)
                                                                                                             cACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCTATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTg                                      <  1:1290721/150‑1 (MQ=255)
                                                                                                                 gATAGCTGTCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACt                                  >  2:1437810/1‑150 (MQ=255)
                                                                                                                    aGCTGTCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCa                               <  1:1437810/150‑1 (MQ=255)
                                                                                                                     gCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGAACTTCAt                              <  2:1224036/150‑1 (MQ=255)
                                                                                                                        gCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGaa                           >  1:966624/1‑150 (MQ=255)
                                                                                                                        gCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCGACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGaa                           >  2:822026/1‑150 (MQ=255)
                                                                                                                         ccccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCa                               <  1:1381028/145‑1 (MQ=255)
                                                                                                                         ccccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCa                               >  2:1381028/1‑145 (MQ=255)
                                                                                                                         ccccGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAAc                          >  1:506021/1‑150 (MQ=255)
                                                                                                                              cTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTTGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAAc                     <  1:1409814/150‑1 (MQ=255)
                                                                                                                                gCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACgg                   <  2:169312/150‑1 (MQ=255)
                                                                                                                                  gTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTa                 >  1:164244/1‑150 (MQ=255)
                                                                                                                                     tCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTg              >  1:1615504/1‑150 (MQ=255)
                                                                                                                                      cTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGc             >  2:976629/1‑150 (MQ=255)
                                                                                                                                       ttCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCg            >  1:1285964/1‑150 (MQ=255)
                                                                                                                                       ttCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCg            >  1:523491/1‑150 (MQ=255)
                                                                                                                                       gcccgcGCCTGACGCCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTg                                      >  1:1352702/6‑124 (MQ=255)
                                                                                                                                       gcccgcGCCTGACGCCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTg                                      <  2:1352702/119‑1 (MQ=255)
                                                                                                                                       gcccgcGCCTGACGCCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCg            >  1:876014/6‑150 (MQ=255)
                                                                                                                                        tCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCgg           <  2:1413379/150‑1 (MQ=255)
                                                                                                                                        cccgcGCCTGACGCCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCgg           >  2:350895/5‑150 (MQ=255)
                                                                                                                                               cTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAAc    <  2:1311936/150‑1 (MQ=255)
                                                                                                                                                tGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAAcc   <  2:412462/150‑1 (MQ=255)
                                                                                                                                                 ggCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCa  >  1:312526/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTCGCTGGCTGGGGACTGCCGACCGTATGGAAGCGGACTTAGTGCCAAAACATGGCATCGAAATTGATTTCATTCGTATCTCTGGTCTGCGTGGAAAAGGTATAAAAGCACTGATAGCTGCCCCGCTGCGTATCTTCAACGCCTGGCGTCAGGCGCGGGCGATTATGAAAGCGTACAAACCTGACGTGGTGCTCGGTATGGGAGGCTACGTGTCAGGTCCAGGTGGTCTGGCCGCGTGGTCGTTAGGCATTCCGGTTGTACTTCATGAACAAAACGGTATTGCGGGCTTAACCA  >  NC_000913/99749‑100043

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