Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 647,575 A→C 5.2% V154G (GTC→GGC)  citX ← apo‑citrate lyase phosphoribosyl‑dephospho‑CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913647,5750AC 5.2% 594.5 / 18.0 168V154G (GTC→GGC) citXapo‑citrate lyase phosphoribosyl‑dephospho‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (77/82);  new base C (5/4);  total (82/86)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 170 total aligned reads displayed.

CGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGCTT  >  NC_000913/647428‑647723
                                                                                                                                                   |                                                                                                                                                    
tgagtacgccctgcTGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGAc                                                                                                                                                     >  1:1207003/8‑149 (MQ=255)
     tagcccagcAGGGCGTACTAATCGATTAATGACGTCGCAAGCTTTGTGGCTTTAGTTGACGTTGAAGGCATCGACAACGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGTTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTg                                                                                                                                                <  2:1058656/148‑1 (MQ=255)
      agcccagcAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTgc                                                                                                                                               >  2:751172/1‑149 (MQ=255)
       gcccagcAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTgcg                                                                                                                                              >  3:1449988/1‑149 (MQ=255)
        cccagcAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTgcgc                                                                                                                                             >  3:898117/1‑149 (MQ=255)
          cagcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTCCACGCGCGCAGACGGCTGCGCttt                                                                                                                                          <  1:1237121/149‑1 (MQ=255)
           agcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCttt                                                                                                                                          <  1:1310306/149‑1 (MQ=255)
              agGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGCCGTCTgcgcgcgtgg                                                                                                                                       <  4:670111/149‑7 (MQ=255)
               gggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTc                                                                                                                                      <  2:1259401/149‑1 (MQ=255)
                ggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCg                                                                                                                                     >  2:1027394/1‑149 (MQ=255)
                  cGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGca                                                                                                                                   <  1:240378/149‑1 (MQ=255)
                    tACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCttt                                                                                                                                          <  2:1105820/140‑1 (MQ=255)
                     aCTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCacaa                                                                                                                                <  2:1286824/149‑1 (MQ=255)
                     aCTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCacaa                                                                                                                                <  2:463900/149‑1 (MQ=255)
                          tCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGc                                                                                                                           >  2:1134209/1‑149 (MQ=255)
                             aTTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAgc                                                                                                                        <  3:436856/149‑1 (MQ=255)
                              ttAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCGCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAgcg                                                                                                                       <  2:1037528/149‑1 (MQ=255)
                                aaTGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAgcggc                                                                                                                     >  3:122377/1‑149 (MQ=255)
                                 aTGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAgcggcg                                                                                                                    >  4:1145938/1‑149 (MQ=255)
                                  tGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGtt                                                                                                                                       <  1:1308125/129‑1 (MQ=255)
                                   gACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAg                                                                                                                  <  1:500848/149‑1 (MQ=255)
                                   gACGTCGCAAGCTTTGTGGTTTTAGTGGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAg                                                                                                                  >  2:379153/1‑149 (MQ=255)
                                    aCGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGa                                                                                                                   >  4:1385209/1‑147 (MQ=255)
                                           aaGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGt                                                                                                          <  4:746092/149‑1 (MQ=255)
                                           aaGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGt                                                                                                          <  2:1067676/149‑1 (MQ=255)
                                                  gtgGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGt                                                                                                   >  3:1121960/1‑149 (MQ=255)
                                                   tgGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTc                                                                                                  <  2:780303/149‑1 (MQ=255)
                                                   tgGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTc                                                                                                  >  3:614928/1‑149 (MQ=255)
                                                     gTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgc                                                                                                <  2:862121/149‑1 (MQ=255)
                                                     gTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgc                                                                                                >  3:1136962/1‑149 (MQ=255)
                                                      ttttAGTTGACTTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgcg                                                                                               >  1:1098536/1‑149 (MQ=255)
                                                      ttttAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgcg                                                                                               >  3:1314406/1‑149 (MQ=255)
                                                       tttAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCgcgg                                                                                              <  2:654320/149‑1 (MQ=255)
                                                          aGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGt                                                                                                          >  1:1189001/1‑134 (MQ=255)
                                                          aGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGcggcgg                                                                                           <  4:122377/149‑1 (MQ=255)
                                                            ttGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGga                                                                                         >  1:1128544/1‑149 (MQ=255)
                                                            ttGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGga                                                                                         >  3:745677/1‑149 (MQ=255)
                                                             tGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCTGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGgag                                                                                        <  1:756041/149‑1 (MQ=255)
                                                             tGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGgag                                                                                        <  4:206319/149‑1 (MQ=255)
                                                                    gCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAAttt                                                                                 <  2:1320852/149‑1 (MQ=255)
                                                                     cAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGGGAGAATTTc                                                                                <  4:42876/149‑1 (MQ=255)
                                                                     cAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTc                                                                                <  4:163961/149‑1 (MQ=255)
                                                                      aGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTca                                                                               >  4:1168566/1‑148 (MQ=255)
                                                                      aGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCg                                                                               >  2:307563/1‑149 (MQ=255)
                                                                          aTCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCtt                                                                           <  3:191593/149‑1 (MQ=255)
                                                                          aTCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCtt                                                                           >  1:394217/1‑149 (MQ=255)
                                                                           tCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTc                                                                          <  1:733992/149‑1 (MQ=255)
                                                                            cGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCg                                                                         <  3:1114477/149‑1 (MQ=255)
                                                                             gACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCgg                                                                        <  4:189205/149‑1 (MQ=255)
                                                                              aCATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTGCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGATAGAATTTCGCCTTCggg                                                                       >  2:900421/1‑149 (MQ=255)
                                                                               cATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGc                                                                      <  2:556604/149‑1 (MQ=255)
                                                                                aTCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCg                                                                     >  4:626402/1‑149 (MQ=255)
                                                                                 tCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGt                                                                    >  3:506265/1‑149 (MQ=255)
                                                                                  cGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTc                                                                   <  2:1434167/149‑1 (MQ=255)
                                                                                    ttCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAg                                                                 >  4:1454270/1‑149 (MQ=255)
                                                                                        gcagTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACa                                                             <  2:1207003/149‑1 (MQ=255)
                                                                                         cagTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACAt                                                            <  3:658351/149‑1 (MQ=255)
                                                                                         cagTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACAt                                                            <  2:1210837/149‑1 (MQ=255)
                                                                                         cagTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACAt                                                            >  2:674729/1‑149 (MQ=255)
                                                                                         cagTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACAGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGCCCGCAGCGAATAGTCGCGGCGGGAGCGAATTTCGCCTTCGGGCGTCAGGACAt                                                            >  2:200844/1‑149 (MQ=255)
                                                                                           gTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCg                                                          <  1:611932/149‑1 (MQ=255)
                                                                                              ccTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGGCAGGACATCGata                                                       >  2:976357/1‑149 (MQ=255)
                                                                                               cTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGatat                                                      >  4:936217/1‑149 (MQ=255)
                                                                                               cTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGatat                                                      <  2:1098536/149‑1 (MQ=255)
                                                                                               cTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGatat                                                      >  2:680211/1‑149 (MQ=255)
                                                                                               cTCCATGCGGTTGAGCAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGACAGGACATAGatat                                                      >  2:1289137/1‑149 (MQ=255)
                                                                                                tCCATGCGGTTGAGTAAATGGGTCAGTTGATTGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCCGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATc                                                     <  4:658517/149‑1 (MQ=255)
                                                                                                tCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATc                                                     >  3:941425/1‑149 (MQ=255)
                                                                                                 ccATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATcc                                                    <  1:575812/149‑1 (MQ=255)
                                                                                                 ccATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATcc                                                    <  2:1242767/149‑1 (MQ=255)
                                                                                                   aTGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGtt                                                                                                                                       >  1:902779/1‑64 (MQ=255)
                                                                                                    tGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCAt                                                 >  3:321597/1‑149 (MQ=255)
                                                                                                      cGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATaa                                               <  4:745677/149‑1 (MQ=255)
                                                                                                       ggTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAAc                                              >  1:34306/1‑149 (MQ=255)
                                                                                                            aGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGccc                                         <  3:189999/149‑1 (MQ=255)
                                                                                                             gTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCg                                        >  2:1461115/1‑149 (MQ=255)
                                                                                                              tAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCga                                       <  3:1007098/149‑1 (MQ=255)
                                                                                                              tAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCga                                       <  4:898117/149‑1 (MQ=255)
                                                                                                                aaTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGCATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCgaga                                     <  4:103044/149‑1 (MQ=255)
                                                                                                                aaTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGAGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCgaga                                     >  1:515559/1‑149 (MQ=255)
                                                                                                                 aTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGagag                                    >  3:579157/1‑149 (MQ=255)
                                                                                                                    ggTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGAt                                 <  1:307563/149‑1 (MQ=255)
                                                                                                                     gTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATg                                >  2:595268/1‑149 (MQ=255)
                                                                                                                         gTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTa                            >  4:553034/1‑149 (MQ=255)
                                                                                                                          ttGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGGGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTAt                           >  1:152603/1‑149 (MQ=255)
                                                                                                                          ttGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTAt                           <  3:318552/149‑1 (MQ=255)
                                                                                                                            gATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGAGAGGAGGCAGTGAATAGTCGCGGCGGGAGAGCATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGt                         >  2:82958/1‑149 (MQ=255)
                                                                                                                                ggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCaa                     <  1:889476/149‑1 (MQ=255)
                                                                                                                                agTTTTGCCACGCGAGCAGACGGCTGCGCTTTGTTCGCGCAACAGGCAGCGGCGAGTCGGCAGTGGATAGTCGAGGCGGCACAGAATTTCTCCTTCGGGCGTCAGGACATCGATGTCCAATACCCGCCCGAGAGGATGACTATATTaaa                     <  4:239891/148‑1 (MQ=255)
                                                                                                                                 gTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAg                    >  1:540702/1‑149 (MQ=255)
                                                                                                                                 gTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAg                    >  2:1156244/1‑149 (MQ=255)
                                                                                                                                 gTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAg                    >  3:1318446/1‑149 (MQ=255)
                                                                                                                                    tttCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTc                 >  4:303935/1‑149 (MQ=255)
                                                                                                                                       ccACGCGCGCTGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTTACTAGTCGCGGCGTGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAt              <  1:1062131/149‑1 (MQ=255)
                                                                                                                                          cgcgcgCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtggt           <  4:614928/149‑1 (MQ=255)
                                                                                                                                          aaagcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtggt           >  4:1260878/4‑149 (MQ=255)
                                                                                                                                           aagcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGt                         <  4:1197017/132‑1 (MQ=255)
                                                                                                                                             gcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCaa                     <  3:1150724/136‑1 (MQ=255)
                                                                                                                                               gcAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGa      <  3:312823/149‑1 (MQ=255)
                                                                                                                                                   aCGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGCtt  <  4:321597/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGCTT  >  NC_000913/647428‑647723

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