Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,673 G→C 5.1% intergenic (+101/‑339) 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,6730GC 5.1% 345.4 / 3.7 143intergenic (+101/‑339)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (65/66);  new base C (7/4);  total (73/70)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATC  >  NC_000913/349525‑349818
                                                                                                                                                    |                                                                                                                                                 
tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACgg                                                                                                                                                   >  1:442713/1‑149 (MQ=255)
  cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTCCGAGGTGCACTGAACTGTAGGTCGGATAAGCCGGac                                                                                                                                                 >  2:787696/1‑148 (MQ=255)
     agagAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGc                                                                                                                                              >  2:1299755/1‑149 (MQ=255)
       agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGt                                                                                                                                            >  2:585015/1‑149 (MQ=255)
       agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGt                                                                                                                                            <  1:1546861/149‑1 (MQ=255)
        gaAGCTCGACAACCTGTTTGCCCGGAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCCGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTc                                                                                                                                           >  1:1348099/1‑149 (MQ=255)
          aGCTCCACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGc                                                                                                                                         <  1:1403356/149‑1 (MQ=255)
             tCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGCTGCCTGATGCGTCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACGGGAGGTCGGATAAGACGGGTGGGGACGCATc                                                                                                                                      >  2:496857/1‑149 (MQ=25)
             tCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATc                                                                                                                                      <  1:790322/149‑1 (MQ=255)
              cGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTGCCATcc                                                                                                                                     >  2:444140/1‑149 (MQ=255)
              agACAACCTGTTTGCCCGCAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATcc                                                                                                                                     <  1:285847/148‑1 (MQ=255)
               gACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGAACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCg                                                                                                                                    >  1:1777512/1‑149 (MQ=255)
                aCAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCgc                                                                                                                                   >  1:571023/1‑148 (MQ=255)
                aCAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGa                                                                                                                                   >  1:554265/1‑149 (MQ=255)
                 cAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGAc                                                                                                                                  >  2:974344/1‑149 (MQ=255)
                 cAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGAc                                                                                                                                  >  1:430592/1‑149 (MQ=255)
                  aaCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCTCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAGGGCGTCGCATCCGACa                                                                                                                                 >  2:82189/1‑149 (MQ=255)
                  aaCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGCACTGTAGGTCGAATAAGACGGCTGGCGTCCCATCCGACa                                                                                                                                 >  1:1396790/1‑149 (MQ=255)
                   aCCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACaa                                                                                                                                <  2:174077/149‑1 (MQ=255)
                    gcTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGACGCCTGATGCTGCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTCGGTCGGATAAGACGGATGGCGTCGCATCCGACAAc                                                                                                                               <  2:1549124/148‑1 (MQ=255)
                    ccTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAAc                                                                                                                               <  1:985680/149‑1 (MQ=255)
                        tttGCCCGTAGCCAGGTGAAATAAATAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGAt                                                                                                                           >  1:906075/1‑149 (MQ=255)
                        tttGCCCGTAGCCAGGTGAAATAAAGAACGCCCGGTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATGAGACGGATGGCGTCGCACACGGCAACCGaa                                                                                                                           >  1:905218/1‑148 (MQ=21)
                        tttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGAt                                                                                                                           >  2:642418/1‑149 (MQ=255)
                         ttGCCCGTAGCCAGGTGAAATAATAAACGCCCGTTGATAGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGCTAAGACGGATGGCGTCGCATCCGACAACCgcgg                                                                                                                          >  1:973571/1‑146 (MQ=255)
                         ttGCCCGTAGCCAGGTGAAATAAAAGACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATg                                                                                                                          >  1:1080280/1‑149 (MQ=255)
                              cGTAGCCAGGTGAAATAAAAGACGCCCGTTGATTGTATTCGGCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGa                                                                                                                     <  2:862270/149‑1 (MQ=255)
                              cGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGCTAAGACGGATGGCGTCGCATCCCACAACCGATGCCGGa                                                                                                                     >  1:1350066/1‑149 (MQ=255)
                                   gcAGGTGAAAGAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                                                                <  2:98436/148‑1 (MQ=255)
                                   ccAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACACCCGATGCCTGCTgcgc                                                                                                                >  2:12761/1‑149 (MQ=255)
                                    catgTGAAATAAAAAACGCCCGTTGATCGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTCTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                               <  1:1032282/146‑1 (MQ=255)
                                    cAGGTGAAATAAAAAACGCCCGTTGGTTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCGGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                               <  2:418358/149‑1 (MQ=255)
                                      ggTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGGCTTATCAGGCCTACGGGGTGGACTGAACTGTAGGTCGGATAAGACGGATGGGGTCGCATCGGACGGCCGAGGCGTGATGcgccgg                                                                                                             >  1:2460/1‑148 (MQ=18)
                                       gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGGTGCCTGATGCGCCGCt                                                                                                            >  1:1015810/1‑149 (MQ=255)
                                           aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACg                                                                                                        <  1:1436027/149‑1 (MQ=255)
                                            aTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                       >  1:1177075/1‑149 (MQ=255)
                                              aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACACCCGATGCCTGATGCGCGGCTGACGTGa                                                                                                     >  2:1803546/1‑149 (MQ=255)
                                                   ccGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGTGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGGTGCCTGATGCGCCGCTGACGTGACTTAt                                                                                                <  1:1180743/148‑1 (MQ=255)
                                                   aCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTAt                                                                                                <  1:1663531/149‑1 (MQ=255)
                                                    cGCCCGGTGGTTGTATGCGACAGCCGATGCCTGATGCGTCGCTGGCGCGAGTTATCATGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCCTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATc                                                                                               <  2:1201214/149‑1 (MQ=17)
                                                     gtccGTTGACTGTATTCGACAGCCGAGCCCTGGTGCGTCGCTGCCGCGACTTATCAGGCCTACGAGGTGCACTAAACTGGAGGTCGGATAATACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                              <  2:554723/147‑1 (MQ=17)
                                                     gctcGTTGATTGCATTCTACAGCCGAGGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGCCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                              <  2:554724/146‑1 (MQ=17)
                                                       ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTACGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCGTGATGCGCCGCTGACGTGGCTTATCAgg                                                                                            >  1:1785970/1‑149 (MQ=17)
                                                             aTTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTGTCAGGCCTACg                                                                                      >  1:956171/1‑149 (MQ=255)
                                                              ttGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCGTACGACAACCGATGCCTGCTGCGCCGCTGCCGTGACTTGTCAGGCCAACgg                                                                                     >  1:1594917/1‑149 (MQ=17)
                                                              ttGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                     >  2:1403477/1‑149 (MQ=255)
                                                               tggaTTCGACAGCCGATGCCTGATGCGTCGCTTACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                    <  1:283912/146‑1 (MQ=255)
                                                               tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                    <  1:1440754/149‑1 (MQ=255)
                                                               tGTATTCGACAGCCGAGGCCGGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                    <  2:479238/149‑1 (MQ=255)
                                                                ggtttCCACAGCCGATGCCTGAGGCGTCGCTCACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACgggg                                                                                   <  2:1277171/146‑1 (MQ=255)
                                                                 gaTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAATCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGt                                                                                  <  1:484708/148‑1 (MQ=255)
                                                                     cGACAGCCGTTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAc                                                                              <  1:891117/149‑1 (MQ=255)
                                                                             gATGCCGTATGCGTCGCTGACGCGACTTATCAGGCCTGCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATATCGCGTTCGGGATGAAATGAACTGt                                                                      <  2:400094/149‑1 (MQ=14)
                                                                               tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                    <  2:1339486/149‑1 (MQ=255)
                                                                                 ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                  >  1:889483/1‑149 (MQ=255)
                                                                                   tGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGGCATGAACTGTAGGTCg                                                                >  2:1261281/1‑149 (MQ=255)
                                                                                   tGATGCGTCGCTGACGCGACTTAGCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                <  1:682805/149‑1 (MQ=255)
                                                                                   tGATGCGTCGCTGACGCGACTTAACAGGCCTACGAGGTCCACTGAAGTGTAGGTCGGATAAGACGAATGGCGGTGCATCCGACAACCGATGACTGATGCGCCGCTGACGTGGCTTATAAGCCAAAGGTGTTGCACTGAGCTATGGGtgg                                                                >  2:1355568/1‑147 (MQ=255)
                                                                                   agATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGCCGGATAAGCCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                <  1:705460/148‑1 (MQ=255)
                                                                                         gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACCAAGCAGGCCCACGGGGGGCACGGAACCGTAGGGGGGGTAAg                                                          >  1:126270/1‑149 (MQ=14)
                                                                                          tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGa                                                         >  1:1730241/1‑149 (MQ=255)
                                                                                          tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                         >  1:303590/1‑149 (MQ=38)
                                                                                            gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACg                                                       <  2:537838/149‑1 (MQ=38)
                                                                                             cTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGc                                                      >  2:1622349/1‑149 (MQ=17)
                                                                                                aCGCGACTGATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                   <  2:1812864/149‑1 (MQ=255)
                                                                                                    gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCATGAACTGTCGGTCGGATAGAACGCATAGCGt                                               >  2:1133595/1‑149 (MQ=11)
                                                                                                    gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                               >  1:1243637/1‑149 (MQ=9)
                                                                                                     tctgaTCAGGCCTACGAGGGGCACTGAACGGTATGTCGGATAAGACGGATGGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTTCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                              <  2:1463787/145‑1 (MQ=14)
                                                                                                     aCTTATCGGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                              <  2:1665353/149‑1 (MQ=255)
                                                                                                     aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCAACGGGGGGCACTGAACGGTAGGTCGGATAAGACGCATAGCGTc                                              >  2:474770/1‑149 (MQ=255)
                                                                                                       ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATCGCGTccc                                            >  2:128174/1‑147 (MQ=255)
                                                                                                           cAGGCCTACGAGTTGCACTGGACCGTAGGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                        <  1:490179/149‑1 (MQ=255)
                                                                                                            aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAGTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTTGGATAAGACGCATAGCGTCGCATCCg                                       >  1:10695/1‑149 (MQ=255)
                                                                                                            aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                       >  2:850035/1‑149 (MQ=255)
                                                                                                            aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTAGCAGGCCTACGGGGTGCACCGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                       >  1:70706/1‑149 (MQ=255)
                                                                                                             ggCCTACGAGGTGCCCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                      <  2:942326/149‑1 (MQ=255)
                                                                                                                 tACGAGGTGCACTTAACTGGCGGTCGAGTAAGCCGCTGTGGGTCGCCTCCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                  <  1:755505/149‑1 (MQ=11)
                                                                                                                 tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTGAACTGTAGGCCGGCTAACCAGAATAGCGTCGCATCCGACAAc                                  >  1:328396/1‑149 (MQ=14)
                                                                                                                    gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATACGACGCATAGCGTCGCATCCGGCAAccca                               >  2:1418001/1‑147 (MQ=255)
                                                                                                                    gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGa                               <  1:963632/149‑1 (MQ=40)
                                                                                                                     gggggCACTGACCTGGAGGTCGGATAAGACGGATGGCGTCGGATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                              <  2:66709/145‑1 (MQ=255)
                                                                                                                       gTGCACTGAACTGTAGGGCGGCTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                            <  1:459034/149‑1 (MQ=255)
                                                                                                                         gCACTGACCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                          <  1:974344/149‑1 (MQ=38)
                                                                                                                              gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgc                     <  2:1714607/149‑1 (MQ=38)
                                                                                                                                gcTGGAGGTCGGATAAGACGCATGGCGTCGTATTCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                   <  2:1573239/148‑1 (MQ=255)
                                                                                                                                   gTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATACGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                >  1:1669712/1‑149 (MQ=38)
                                                                                                                                   gTAGGTCGGATAAGACGGATGGCAACGGATGTGTCAACCGATGCCTGATGCGCCGCGGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAAACGCATAGCGTCGCATCCGACAACCGATGCCTGATGgcccgc                >  2:1749941/1‑149 (MQ=11)
                                                                                                                                      ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa             <  1:459851/149‑1 (MQ=37)
                                                                                                                                       gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGTACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc            <  2:755565/149‑1 (MQ=38)
                                                                                                                                        tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCAATGCCTGATGCCCCGCGGAcg           >  2:680347/1‑149 (MQ=255)
                                                                                                                                           gATACGACGCATGGGGTGGCATCCGACAACCGATGCCGGATGCGCCGCTGACGTGACATATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGTATAGCGTTGCATCCGACAACCGCTGCCTGATGCGCCGCTGACGCGa        <  2:1355298/149‑1 (MQ=11)
                                                                                                                                            aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATGCGCCGCTGACGCGAc       <  1:1001703/149‑1 (MQ=37)
                                                                                                                                            aTAAGACGGATGGCGTCGCATCCGACAACCGATGCATGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc       <  1:1691477/149‑1 (MQ=37)
                                                                                                                                              aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt     <  1:1204792/149‑1 (MQ=38)
                                                                                                                                                gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt   <  2:889483/149‑1 (MQ=37)
                                                                                                                                                gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt   <  1:1610724/149‑1 (MQ=37)
                                                                                                                                                gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCCCCGCCGACACGACCTAt   >  2:1437005/1‑149 (MQ=255)
                                                                                                                                                gACGGATGGCGCCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACACCCGCTGCCTGATGCGCCGCTGCCGCGACtcac   >  2:249869/1‑146 (MQ=11)
                                                                                                                                                 aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCCATGCCTGACGCGCCGCTGACACGACttcta  >  2:1159436/1‑146 (MQ=255)
                                                                                                                                                 aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCACAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc  >  1:415212/1‑149 (MQ=38)
                                                                                                                                                    |                                                                                                                                                 
TACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATC  >  NC_000913/349525‑349818

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: