Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,073 T→C 5.9% intergenic (‑70/+44) 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,0730TC 5.9% 237.9 / 4.8 86intergenic (‑70/+44)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base T (29/50);  new base C (3/2);  total (33/53)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

AACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCAC  >  NC_000913/644927‑645218
                                                                                                                                                  |                                                                                                                                                 
aaCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGACTACTAGCTATTAAATCATAAAAAAGCCCGCCAGGTGTGCTTa                                                                                                                                                 >  2:1137207/1‑149 (MQ=255)
   gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                              >  1:519981/1‑149 (MQ=255)
   gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGGGGTCTTACCt                                                                                                                                              >  1:362350/1‑149 (MQ=255)
        ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                         <  2:735255/149‑1 (MQ=255)
        ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                         >  1:452877/1‑149 (MQ=255)
        ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGATTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                         <  1:268960/149‑1 (MQ=255)
        ggCGTAGGCAGACGAGCAGGGACGCAGAGCAACCGCGCGCAACGCGGCGTTACGCCAGAATGCGTTCATAATGGGGGTGGATCTCATACATACACTGAATACTATCTATTAACTCATACATAAGCCCGCCAGGTAAGCCCACCTggcgg                                                                                                                                         <  2:703783/149‑1 (MQ=255)
         gCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTAAGCCCACCTggcggg                                                                                                                                        <  1:1095031/149‑1 (MQ=255)
          cGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgggc                                                                                                                                       <  1:298146/149‑1 (MQ=255)
            ttGGCAGACCTGAATCGAAGCAGACACACCGGGAGGATCTCGGCGTTACGCCTGAATGCTTTCATAATCGGATTGGAACTCATACATACACTGAATAATATCTATTAAATCATAAAACAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                     <  2:1154450/149‑1 (MQ=255)
            ttGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGt                                                                                                                                     <  1:199501/149‑1 (MQ=255)
              ggCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtga                                                                                                                                   >  2:180910/1‑149 (MQ=255)
                  gaagagcagGGAAGCCGAGAACCCGCGAGCAACGCGGCGTTACGCCAGTATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATCCTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                               <  2:1136742/149‑1 (MQ=255)
                  gaagagaagGGAAGCAGAGAACCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTgatgat                                                                                                                               <  1:437943/149‑1 (MQ=255)
                   aagagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGattatt                                                                                                                              >  1:591204/1‑149 (MQ=255)
                    agagaagTGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACAATCTATTAAATAATAAAAAAGCCCCCCACGTGGGCTTACCTGGCGGGCGTGATGAttt                                                                                                                             >  2:269741/1‑149 (MQ=255)
                     gagaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTa                                                                                                                            >  1:228295/1‑149 (MQ=255)
                       taagGGAATCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAtt                                                                                                                          <  1:954941/148‑1 (MQ=255)
                       gaagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTAtt                                                                                                                          <  1:431326/149‑1 (MQ=255)
                        aagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                         <  2:1164543/149‑1 (MQ=255)
                        aagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                         <  2:828441/149‑1 (MQ=255)
                         agGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCa                                                                                                                        <  1:1062536/149‑1 (MQ=255)
                           ggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGc                                                                                                                      <  1:755244/149‑1 (MQ=255)
                               gCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGttt                                                                                                                  <  1:1171717/149‑1 (MQ=255)
                                cAGAGAAACCGCGCGCAACGCGGCGTTACGCCAGAAGGCTGTCATAATGGGTGTGGATCTCATACATAAACTGATTACTATCTATTAAAGCATAAAAAAGCCCGCCAGGGGGGCTTACCGGGCGGGCTTGATGCTTTATTCTGCGTTTg                                                                                                                 >  2:927719/1‑149 (MQ=255)
                                 agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                                <  2:54191/149‑1 (MQ=255)
                                    gaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGa                                                                                                             <  1:200777/149‑1 (MQ=255)
                                     aaaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGaa                                                                                                            <  1:942038/149‑1 (MQ=255)
                                      aaCCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAAc                                                                                                           <  1:421099/149‑1 (MQ=255)
                                        ccGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGt                                                                                                         <  1:1058990/149‑1 (MQ=255)
                                         cgcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTa                                                                                                        >  2:429917/1‑149 (MQ=255)
                                          gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                       <  2:46707/149‑1 (MQ=255)
                                          gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                       <  2:715815/149‑1 (MQ=255)
                                          gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                       <  1:500851/149‑1 (MQ=255)
                                          gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                       <  1:76626/149‑1 (MQ=255)
                                          gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAAAACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                       <  1:505306/149‑1 (MQ=255)
                                            gAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTa                                                                                                     <  1:385461/149‑1 (MQ=255)
                                               cAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTTAATCATAAAAACGCCAGCGAGGTGGGCTTACCTGGCGGGCATGATGATTTATTCAGCGTTTGGCGAACGTATTAGtt                                                                                                  >  1:221898/1‑149 (MQ=255)
                                                aaCGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTATTAAAGCGTTTGGCGAACGTATTAGTTc                                                                                                 >  1:844915/1‑149 (MQ=255)
                                                 aCGCGGCGTTACGCCAGAATGCTTTCATAATGGGAGTGCAACTCATACATACACTGAATACAATCTATTAAATCATAAAAAATCCCGCCAGGTAGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTATGGCGAACGTATTAGTTcc                                                                                                >  2:924388/1‑149 (MQ=255)
                                                  cgtggggTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCca                                                                                               <  2:85730/143‑1 (MQ=255)
                                                     ggcgTTACGCCAGAATGCTTTCATAATGGGTGTTGAACTCATACATACACTGAATACTAGCTATTAAATCATAAAAAAGACCGCCATGTGGGCTTACCTTGCGGGCGTGATTATTTATTCAGCGTTTGGCTAAAGTATTAGTTCCACAt                                                                                            <  2:20508/149‑1 (MQ=255)
                                                       ggTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACCTACACTGAATACTATCTATTAAATCATAACCCAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATgg                                                                                          <  1:130470/148‑1 (MQ=255)
                                                        gTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGc                                                                                         >  1:30633/1‑149 (MQ=255)
                                                                  aaTGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgc                                                                               <  1:227118/149‑1 (MQ=255)
                                                                     gcgatCATAATTCGTGTTGAACTCATACATACCCTTTATACTATCTTTTAATTCAGAAAAAAGCCCCGCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGATTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccag                                                                            <  2:1119880/145‑1 (MQ=255)
                                                                       tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                          >  2:721135/1‑149 (MQ=255)
                                                                       tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcc                                                                          >  1:484295/1‑149 (MQ=255)
                                                                         tCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaa                                                                        <  1:186493/149‑1 (MQ=255)
                                                                          cATTATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaac                                                                       <  1:568109/149‑1 (MQ=255)
                                                                          cATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaac                                                                       <  2:452877/149‑1 (MQ=255)
                                                                           aTAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaacc                                                                      <  1:468285/149‑1 (MQ=255)
                                                                             aaTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTCGCGGGCGTGATGATTTATTCAGCGTTTGGCGACCGTATTAGTTCCACATGGCGAGAATCGGCCAGccaaccaa                                                                    >  2:634104/1‑149 (MQ=255)
                                                                             aaTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATAATAAAAAAGCAAGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGccaaccaa                                                                    <  2:709383/149‑1 (MQ=255)
                                                                              aTGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCcaaccaac                                                                   <  1:429917/149‑1 (MQ=255)
                                                                                        aaCTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATag                                                         <  2:1095210/149‑1 (MQ=255)
                                                                                        aaCTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGGGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATag                                                         <  2:789818/149‑1 (MQ=255)
                                                                                             atacatacACTGAATACTATCTATTAAATCATAATAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGAGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCACCCAGCAGCAACATTGAGATg                                                    >  2:237716/1‑149 (MQ=255)
                                                                                             atacatacACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGAATTAGTTCCACATGGCCAGCATCGGCCAGCCAACCAACAGCAACATAGAGATg                                                    >  2:61404/1‑149 (MQ=255)
                                                                                                  tacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGTCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAAt                                               <  2:645602/149‑1 (MQ=255)
                                                                                                  tacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAAt                                               <  2:305508/149‑1 (MQ=255)
                                                                                                      cTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAcc                                           >  1:194453/1‑149 (MQ=255)
                                                                                                       tGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCAccc                                          >  2:372066/1‑149 (MQ=255)
                                                                                                          aTACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGa                                       >  1:959615/1‑149 (MQ=255)
                                                                                                           tACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCATACAAGGCGGGCGTTATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAAACGGCCAGCCAACCAACAGCAACATAGAGATGTAAATAACCCAGaa                                      >  2:216433/1‑149 (MQ=255)
                                                                                                             ctatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGa                                    <  1:339007/149‑1 (MQ=255)
                                                                                                             ctatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGa                                    <  1:221140/149‑1 (MQ=255)
                                                                                                                  tatTAAATGATAAAACAGCCCGCCAGGTGGGCTTACATGGCGAGCGTGATGATTTATTCAGCGTGTGGCGAACTTATTAGATCCAGATGGCGATAATCTGCCAGCCAACGAACAGCAACATAGAGATGTAATTCACCCCGAAGATTgcg                               <  2:437398/149‑1 (MQ=255)
                                                                                                                  tatTAAATCATAACAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                               <  2:143065/149‑1 (MQ=255)
                                                                                                                  tatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                               <  1:337360/149‑1 (MQ=255)
                                                                                                                  tatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcg                               <  2:289500/149‑1 (MQ=255)
                                                                                                                   atTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTgcgc                              >  1:132623/1‑149 (MQ=255)
                                                                                                                    ttAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGcc                             >  2:882347/1‑149 (MQ=255)
                                                                                                                    ttAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGcc                             >  1:1158010/1‑149 (MQ=255)
                                                                                                                       aaTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAg                          >  2:66744/1‑149 (MQ=255)
                                                                                                                         tCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAgtc                        <  1:269741/149‑3 (MQ=255)
                                                                                                                          cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCCGCCCACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGCCg                       >  2:835202/1‑149 (MQ=255)
                                                                                                                            tcatcacgCCCGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGcc                     >  2:430580/8‑149 (MQ=39)
                                                                                                                              aaaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAg                   <  1:97252/149‑1 (MQ=255)
                                                                                                                               tcacgCCCGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGCACGTATTAGTTCCACATGGCGAGAATCGGCCCGCCCACCAACAGCAACAAACAGCAGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGt                  >  2:100268/5‑149 (MQ=255)
                                                                                                                               aaaaGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCACCCAACAGCAACATAGAGATGTAAATCACCCCGAAGAATGCGCCAAGACGCCAGt                  >  2:962185/1‑149 (MQ=255)
                                                                                                                                  tgCCCGCCAGGTGGGCTTACCTTGCGGGCATGATGATTTATTAAGTGTTTGGAGAACGTATTTGTTAGACATGGCGAGAATCTTCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAACATTGCGCCAACACGCCAGTAAt               <  2:221898/148‑1 (MQ=255)
                                                                                                                                     ccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt            >  1:190716/1‑149 (MQ=255)
                                                                                                                                     ccGCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCCCATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGCTTGCGCCACGACGCCAGTAATCtt            >  2:952319/1‑149 (MQ=255)
                                                                                                                                            gTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGAttt     <  2:519981/149‑1 (MQ=255)
                                                                                                                                               ggCGTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACAGGGCGAGAATCCGCCAGCCAACCAACAGCAACCTAGAGATGTAAATCACCCCGACGCTTGCGCCAAGACGCCAGTAATCTTTTGATTTcac  <  2:659914/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTTCAC  >  NC_000913/644927‑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: