Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,070 G→A 5.0% intergenic (‑67/+47) rna ← / ← citT ribonuclease I/citrate/succinate antiporter; citrate carrier

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,0700GA 5.0% 386.9 / 5.6 121intergenic (‑67/+47)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base G (49/64);  new base A (3/3);  total (53/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 125 total aligned reads displayed.

TGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCAC  >  NC_000913/644924‑645218
                                                                                                                                                  |                                                                                                                                                    
tGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCt                                                                                                                                                   >  2:691518/1‑150 (MQ=255)
  cAACGCTAAGGCGTTGGCAGAAGATAAGGGAAGCAGAGAAACCGAGAGCAACGAGGCGTTACGCAAGAATGCTTTCATAATGGGTGTGGAACTCATAAATACACTGAATACTATCTATAAAATAATAAAAAAGCCCTCAATTTGGGCTTa                                                                                                                                                 >  2:919313/1‑150 (MQ=255)
  cAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTa                                                                                                                                                 >  1:515091/1‑150 (MQ=255)
  cAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTa                                                                                                                                                 >  1:857197/1‑150 (MQ=255)
    aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                               <  1:1846477/150‑1 (MQ=255)
       cTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTgg                                                                                                                                            >  1:63791/1‑150 (MQ=255)
          aGGCGTTGGCAGAAGAGAAGGGAAGCAGTGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGCGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                         <  2:1160408/150‑1 (MQ=255)
          aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                         <  2:394756/150‑1 (MQ=255)
          aGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgg                                                                                                                                         <  1:417135/150‑1 (MQ=255)
          aGGCGGTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcgg                                                                                                                                         <  2:1550733/150‑1 (MQ=255)
           ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcggg                                                                                                                                        <  1:1092445/150‑1 (MQ=255)
           ggCGGTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcggg                                                                                                                                        <  1:354624/150‑1 (MQ=255)
            gCGTTGGCAGAAGAGATGGTTAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTCTCTATTAACTCCTAAAAAAGCCCGCCAGGTGTGCTTACCTggcgggc                                                                                                                                       <  1:1995848/150‑1 (MQ=255)
             cGTTGGCAGAAGAGAAGGGTATCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgggcg                                                                                                                                      <  2:308496/150‑1 (MQ=255)
             cGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgggcg                                                                                                                                      <  2:505943/150‑1 (MQ=255)
              gTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                     <  1:1824991/150‑1 (MQ=255)
              gTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                     <  1:1445165/150‑1 (MQ=255)
               ttGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtg                                                                                                                                    <  2:227257/150‑1 (MQ=255)
               gtGGCAGAAGAGCAGGGTAGCAGTGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATTATGGGTGTGGCACTCATACATACACTGAATACTATCTATTAAATCCTAAAAAAGCCCGGCAGGGGGGCTTACCTGGCGGGCGtg                                                                                                                                    <  2:366166/149‑1 (MQ=255)
                 ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgat                                                                                                                                  >  2:540033/1‑150 (MQ=255)
                  gCAGAAGAGATGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgatg                                                                                                                                 <  2:1597666/150‑1 (MQ=255)
                   cAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGTAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgatga                                                                                                                                <  2:489383/150‑1 (MQ=255)
                    agaagagaagGGAAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGACTGCTTGCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                               <  1:283827/150‑1 (MQ=255)
                      aagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGAttt                                                                                                                             <  1:915028/150‑1 (MQ=255)
                         agaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAtt                                                                                                                          >  2:1965343/1‑150 (MQ=255)
                         agaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAtt                                                                                                                          <  2:223843/150‑1 (MQ=255)
                          gtagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                         <  1:857048/148‑1 (MQ=255)
                          gaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                         <  2:877932/150‑1 (MQ=255)
                             gggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGGGGGCTTACCTGGCGGGCGTGATGATTTATTCAGc                                                                                                                      <  1:940920/150‑1 (MQ=255)
                                   catagaAACCGCGAGCAACGCGGCGTTACGACAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGAGGGATTACCTGGCGGGCGTGAATATTTATTCAGCTTTTgg                                                                                                                >  1:1455689/4‑150 (MQ=255)
                                   cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                                <  2:1215364/150‑1 (MQ=255)
                                   cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                                <  2:255606/150‑1 (MQ=255)
                                   cAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                                <  1:1704002/150‑1 (MQ=255)
                                    agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                               <  1:1766452/150‑1 (MQ=255)
                                    agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGc                                                                                                               >  1:1679447/1‑150 (MQ=255)
                                           ccGCGCGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACGCCTACAGACACTGAAGACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTa                                                                                                        >  1:969526/1‑150 (MQ=255)
                                             gcgAGCAACGCGGCGTTTCGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACAGACACTGAATACTTTCTATTAAATCATAAAAAAGCCCGCCAGGTGCGCTTACCTGGCGGGCTTGATGATTTATTAAGCGTTTGGCGAACGTAtt                                                                                                      >  1:1406088/1‑150 (MQ=255)
                                               gAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTat                                                                                                    >  2:256081/1‑149 (MQ=255)
                                                   aaCGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTcc                                                                                                <  1:472420/150‑1 (MQ=255)
                                                    aCGCGGCGTTACGCCATAATGCTTTCATAATGGGTGTGGAACTCATACATACTCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCca                                                                                               <  1:706618/150‑1 (MQ=255)
                                                     cgcgGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCcac                                                                                              >  1:1855816/1‑150 (MQ=255)
                                                           gTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCg                                                                                        <  2:182618/150‑1 (MQ=255)
                                                             tACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCgag                                                                                      <  2:1529353/150‑1 (MQ=255)
                                                               cGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGaa                                                                                    <  2:906868/150‑1 (MQ=255)
                                                               cGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGaa                                                                                    >  1:1339462/1‑150 (MQ=255)
                                                                  cAGAATGCTTTCCTAAGTGGTGGGGAACTCATACATACAATGAATACTATCTATTAAACCAAAAAAAAGCCCGCCAGGGGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCg                                                                                 <  2:371604/150‑1 (MQ=255)
                                                                     aaTGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgcc                                                                              <  1:237830/150‑1 (MQ=255)
                                                                         cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                          <  1:474325/150‑1 (MQ=255)
                                                                         cTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                          <  2:2655/150‑1 (MQ=255)
                                                                         cTTTCATAATGCGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGAGGGCTTACCTCGCGGGCGTGATGAATTATTCAGCGTTAGGCGAACGTATTAGTTAACCATCGCGAGAATCGgccagcc                                                                          >  2:777441/1‑150 (MQ=255)
                                                                          tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCTGGCGTGATGATTTATTCAGCGTATGCCTAACGTATTAGTTCCACATGGCGAGAATCGgccagcca                                                                         >  1:284088/1‑150 (MQ=255)
                                                                          tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcca                                                                         >  2:1164949/1‑150 (MQ=255)
                                                                           ttCATTATGGGTGTGGAACTCATACATACACTGACTACTATCTCTTAAATCATACAAAAGCCCGCCACGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaa                                                                        <  2:531066/150‑1 (MQ=255)
                                                                           ttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaa                                                                        >  1:198842/1‑150 (MQ=255)
                                                                            tCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaac                                                                       >  2:1868754/1‑150 (MQ=255)
                                                                                aaTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTTGTTCCACATGGCGAGAATCGGCCAGCcaaccaac                                                                   >  2:1106583/1‑150 (MQ=255)
                                                                                 aTGGGTGTGGAACTCATACATACACTGATTACTATCTATTAAATCATAAAAAAGCCCCCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACa                                                                  <  1:1691845/150‑1 (MQ=255)
                                                                                  tGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCcaaccaacnn                                                                 <  2:7625/150‑3 (MQ=255)
                                                                                  tGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAg                                                                 <  1:1882622/150‑1 (MQ=255)
                                                                                       gtgGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACa                                                            <  1:716274/150‑1 (MQ=255)
                                                                                       gtgGAACCCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACa                                                            >  2:1565246/1‑150 (MQ=255)
                                                                                         ggAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATa                                                          >  1:1932740/1‑150 (MQ=255)
                                                                                          gAACTCATACATACACTGAATACTATCTCTTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATag                                                         <  1:1766584/150‑1 (MQ=255)
                                                                                          gAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATag                                                         <  2:1488303/150‑1 (MQ=255)
                                                                                            aCTCATACATACACTGAATACTATCTATTAAATCATAAAAAATCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATagag                                                       <  2:1106517/150‑1 (MQ=255)
                                                                                             cTCATACATACACTGAATACTATCTATTAAATCATAAAAAACCCCGCCAGGCGGGCTTACCTGGCGGGTGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAAGCGGCCAGCCAACCAACAGCAACATAgaga                                                      <  2:298456/150‑1 (MQ=255)
                                                                                               catacataCACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                    >  1:1473637/1‑150 (MQ=255)
                                                                                                 tacatacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACATCAACATAGAGATGTa                                                  >  2:1451544/1‑150 (MQ=255)
                                                                                                   catacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTaaa                                                >  2:1693990/1‑150 (MQ=255)
                                                                                                   catacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTaaa                                                >  1:562651/1‑150 (MQ=255)
                                                                                                    atacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAAt                                               <  1:691518/150‑1 (MQ=255)
                                                                                                        acTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcc                                           >  2:1185402/1‑150 (MQ=255)
                                                                                                           gAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcaccc                                        >  2:196421/1‑149 (MQ=255)
                                                                                                                 tatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGAtt                                  >  1:680811/1‑150 (MQ=255)
                                                                                                                 tatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTg                                 <  2:238201/150‑1 (MQ=255)
                                                                                                                   tctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgc                                <  2:1001248/150‑1 (MQ=255)
                                                                                                                    ttatTAAATCATAAAAAAGACCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGACAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                               <  2:746649/149‑1 (MQ=255)
                                                                                                                    ctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                               <  1:1302463/150‑1 (MQ=255)
                                                                                                                         aaaTGATAAAAATGCCCGACAGGTGGGCTTACGTTGCGTGCGTGATGATTTATTCAGCTTTTGGCGAACGTATTATTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTTCGCCAAg                          <  1:919313/150‑1 (MQ=255)
                                                                                                                            tCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGCTTTATTCAGCGTTTGGCTAACGTATTACTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCGCCCCGAAGATTGCGCCacgacg                       >  1:241552/1‑150 (MQ=255)
                                                                                                                            tCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTAGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACg                       >  2:591915/1‑150 (MQ=255)
                                                                                                                             cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGc                      >  1:1377393/1‑150 (MQ=255)
                                                                                                                              aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCACCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAATACGcc                     >  2:1728745/1‑150 (MQ=255)
                                                                                                                               tAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTTGTTACACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGAGGTAAATCACCCCGAAGATTGCGCCAAGACGCCa                    >  2:1776514/1‑150 (MQ=255)
                                                                                                                                aaaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAg                   >  2:379460/1‑150 (MQ=255)
                                                                                                                                aaaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAg                   >  2:1272753/1‑150 (MQ=255)
                                                                                                                                 aaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGt                  >  2:1829363/1‑150 (MQ=255)
                                                                                                                                    aaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAAt               >  2:1327710/1‑150 (MQ=255)
                                                                                                                                      gCCCGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGGTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCCACCAACCGCTACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCt             >  2:1796190/1‑150 (MQ=255)
                                                                                                                                       cccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt            <  2:761646/150‑1 (MQ=255)
                                                                                                                                       cccGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt            >  1:835500/1‑150 (MQ=255)
                                                                                                                                         cGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtttt          <  2:140464/150‑1 (MQ=255)
                                                                                                                                         cGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtttt          <  1:1892498/150‑1 (MQ=255)
                                                                                                                                          gCCCGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTg         <  2:969526/150‑1 (MQ=255)
                                                                                                                                           ccAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGa        <  1:1581962/150‑1 (MQ=255)
                                                                                                                                               ggtggCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTc    <  2:1702238/147‑1 (MQ=255)
                                                                                                                                               gTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTc    <  2:1385995/150‑1 (MQ=255)
                                                                                                                                                tGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTca   <  1:540033/150‑1 (MQ=255)
                                                                                                                                                gtggCTTACGTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTca   <  1:474256/148‑1 (MQ=255)
                                                                                                                                                 gggCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTcac  <  1:1868754/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCAC  >  NC_000913/644924‑645218

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: