Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 350,953 G→T 6.9% L121L (CTG→CTT hisH → imidazole glycerol phosphate synthase subunit HisH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947350,9530GT 6.9% 226.0 / 2.5 71L121L (CTG→CTThisHimidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base T (3/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCAGGACCGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGTGGTGGGTCGCGGTCA  >  NC_002947/350806‑351099
                                                                                                                                                   |                                                                                                                                                  
ccAGGACCGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAGAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGaa                                                                                                                                                  <  1:791373/150‑1 (MQ=255)
ccAGGACCGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGaa                                                                                                                                                  <  1:212175/150‑1 (MQ=255)
       cGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGc                                                                                                                                             >  1:1044549/1‑148 (MQ=255)
       cGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGc                                                                                                                                             <  2:1044549/148‑1 (MQ=255)
       cGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCg                                                                                                                                           <  2:24521/150‑1 (MQ=255)
         cccGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGca                                                                                                                                         >  2:1218962/1‑150 (MQ=255)
            gTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACAt                                                                                                                                      >  2:895034/1‑150 (MQ=255)
            gTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACAt                                                                                                                                      >  1:1138874/1‑150 (MQ=255)
               ccTCGGCATCTGTGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATggg                                                                                                                                   <  1:115896/150‑1 (MQ=255)
                 tCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCt                                                                                                                                 >  2:418900/1‑150 (MQ=255)
                     cATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGaa                                                                                                                             <  1:565712/150‑1 (MQ=255)
                      aTCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAAc                                                                                                                            <  1:895034/150‑1 (MQ=255)
                      aTCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAAc                                                                                                                            <  1:829921/150‑1 (MQ=255)
                         tGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGaa                                                                                                                         <  2:98820/150‑1 (MQ=255)
                             tCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCa                                                                                                                     <  1:1299723/150‑1 (MQ=255)
                                gTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCt                                                                                                                                 >  2:899201/1‑135 (MQ=255)
                                gTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCt                                                                                                                                 <  1:899201/135‑1 (MQ=255)
                                 tATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCa                                                                                                                 >  2:436941/1‑150 (MQ=255)
                                     cAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGGCACCTTCAGGTGctcgc                                                                                                                                         <  1:981114/121‑5 (MQ=255)
                                     cAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGGCACCTTCAGGTGctcgc                                                                                                                                         >  2:981114/1‑117 (MQ=255)
                                              ctCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGc                                                                                                                                             >  2:1218149/1‑109 (MQ=255)
                                              ctCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGc                                                                                                                                             <  1:1218149/109‑1 (MQ=255)
                                               tCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGaccacc                                                                                                   >  2:52122/1‑150 (MQ=255)
                                                     acaGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCtgt                                                                                             <  1:429102/150‑1 (MQ=255)
                                                     acaGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCtgt                                                                                             <  1:1218962/150‑1 (MQ=255)
                                                       agggAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTgg                                                                                           >  1:1425850/4‑150 (MQ=255)
                                                        gTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGc                                                                                          >  2:1351127/1‑150 (MQ=255)
                                                        gTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGc                                                                                          >  1:714678/1‑150 (MQ=255)
                                                         tGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCa                                                                                                                     >  1:366921/1‑122 (MQ=255)
                                                         tGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCa                                                                                                                     <  2:366921/122‑1 (MQ=255)
                                                            ggAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGa                                                                                      <  1:941286/150‑1 (MQ=255)
                                                                         gTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATcgc                                                                         <  2:1138874/150‑1 (MQ=255)
                                                                            gACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGa                                                                                                        <  2:1048630/116‑1 (MQ=255)
                                                                            gACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGa                                                                                                        >  1:1048630/1‑116 (MQ=255)
                                                                               tGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGt                                                                   >  1:283372/1‑150 (MQ=255)
                                                                                gCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGtt                                                                  <  1:716319/150‑1 (MQ=255)
                                                                                  aTTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTc                                                                >  1:155275/1‑150 (MQ=255)
                                                                                  aTTGGCCTGTTCCCTGGCCAGGTGCGCTTTTCCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATggg                                                                                                                                   >  2:1278508/1‑83 (MQ=255)
                                                                                  aTTGGCCTGTTCCCTGGCCAGGTGCGCTTTTCCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATggg                                                                                                                                   <  1:1278508/83‑1 (MQ=255)
                                                                                        cTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACAACATCCCTGATCGCGCACGTTTCTACTTc                                                          >  1:1307533/1‑150 (MQ=255)
                                                                                         tGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGAc                                                                                     <  2:1368037/122‑1 (MQ=255)
                                                                                         tGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGAc                                                                                     >  1:1368037/1‑122 (MQ=255)
                                                                                          gTTCCCTGGCCAGGGGCGCTTTTTCGTCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGt                                                        >  2:766929/1‑150 (MQ=255)
                                                                                            tCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGc                                                      <  1:52122/150‑1 (MQ=255)
                                                                                             cccTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGca                                                     >  2:1304607/1‑150 (MQ=255)
                                                                                              ccTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGGAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGcac                                                    >  1:1153007/1‑150 (MQ=255)
                                                                                               cTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGcaca                                                   <  2:1026759/150‑1 (MQ=255)
                                                                                                  gCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCtgcgg                                                                                           >  1:1409641/1‑104 (MQ=255)
                                                                                                  gCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGCgg                                                                                           <  2:1409641/107‑1 (MQ=255)
                                                                                                      ggTGCGCTGTGTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGTAACGAAGTCAGCCAGACCATCGACCACCCGCTGTCGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGctacta                                            <  2:404080/150‑1 (MQ=255)
                                                                                                       gTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCtactac                                           <  1:1351127/150‑1 (MQ=255)
                                                                                                          cgcTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATc                                        <  2:1056188/150‑1 (MQ=255)
                                                                                                           gcTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCtactac                                           >  2:1191822/1‑146 (MQ=255)
                                                                                                           gcTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCtactac                                           <  1:1191822/146‑1 (MQ=255)
                                                                                                                    gCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGAAACGTACCGGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACCGCTACTACATCAATGCCGGCa                              >  1:798853/1‑150 (MQ=255)
                                                                                                                     cAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCaa                             <  1:1036010/150‑1 (MQ=255)
                                                                                                                      aaaGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAg                            >  2:1121339/1‑150 (MQ=255)
                                                                                                                         gACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCg                         <  2:1425850/150‑1 (MQ=255)
                                                                                                                          aCCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCgg                        <  2:714678/150‑1 (MQ=255)
                                                                                                                               cAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGAt                                                                            <  2:213620/93‑1 (MQ=255)
                                                                                                                               cAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGAt                                                                            >  1:213620/1‑93 (MQ=255)
                                                                                                                                    agagGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAggtggt              <  1:1088425/150‑1 (MQ=255)
                                                                                                                                           cggCACCTTCAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGTGGTGGGTcgc       >  1:519970/3‑150 (MQ=255)
                                                                                                                                            aGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGTGGTGGGTcgcg      <  1:418900/150‑1 (MQ=255)
                                                                                                                                             gCACCTTCAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGtggtgg            <  2:1216862/143‑1 (MQ=255)
                                                                                                                                             gCACCTTCAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGtggtgg            >  1:1216862/1‑143 (MQ=255)
                                                                                                                                                ccTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGTGGTGGGTCGCGGTCa  >  2:134483/1‑150 (MQ=255)
                                                                                                                                                 cTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAg                                                  <  2:602783/101‑1 (MQ=255)
                                                                                                                                                 cTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAg                                                  >  1:602783/1‑101 (MQ=255)
                                                                                                                                                 cTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCACCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACGTCAATGc                                   >  2:1019696/1‑116 (MQ=255)
                                                                                                                                                 cTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCACCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACGTCAATGc                                   <  1:1019696/116‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGGACCGCCCGTTCCTCGGCATCTGCGTCGGTATGCAGGCGCTGCTCGAGCACAGTGAGGAAAACGCCGGCGTCGACTGCATTGGCCTGTTCCCTGGCCAGGTGCGCTTTTTCGGCAAAGACCTGCAAGAAGAGGGCGAGCACCTGAAGGTGCCGCACATGGGCTGGAACGAAGTCAGCCAGACCATCGACCACCCGCTGTGGCACGACATCCCTGATCGCGCACGTTTCTACTTCGTGCACAGCTACTACATCAATGCCGGCAAGCCGGGCCAGGTGGTGGGTCGCGGTCA  >  NC_002947/350806‑351099

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