Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,761 A→C 7.2% intergenic (+189/‑251) 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,7610AC 7.2% 189.5 / 3.6 76intergenic (+189/‑251)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base C (2/4);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA  >  NC_000913/349614‑349909
                                                                                                                                                   |                                                                                                                                                    
gTCGCTGACGCGACTTATCAGGCCTACGAGGGGCACTGAACTGTAGGTCGTATAAGACGGATGGCGTCGCATCCGACAACCGATGCCGGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTGAACTGTAGGTCGGATAAGa                                                                                                                                                    <  1:1564039/150‑1 (MQ=255)
 tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGACCTGTAGGTCGGATAATACGGATGGCTTCGCATACGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                                                                                                                   <  1:1113063/150‑1 (MQ=255)
 tCGCTGACGCGACTTATCAGGCCAACGAGGCGCACTGAACTGTAGGCCGGATAAGACCGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                                                                                                                   <  2:1266546/150‑1 (MQ=255)
   gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGc                                                                                                                                                 <  2:484257/150‑1 (MQ=255)
      gACGCGACTTATCAGGCCTACGAGGGGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                                                                                                              <  1:1485699/150‑1 (MQ=255)
       aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGctatg                                                                                                                                             <  2:598840/150‑5 (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)
                                         tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                                     <  2:138343/108‑1 (MQ=40)
                                         tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                                     >  1:138343/1‑108 (MQ=40)
                                            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:379507/1‑150 (MQ=32)
                                                                                                                     tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                               <  1:1344027/150‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
GTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA  >  NC_000913/349614‑349909

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: