Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,642 A→G 5.8% intergenic (+70/‑370) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,6420AG 5.8% 253.3 / 3.2 97intergenic (+70/‑370)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base G (6/2);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAT  >  NC_000913/349494‑349790
                                                                                                                                                    |                                                                                                                                                    
ggAGCTGTACGAAGCCATCAACTACTACCAGTACAAAGGGAAGCTCGACAACCTGATTGCCCGTAGACAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGa                                                                                                                                                      <  2:95068/148‑1 (MQ=38)
cGAGCCGTCCGAAACCACCAACTCCTCCCAGTTCGAAGATAAGCCCGACACCCTGTTCGCCCGGAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGa                                                                                                                                                      <  2:1012878/149‑1 (MQ=37)
   gCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTGGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGt                                                                                                                                                   <  2:86732/149‑1 (MQ=255)
       tACGAAAGGATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGGTTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGCGGTGCAc                                                                                                                                               >  1:946455/1‑149 (MQ=255)
       tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc                                                                                                                                               >  1:1275034/1‑149 (MQ=255)
       tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCAc                                                                                                                                               >  1:774275/1‑149 (MQ=255)
       tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCAc                                                                                                                                               >  1:1293682/1‑149 (MQ=255)
       tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTGGATTGTATTCGACAGCCGATGCCTGATGCGGCGCTGACGCGACTTATCAGGCGTACGAGGTGCAc                                                                                                                                               >  1:643167/1‑149 (MQ=255)
           aaaGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAACCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGaa                                                                                                                                           <  1:505072/149‑1 (MQ=255)
           aaaGCATCAACTACCACCAGTACGAACAGAAGCTCGACAACCTGTTTGCCCGTAGCCAAGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGaa                                                                                                                                           <  2:1050714/149‑1 (MQ=255)
                aTCAACTCCGACCAGTACGAAGAGAAGCTCGACAACCTGTTCGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCTTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTa                                                                                                                                      <  2:897375/149‑1 (MQ=255)
                 tCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGGCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCTATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGATGCACTGAACTGtgg                                                                                                                                     >  2:858323/1‑147 (MQ=255)
                 tCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAg                                                                                                                                     >  2:1533/1‑149 (MQ=255)
                    actactACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTc                                                                                                                                  >  1:185433/1‑149 (MQ=255)
                     ctactaCCCGTACGAACGGAATCTCGACACCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCCACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCg                                                                                                                                 <  2:871665/149‑1 (MQ=37)
                          cccAGTACGGAGAGACGCTCGACAGCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGTCAGGCGATGCCTGCTGCGTCGCTGACGCGACTTATCAGGCCTACGCGGTGCACTGAACTGTAGGTCGGATaa                                                                                                                            <  2:679186/148‑1 (MQ=25)
                          acCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATaa                                                                                                                            <  2:619310/149‑1 (MQ=255)
                          acCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATaa                                                                                                                            <  1:262403/149‑1 (MQ=255)
                          acCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGCTGATTGTATTCGACAGCCGATGCCTGATGCGACGCTGACGCGACTAATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATaa                                                                                                                            >  1:1069388/1‑149 (MQ=255)
                            cAGAACGAAGAGAAGCTCGACAACCTTTTTGCCTGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                          <  1:853754/149‑1 (MQ=255)
                              gTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGCCGGAGAAGACg                                                                                                                        >  2:1087164/1‑149 (MQ=255)
                               tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACgc                                                                                                                       >  2:796064/1‑148 (MQ=255)
                               tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACgg                                                                                                                       >  1:552889/1‑149 (MQ=255)
                               tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACgg                                                                                                                       >  1:1219725/1‑149 (MQ=255)
                                 cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAt                                                                                                                     <  1:858803/149‑1 (MQ=255)
                                 cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAt                                                                                                                     <  1:1094834/149‑1 (MQ=255)
                                 cGAAGAGAAGCGCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAt                                                                                                                     <  1:820220/149‑1 (MQ=255)
                                  gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATg                                                                                                                    <  1:1287542/149‑1 (MQ=255)
                                    aggcaaccTCCACAACCTGTTTCCCCGTAGCCAAGAGAAAAAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGc                                                                                                                  <  2:733635/142‑1 (MQ=25)
                                      agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGt                                                                                                                <  1:264753/149‑1 (MQ=255)
                                      agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGt                                                                                                                <  1:117217/149‑1 (MQ=255)
                                       gaAGCTCGCCAACCTGTTGGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTc                                                                                                               <  1:1089353/149‑1 (MQ=255)
                                           cTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCCCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAt                                                                                                           >  1:300246/1‑149 (MQ=255)
                                           cTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAt                                                                                                           >  1:942322/1‑149 (MQ=255)
                                           cTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAt                                                                                                           >  1:527642/1‑149 (MQ=255)
                                             caaCAACCTGTTTGCCCGTAGCCAGGTGAAAAAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCGGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATcc                                                                                                         <  1:506760/147‑1 (MQ=255)
                                               aaaaCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGa                                                                                                       <  2:873405/147‑1 (MQ=255)
                                               aCAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGa                                                                                                       <  1:435730/149‑1 (MQ=255)
                                                cAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGAc                                                                                                      <  1:96034/149‑1 (MQ=255)
                                                     tGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCg                                                                                                 <  2:943767/149‑1 (MQ=255)
                                                      gTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTTGGATAAGACTGATGGCGTCGCATCCGACAACCGa                                                                                                <  1:654944/149‑1 (MQ=255)
                                                         tGCCCGTAGCCAGGTGAACTAAAACACGCCCGTTGATTAAACACATGGACCGGTTGCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGc                                                                                             <  1:297765/149‑1 (MQ=17)
                                                            ccGTCGCCAGGTGAAATAAAAACCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAGACGGATGGCGTCGCATCCGACAACCGATGCCTGa                                                                                         <  1:452418/149‑1 (MQ=255)
                                                            ccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTg                                                                                          >  1:520774/1‑149 (MQ=255)
                                                             cGTAGCCAGGTGAACTAAAAAACGCCCGTTGATTGTATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGa                                                                                         <  1:452530/149‑1 (MQ=255)
                                                             cGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGa                                                                                         <  1:1085312/149‑1 (MQ=255)
                                                               tAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGGTGCCTGATg                                                                                       >  1:282761/1‑149 (MQ=255)
                                                               tAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATg                                                                                       >  1:431190/1‑149 (MQ=255)
                                                               tAGCCAGGTGAAATAAAAAACGCCCGTTGATTGGATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATg                                                                                       >  1:947992/1‑149 (MQ=255)
                                                                  ccACGTCAAACAAAAACCGCCCCTTGATTGGATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCCACGACGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                                    <  2:754755/149‑1 (MQ=17)
                                                                    aGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                                  >  2:1173063/1‑149 (MQ=255)
                                                                       tGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGACGCCGCCTGCGCCGCTg                                                                               >  2:12530/1‑149 (MQ=17)
                                                                       tGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGcgccgcgg                                                                               >  2:1070848/1‑147 (MQ=255)
                                                                        gaaagaaaAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATCAGACGGAAGGCGGCGCATCCGACCACCGATGCCTGATGCGCCGCTgg                                                                              >  2:1238923/1‑148 (MQ=18)
                                                                        gAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                              <  2:230869/149‑1 (MQ=255)
                                                                        gAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                              >  1:221589/1‑149 (MQ=255)
                                                                           aTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGAAAGATGCGACGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGTCGGATAAGGCAGATGGCGTCGCATCCGACAGCCGAGGCCTGATGCGCCGCTGGCGt                                                                           >  2:222212/1‑149 (MQ=17)
                                                                                    gCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                  >  1:1116612/1‑149 (MQ=255)
                                                                                    gCCCGTTGATTGGATTCGACCGCCGATGCCTGATGCGTCGCTGACGCGACTTAGCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                  <  1:69699/149‑1 (MQ=255)
                                                                                     ccggTTGATGGTATGCGACAGGCGATGCGTGATGCGAGGCGGACGCGAAGTATCAGGCGTAAGATGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAg                                                                 <  2:179097/146‑1 (MQ=14)
                                                                                     cccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAg                                                                 <  1:1204517/149‑1 (MQ=255)
                                                                                        gTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTTTCAGGcc                                                              >  1:346884/1‑149 (MQ=255)
                                                                                         ttGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCGGCCGGAACTTATCGGGCCt                                                             >  1:316623/1‑149 (MQ=17)
                                                                                           gATTGTATTCGACAGCCGATGCCTGATGCGTCGCAGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                           <  1:858323/149‑1 (MQ=255)
                                                                                             ttGGATTAGACCGCCGAGGCGTTAGGCGGCGCGGGCGCGACGTATCAGGCATCCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACgg                                                         <  1:388460/149‑1 (MQ=14)
                                                                                              tGTATTCGACAGCCGCGGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                        <  1:1049993/149‑1 (MQ=255)
                                                                                              tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGCCGTGACTTATCAGGCCTACggg                                                        >  1:67884/1‑149 (MQ=255)
                                                                                              tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACCTGACTTATCAGGCCGACggg                                                        >  1:1015972/1‑149 (MQ=255)
                                                                                                tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGt                                                      <  1:1533/149‑1 (MQ=255)
                                                                                                 aTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                     <  1:294365/149‑1 (MQ=255)
                                                                                                  ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                                                                                                                                >  1:817697/1‑73 (MQ=11)
                                                                                                  ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                                                                                                                                <  2:817697/73‑1 (MQ=11)
                                                                                                  ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGGCAACCGATGCCTGACGCGCCGCTGACGCGGCGTATCAGGCCTCCggggaga                                                    >  2:132056/1‑146 (MQ=11)
                                                                                                   tCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                   <  1:1061703/149‑1 (MQ=255)
                                                                                                     gACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGCGGTGGACt                                                 >  1:907307/1‑149 (MQ=255)
                                                                                                     gACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTCATGCACCGGGGACGCGACTTATCAGGCCTACGGGGTGGACt                                                 >  1:358347/1‑149 (MQ=14)
                                                                                                      aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATCCCTGCTGCGCCACCGACGGAGCTTACCACGCCTCCGGGGTGcccct                                                >  1:195164/1‑145 (MQ=14)
                                                                                                      aCAGCCGATGCCTGATGCGTCGCAGACGCGACTTATCAGGCCTACGAGGTGCATTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTg                                                >  1:714102/1‑149 (MQ=255)
                                                                                                         gCCGATGCCTGAGGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                             <  1:609714/149‑1 (MQ=255)
                                                                                                          ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACt                                            >  1:1121393/1‑149 (MQ=255)
                                                                                                           cGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTg                                           <  2:520774/149‑1 (MQ=255)
                                                                                                                ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGGTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTGTCAGGCCTACGGGGTGCACTGAACTGTAGGt                                      >  1:370858/1‑149 (MQ=255)
                                                                                                                ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCCACGGGGGGCACGGCACTGTAggg                                      >  2:864004/1‑148 (MQ=14)
                                                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCGGGCCTCCGGGGTGCACTGAACTGTAGGTCgg                                   >  1:549058/1‑149 (MQ=255)
                                                                                                                    aTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGCACGTCGGATCAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCGGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                  >  1:630181/1‑149 (MQ=14)
                                                                                                                     tgggTCGCTGGAGGGACGTATGAGGCCTACGGGGTGCACTGAACTGTGGGTTGGATAAGACGGATGGGGTCGCATACCACAACCGATGCCTGATGCGTCACTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                 <  1:554048/146‑1 (MQ=11)
                                                                                                                     tGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                 <  1:711100/149‑1 (MQ=255)
                                                                                                                         tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGCTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                             >  1:94558/1‑149 (MQ=255)
                                                                                                                         tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTCCGGGGTGCACTGAACTGTCGGCCGGATAAGa                             >  2:378876/1‑149 (MQ=255)
                                                                                                                         tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGa                             >  2:424598/1‑149 (MQ=255)
                                                                                                                           gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                           >  1:381877/1‑149 (MQ=255)
                                                                                                                                 gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGc                     >  2:1287077/1‑149 (MQ=12)
                                                                                                                                 gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGc                     >  2:1099183/1‑149 (MQ=12)
                                                                                                                                      ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGc                >  1:270304/1‑149 (MQ=255)
                                                                                                                                      tgaccAGGCCTCCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                <  1:1170267/145‑1 (MQ=255)
                                                                                                                                           aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg           >  1:981502/1‑149 (MQ=255)
                                                                                                                                            ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa          <  2:947992/149‑1 (MQ=17)
                                                                                                                                            ggCCTACGAGGTGCACTGAACGGGAGGGCGGCTAAGACGGATGGCGACGAAAACGACAACCGATGCCTGATGCGGCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGa          <  1:1212861/149‑1 (MQ=14)
                                                                                                                                             gCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGAc         >  2:1274171/1‑149 (MQ=255)
                                                                                                                                                    aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt  >  1:1055187/1‑149 (MQ=39)
                                                                                                                                                    |                                                                                                                                                    
CGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAT  >  NC_000913/349494‑349790

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