Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,529 G→A 7.5% R85R (CGC→CGT citD ← citrate lyase, acyl carrier (gamma) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,5290GA 7.5% 509.9 / 22.3 159R85R (CGC→CGTcitDcitrate lyase, acyl carrier (gamma) subunit
Reads supporting (aligned to +/- strand):  ref base G (69/78);  new base A (9/3);  total (78/81)
Fisher's exact test for biased strand distribution p-value = 7.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 160 total aligned reads displayed.

AGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTC  >  NC_000913/650383‑650670
                                                                                                                                                  |                                                                                                                                             
aGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGcgcg                                                                                                                                             >  1:732012/1‑149 (MQ=255)
 gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGt                                                                                                                                            <  3:803709/149‑1 (MQ=255)
 gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGt                                                                                                                                            <  3:1241116/149‑1 (MQ=255)
     ggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCca                                                                                                                                        <  1:44737/149‑1 (MQ=255)
       tAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                                                                      >  2:1103966/1‑149 (MQ=255)
       tAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                                                                      <  2:842230/149‑1 (MQ=255)
        aGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                                                                                     <  2:298968/149‑1 (MQ=255)
         gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagcag                                                                                                                                    <  3:466239/149‑1 (MQ=255)
          aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                                                                                                                                   <  2:777450/149‑1 (MQ=255)
          aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                                                                                                                                   <  3:286639/149‑1 (MQ=255)
          aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                                                                                                                                   <  4:235617/149‑1 (MQ=255)
            gAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGc                                                                                                                                 >  1:1093867/1‑149 (MQ=255)
              aGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCtt                                                                                                                               <  2:383338/149‑1 (MQ=255)
               ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                                                              <  4:1162/149‑1 (MQ=255)
                 aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                                                            <  3:1207879/149‑1 (MQ=255)
                 aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                                                            <  2:19438/149‑1 (MQ=255)
                 aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                                                            <  2:1487734/149‑1 (MQ=255)
                 aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                                                            <  4:395548/149‑1 (MQ=255)
                  gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAg                                                                                                                           <  1:54390/149‑1 (MQ=255)
                    tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                                                                         <  3:324982/149‑1 (MQ=255)
                    tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                                                                         <  2:509456/149‑1 (MQ=255)
                      cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                                                                                       <  1:1103966/149‑1 (MQ=255)
                       tGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTg                                                                                                                      >  2:58480/1‑149 (MQ=255)
                       tGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTg                                                                                                                      <  2:1131276/149‑1 (MQ=255)
                         aCCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCa                                                                                                                    <  2:29833/149‑1 (MQ=255)
                             tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                                                              <  1:304444/135‑1 (MQ=255)
                                  gCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATg                                                                                                           <  4:1472918/149‑1 (MQ=255)
                                   cATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                                                                          >  2:1568213/1‑149 (MQ=255)
                                    aTTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCa                                                                                                         <  2:411850/149‑1 (MQ=255)
                                         cACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCa                                                                                                    <  2:155794/149‑1 (MQ=255)
                                               gCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgcgc                                                                                              >  1:604251/1‑149 (MQ=255)
                                                 acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                                                                                       <  4:1503708/122‑1 (MQ=255)
                                                   aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                                                                          <  4:474003/149‑1 (MQ=255)
                                                      caacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTc                                                                                       <  3:149031/149‑1 (MQ=255)
                                                        acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                                                                                     >  3:1017267/1‑149 (MQ=255)
                                                         caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                                                                    <  4:202184/149‑1 (MQ=255)
                                                            gctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGAc                                                                                 >  1:571806/1‑149 (MQ=255)
                                                               gcggcgGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAtt                                                                              <  1:1257796/149‑1 (MQ=255)
                                                                  gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg                                                                           <  4:242494/149‑1 (MQ=255)
                                                                            gTTTTACGTTGTCGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCa                                                                                                         >  3:602247/1‑109 (MQ=255)
                                                                            gTTTTACGTTGTCGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCa                                                                                                         <  4:602247/109‑1 (MQ=255)
                                                                               ttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc                                                              >  2:1406858/1‑149 (MQ=255)
                                                                                   gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGAc                                                                                 >  1:1560896/1‑126 (MQ=255)
                                                                                   gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                          >  3:159931/1‑149 (MQ=255)
                                                                                   gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                          >  4:1441838/1‑149 (MQ=255)
                                                                                   gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                          >  2:721807/1‑149 (MQ=255)
                                                                                   gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATACAGCTGTACGCCGCGCACGt                                                          <  1:1477928/149‑1 (MQ=255)
                                                                                    ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCAc                                                            <  3:1570339/146‑1 (MQ=255)
                                                                                    ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGtt                                                         >  2:155620/1‑149 (MQ=255)
                                                                                       ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa                                                      >  2:1467798/1‑149 (MQ=255)
                                                                                       ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa                                                      <  4:159931/149‑1 (MQ=255)
                                                                                         gCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGc                                                    <  3:569691/149‑1 (MQ=255)
                                                                                           aGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAtt                                                                              >  1:848112/1‑121 (MQ=255)
                                                                                            gCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCggg                                                 >  4:129705/1‑149 (MQ=255)
                                                                                            gCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCggg                                                 <  2:955933/149‑1 (MQ=255)
                                                                                                  cGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGaa                                           <  3:754713/149‑1 (MQ=255)
                                                                                                   ggAAAACATTGGCAATCCTCCCATGGCAGAGCCCGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAAc                                          >  1:455821/1‑149 (MQ=255)
                                                                                                    gAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGtt                                                         <  3:1212178/133‑1 (MQ=255)
                                                                                                    gAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACg                                         >  1:768523/1‑149 (MQ=255)
                                                                                                       aTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTcc                                      <  3:341084/149‑1 (MQ=255)
                                                                                                          aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                                                                        <  2:596808/80‑1 (MQ=255)
                                                                                                          aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                                                                        >  1:596808/1‑80 (MQ=255)
                                                                                                           ttGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc                                                                          >  4:587999/1‑109 (MQ=255)
                                                                                                           ttGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc                                                                          <  3:587999/109‑1 (MQ=255)
                                                                                                           ttGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGaa                                  >  2:20397/1‑149 (MQ=255)
                                                                                                           ttGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGaa                                  >  1:1278330/1‑149 (MQ=255)
                                                                                                               cAATCCTCCCATGGCAGAGCCGGGATATCGCTGGCGAGAGCTTGTAGGGCTACCAGTCGTGCACGTAAAATGTAGTCCAGTTAGCCTTTGTCATCGAAATTCAGCTGTACTGCTAGTACGTTGTAGCGGGTGAGAACGTCCATAAttgt                              <  2:658962/149‑1 (MQ=255)
                                                                                                               cAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtggt                              >  2:881459/1‑149 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGTga                                              >  1:307411/1‑132 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGTga                                              <  2:307411/132‑1 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTcc                                      >  2:373607/1‑140 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCa                                     <  4:195962/141‑1 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCa                                     >  3:195962/1‑141 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtggtg                             <  4:90587/149‑1 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtggtg                             <  2:885611/149‑1 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtggtg                             >  2:801128/1‑149 (MQ=255)
                                                                                                                aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtggtg                             <  2:1278330/149‑1 (MQ=255)
                                                                                                                  tcctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATggtggt                           <  1:16370/149‑1 (MQ=255)
                                                                                                                  tcctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATggtggt                           >  1:347886/1‑149 (MQ=255)
                                                                                                                      cccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGa                       >  1:763552/1‑149 (MQ=255)
                                                                                                                       ccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGaa                      >  1:248194/1‑149 (MQ=255)
                                                                                                                         aTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAAtt                    <  1:1467798/149‑1 (MQ=255)
                                                                                                                           ggCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGc                  <  3:211805/149‑1 (MQ=255)
                                                                                                                            gCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCa                 <  3:1025273/149‑1 (MQ=255)
                                                                                                                            gCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCa                 >  3:231742/1‑149 (MQ=255)
                                                                                                                            gCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTACa                 >  1:894003/1‑149 (MQ=255)
                                                                                                                             cAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCAt                <  1:155620/149‑1 (MQ=255)
                                                                                                                              agagCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATc               >  1:274923/1‑149 (MQ=255)
                                                                                                                                agCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                                                                         <  2:535396/89‑1 (MQ=255)
                                                                                                                                 gCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAATTGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAAc                                          <  1:1387783/119‑1 (MQ=255)
                                                                                                                                 gCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAAc                                          >  2:1387783/1‑119 (MQ=255)
                                                                                                                                      gATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACt       <  1:514096/149‑1 (MQ=255)
                                                                                                                                         cctgCTGGCACGCGCCAGCGGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt    <  4:1413663/146‑1 (MQ=255)
                                                                                                                                         cctgCTGGCACGCGCCAGCGGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt    >  3:1413663/4‑149 (MQ=255)
                                                                                                                                         cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt    >  3:511562/4‑149 (MQ=255)
                                                                                                                                         cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt    >  1:587554/4‑149 (MQ=255)
                                                                                                                                          ctgCTGGCACGGGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt   >  3:915849/3‑149 (MQ=255)
                                                                                                                                          ctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGc                         <  4:743353/125‑1 (MQ=255)
                                                                                                                                          ctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt   >  1:972999/3‑149 (MQ=255)
                                                                                                                                           tgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTtc  >  1:516697/2‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
AGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTC  >  NC_000913/650383‑650670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: