Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,768 T→A 5.9% intergenic (+196/‑244) 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,7680TA 5.9% 187.4 / 2.6 74intergenic (+196/‑244)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base A (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

GACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCC  >  NC_000913/349620‑349916
                                                                                                                                                    |                                                                                                                                                    
gACGCGACTTATCAGGCCTACGAGGGGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                                                                                                                     <  1:1485699/150‑1 (MQ=255)
   gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                                                                                                                  <  2:1286451/150‑1 (MQ=40)
      aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                                                                               <  1:1593013/150‑1 (MQ=39)
           tCAGTCCTACGATGTGCAGTGAACTGGAGGACGCTTAACACGGAGTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                                                          <  1:85049/150‑1 (MQ=11)
           tCAGGCCTACGAGGTGCACTGAACTGTAGTTCGGATAATACGGATTGCGTCGTATCAGACAACCGACGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTAGGATAAGACGCATAGCGTCGCATcc                                                                                                                                          <  2:1013224/150‑1 (MQ=11)
              ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAc                                                                                                                                       <  2:648803/150‑1 (MQ=14)
              ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                                                       <  2:49639/150‑1 (MQ=39)
               gCCTACGAGGTGCACTGACCGGTCGGGCGTAAAAGCCGGATGGCGGCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAcc                                                                                                                                      <  1:432402/150‑2 (MQ=14)
                 cTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                                                                                                                    <  2:390352/150‑1 (MQ=39)
                  atcGAGGTGCACTCAACTGTAGGTCGGATACGACGGAGGGCGTCGCATCCGTCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                                                                                                                   <  2:369491/148‑1 (MQ=255)
                          gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                                                           <  2:974667/150‑1 (MQ=37)
                             cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg                                                                                                                        >  2:9520/1‑150 (MQ=37)
                             cTGAACTGTAGGTCGGATAAGACGGAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg                                                                                                                        <  2:1349032/150‑1 (MQ=37)
                               gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcc                                                                                                                      >  1:458904/1‑149 (MQ=37)
                               gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg                                                                                                                      >  2:496502/1‑150 (MQ=21)
                                 aCTGTAGGTCGTATAAGACGGATGTCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                                    <  2:308234/150‑1 (MQ=255)
                                 aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                                    >  2:1344027/1‑150 (MQ=37)
                                 aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                                    <  2:534365/150‑1 (MQ=17)
                                  cTGTAGGTCGCATAACACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                                                                   <  1:835337/150‑1 (MQ=255)
                                      aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                               >  2:1550198/1‑150 (MQ=37)
                                       ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATCAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                                                                                                              >  1:912989/1‑150 (MQ=37)
                                        gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATAAGGCATACCGTCTTCAATGAACTGTAGGTCTTATAATACGCATAGCGTCGCATCCGACAACCGATTCCTTATGCGCCGCTGAcg                                                                                                             >  2:337235/1‑150 (MQ=11)
                                          cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                                                                           <  1:967921/150‑1 (MQ=37)
                                           ggATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                                                          >  1:725223/1‑150 (MQ=37)
                                                    ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg                                                                                                 >  1:1600203/1‑150 (MQ=35)
                                                           gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCTGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                          <  2:981112/150‑1 (MQ=17)
                                                           gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                          <  2:1096812/150‑1 (MQ=32)
                                                               cATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTg                                                                                      >  1:1500470/1‑150 (MQ=32)
                                                                    taCAACCGATGCCTGAGGCGCCGCTGAAGTGACTTATCAGGCATACGGGGTTCACTGAACTGTTGGTCGTAAAAGACGCATAGGGTCGCATGCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                                 <  2:147277/149‑1 (MQ=11)
                                                                    gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGTATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                                 <  1:1550198/150‑1 (MQ=21)
                                                                       aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                              <  2:1220678/150‑1 (MQ=34)
                                                                                 tGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                    <  2:459981/150‑1 (MQ=21)
                                                                                  gATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATTCCTGATTCGCCGCTGACGCGACTAATCAGGCCTACGGGGTGCACGGAACTGTAAGTCGGa                                                                   >  2:833087/1‑150 (MQ=11)
                                                                                   aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                                                                                                                            >  2:1706102/1‑92 (MQ=33)
                                                                                    tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                                 >  2:1398037/1‑150 (MQ=34)
                                                                                        ccgcTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCCACTTATCAGACCTACGGCGGGCACTGAACTGTAGGTCGCATAAGAc                                                             >  1:96079/1‑150 (MQ=11)
                                                                                              aCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACTCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGCCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                                       >  2:1319255/1‑150 (MQ=21)
                                                                                               cGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGCCTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTGGGATAAAACCCCTAGc                                                      >  1:1658541/1‑150 (MQ=11)
                                                                                               cGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCAGAAGCCCCGCTGACTCGACTAATAAGCCAAACGGCGAGCACTAAACAGAAGGACAGAAAAGACCCATAGc                                                      >  2:1042245/1‑150 (MQ=11)
                                                                                                      tATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATTCGCCGCTGACGCGACTTATAAGGCCAACGGGGTGCACTGAACTTTAGGTCGCATAAGACGCATACCGTCGCaa                                               >  1:1256706/1‑149 (MQ=11)
                                                                                                      tATCACGCCGACGGGGTGCCCTCACCTGGTGGTCTGATAAGCCGCATCGGCTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                               <  1:300936/150‑1 (MQ=11)
                                                                                                       aTCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                                              >  2:735839/1‑150 (MQ=32)
                                                                                                        tCAGGCCTACGGGGGGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                             <  2:1410886/150‑1 (MQ=17)
                                                                                                         cAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                            >  2:446855/1‑150 (MQ=32)
                                                                                                          aGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                           <  1:9520/150‑1 (MQ=32)
                                                                                                           ggCCTACGGGGTGCACTGAACTGTAGTTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                          <  2:833335/150‑1 (MQ=17)
                                                                                                            gCCTCCGGGGCGCACTGAACTGTCGGGCGCACAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAcc                                         <  1:131625/150‑2 (MQ=11)
                                                                                                               tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                      <  1:1344027/150‑1 (MQ=32)
                                                                                                               tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                      >  1:379507/1‑150 (MQ=32)
                                                                                                                aCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAAGCATAGCGTCGCATCCGACAAccc                                     >  2:1616461/1‑149 (MQ=17)
                                                                                                                   gggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                  >  2:409454/1‑150 (MQ=35)
                                                                                                                     gTGCTCTGACCTGTAGGTCGCATAAGACGCATTGGGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                <  2:640894/150‑1 (MQ=11)
                                                                                                                       gCACTGTACTGTAGGTCGTATAATACGCATTGCGTAGCATCCGACAACAGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                              <  2:201599/150‑1 (MQ=11)
                                                                                                                         aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCTTCGCATCCGACAACAGATGCCAGAt                            >  2:1069986/1‑150 (MQ=255)
                                                                                                                            gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg                         <  2:1600203/150‑1 (MQ=37)
                                                                                                                            gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATccg                         >  1:1456363/1‑150 (MQ=255)
                                                                                                                              ccTGTAGGTCGGATAAGACGCATCGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                       <  1:998592/149‑1 (MQ=21)
                                                                                                                               cTGTCGGTCTGATAAGACGCATTGTTTCGCATCCGACAACCGATGCCTGATGCGACGCTGACGCGTCTGATCTTGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                      <  2:282139/150‑1 (MQ=11)
                                                                                                                                 gTAGGTCGGATAAGACGCATAGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                    <  2:55274/150‑1 (MQ=37)
                                                                                                                                     gTCGGATGCGACGCTATGCGTCTTATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATTCCTGATGCGCCGCTGAcc                >  2:92277/1‑149 (MQ=11)
                                                                                                                                     gTCGGATGCGACGCTATGCGTCTTATACGACAACCCATGCCAGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGGGCACAGAACTGTAGGTCGGATAAGACGCAAAGCGGCGCCTCCGACAACCCATGCCACATGCGCCCCTGAcc                >  2:673279/1‑149 (MQ=2)
                                                                                                                                     gTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg                >  2:1418170/1‑150 (MQ=37)
                                                                                                                                     gTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcc                >  2:75851/1‑149 (MQ=37)
                                                                                                                                     gTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGAATCCTCATGCGCCGCTGAcg                >  2:1418174/1‑150 (MQ=255)
                                                                                                                                     atCGGATAAGACGCATAGCGTCGCATCCGACAACCAATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg                >  1:157154/2‑150 (MQ=37)
                                                                                                                                      ccGGATAAGACGCATCCCCTCTAATCCGAAAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTTCGGGTTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgc               <  2:1163000/149‑1 (MQ=11)
                                                                                                                                         gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc            >  1:1034030/1‑150 (MQ=37)
                                                                                                                                         caTACGACTCATCGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc            <  1:270856/149‑1 (MQ=255)
                                                                                                                                          aTACGACGCACCGCGTCGCTTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTCTCACGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt           <  1:337235/150‑1 (MQ=11)
                                                                                                                                            aaGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTa         <  1:496502/150‑1 (MQ=17)
                                                                                                                                              gACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc       >  1:1702062/1‑150 (MQ=37)
                                                                                                                                                 gCATAGCGTCGAATCCGACAACCGATGCCTGATGCGCCGCTGACGCTACTTATCAGGCCTACGGGGTGCTCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg    <  2:912989/150‑1 (MQ=255)
                                                                                                                                                   aTAGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc  <  1:1559298/150‑1 (MQ=21)
                                                                                                                                                   aTAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc  <  1:1398037/150‑1 (MQ=34)
                                                                                                                                                    |                                                                                                                                                    
GACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCC  >  NC_000913/349620‑349916

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