Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,069 G→A 6.7% intergenic (‑66/+48) 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,0690GA 6.7% 226.0 / 4.3 87intergenic (‑66/+48)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base G (29/52);  new base A (4/2);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTT  >  NC_000913/644921‑645215
                                                                                                                                                    |                                                                                                                                                  
gCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCCTACATACACTGAAAACTATCTATTAAATCATAAAAAAGCCCGCCAGGTg                                                                                                                                                    >  2:339007/1‑149 (MQ=255)
gCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTg                                                                                                                                                    >  2:221140/1‑149 (MQ=255)
   tGCAACGCTAAGGCGTTGGCAGATGAGAAGGGAACCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                 <  2:1175534/149‑1 (MQ=255)
   tGCAACGCTAAGGCGGTGGCCGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGGTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGc                                                                                                                                                 <  2:560426/149‑1 (MQ=255)
      aaCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGACTACTAGCTATTAAATCATAAAAAAGCCCGCCAGGTGTGCTTa                                                                                                                                              >  2:1137207/1‑149 (MQ=255)
         gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                           >  1:519981/1‑149 (MQ=255)
         gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGGGGTCTTACCt                                                                                                                                           >  1:362350/1‑149 (MQ=255)
              ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                      >  1:452877/1‑149 (MQ=255)
              ggCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTggcgg                                                                                                                                      <  2:735255/149‑1 (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:828441/149‑1 (MQ=255)
                              aagGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTc                                                                                                                      <  2:1164543/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                                                                                                    <  1:76626/149‑1 (MQ=255)
                                                gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                    <  2:715815/149‑1 (MQ=255)
                                                gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                    <  2:46707/149‑1 (MQ=255)
                                                gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTAt                                                                                                    <  1:500851/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:221140/149‑1 (MQ=255)
                                                                                                                   ctatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGa                                 <  1:339007/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)
                                                                                                                                                    |                                                                                                                                                  
GCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTTGATTT  >  NC_000913/644921‑645215

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: