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% 330.2 / 7.8 96V154G (GTC→GGC) citXapo‑citrate lyase phosphoribosyl‑dephospho‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (35/56);  new base C (3/2);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGC  >  NC_000913/647430‑647721
                                                                                                                                                 |                                                                                                                                                  
cgcCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACg                                                                                                                                                  <  2:682278/148‑1 (MQ=255)
       ccagcAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTgcgc                                                                                                                                           >  2:235871/1‑148 (MQ=255)
         agcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCtt                                                                                                                                         <  2:187517/148‑1 (MQ=255)
         agcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCtt                                                                                                                                         <  2:256184/148‑1 (MQ=255)
         agcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCtt                                                                                                                                         <  2:798816/148‑1 (MQ=255)
         agcagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCtt                                                                                                                                         <  2:804355/148‑1 (MQ=255)
           cagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGCCGTCTgcgcgcgtg                                                                                                                                      <  3:909513/149‑6 (MQ=255)
           cagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGt                                                                                                                                      <  3:892838/149‑1 (MQ=255)
           cagGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGt                                                                                                                                      >  3:819866/1‑149 (MQ=255)
            agGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGt                                                                                                                                      <  1:703257/148‑1 (MQ=255)
             gggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGtt                                                                                                                                     <  1:235871/148‑1 (MQ=255)
             gggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGtt                                                                                                                                     <  1:925765/148‑1 (MQ=255)
             gggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTc                                                                                                                                    <  3:504186/149‑1 (MQ=255)
              ggCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTc                                                                                                                                    >  1:794743/1‑148 (MQ=255)
                   aCTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCacaa                                                                                                                              <  3:310564/149‑1 (MQ=255)
                    cTCATCGATTAATGACGTCGCGAGCTTTGTGGTTTTGGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCacaac                                                                                                                             <  3:17611/149‑1 (MQ=255)
                       aTCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAgg                                                                                                                          >  3:172218/1‑149 (MQ=255)
                       aTCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAgg                                                                                                                          >  3:894983/1‑149 (MQ=255)
                         cGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGc                                                                                                                         >  1:197613/1‑148 (MQ=255)
                              aaTGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAgcggc                                                                                                                   <  3:744164/149‑1 (MQ=255)
                                   cGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAggc                                                                                                              <  3:175868/149‑1 (MQ=255)
                                   cGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                                                                                               <  2:476220/148‑1 (MQ=255)
                                   cGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                                                                                               >  1:595444/1‑148 (MQ=255)
                                       gCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAggcggc                                                                                                           <  2:250794/148‑1 (MQ=255)
                                        cAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCa                                                                                                          >  1:300092/1‑148 (MQ=255)
                                           gCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGa                                                                                                      <  4:408848/149‑1 (MQ=255)
                                            cTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCAAAACAGGCAGCGGCGAGGCGGCAGTGa                                                                                                      <  2:197613/148‑1 (MQ=255)
                                              ttGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGTGAGGCGGCAGTGAAt                                                                                                    <  1:320621/148‑1 (MQ=255)
                                              ttGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATa                                                                                                   <  4:289145/149‑1 (MQ=255)
                                              ttGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATa                                                                                                   <  3:320674/149‑1 (MQ=255)
                                                 tgGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTc                                                                                                <  3:193424/149‑1 (MQ=255)
                                                   gTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgc                                                                                              >  3:604744/1‑149 (MQ=255)
                                                    ttttAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTcgcg                                                                                             >  3:497116/1‑149 (MQ=255)
                                                       tAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCgcggc                                                                                           <  1:177034/148‑1 (MQ=255)
                                                        aGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCgcggcg                                                                                          <  2:86993/148‑1 (MQ=255)
                                                         gTTGACGTTGCAGGCATCGACATCGTTCATCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGCCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCggg                                                                                        >  3:777813/1‑149 (MQ=255)
                                                          ttGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCggg                                                                                        >  2:18667/1‑148 (MQ=255)
                                                           tGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGga                                                                                       >  1:851890/1‑148 (MQ=255)
                                                            gACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCGGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGggga                                                                                     <  4:773742/149‑1 (MQ=255)
                                                             aCGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGgaga                                                                                     <  2:478597/148‑1 (MQ=255)
                                                              cGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGgagaga                                                                                   <  4:898012/149‑1 (MQ=255)
                                                               gTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGaa                                                                                  >  3:403796/1‑149 (MQ=255)
                                                               gTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCAGGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGgagaga                                                                                   <  1:344306/148‑1 (MQ=255)
                                                                   cAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAAttta                                                                              <  4:194021/149‑2 (MQ=255)
                                                                   cAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTc                                                                              <  3:393635/149‑1 (MQ=255)
                                                                     ggCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAAttt                                                                               <  1:845119/146‑1 (MQ=255)
                                                                     ggCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAAttt                                                                               >  2:845119/1‑146 (MQ=255)
                                                                         tCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCtt                                                                         <  1:649191/148‑1 (MQ=255)
                                                                             cATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCggg                                                                     <  1:816477/148‑1 (MQ=255)
                                                                             cATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCggg                                                                     >  1:37020/1‑148 (MQ=255)
                                                                               tCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGc                                                                            <  2:572244/139‑1 (MQ=255)
                                                                               tCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGc                                                                            >  1:572244/1‑139 (MQ=255)
                                                                               tCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCg                                                                   >  1:250087/1‑148 (MQ=255)
                                                                                cGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGt                                                                  <  2:163551/148‑1 (MQ=255)
                                                                                 gTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCa                                                                >  3:342675/1‑149 (MQ=255)
                                                                                  ttCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCAACGCGCGCAGACGGCTGCGCGTTGTTCGCACAACAGGGAGCGGCGAGGCGGCAGTGAATAGTAGCGGCGGCAGAGAATTTCGCCTTCGGGCGTCAg                                                               <  4:831391/149‑1 (MQ=255)
                                                                                        agTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGTATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATc                                                         <  3:545313/149‑1 (MQ=255)
                                                                                             cTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGatat                                                    >  4:359936/1‑149 (MQ=255)
                                                                                               ccATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATc                                                   >  2:426352/1‑148 (MQ=255)
                                                                                               ccATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATc                                                   >  1:462642/1‑148 (MQ=255)
                                                                                                cATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATccc                                                 >  3:628820/1‑149 (MQ=255)
                                                                                                 aggCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGGAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATccc                                                 <  1:296172/146‑1 (MQ=255)
                                                                                                  tGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCa                                                >  1:624591/1‑148 (MQ=255)
                                                                                                   gCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCAt                                               >  2:70576/1‑148 (MQ=255)
                                                                                                          aGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGccc                                       <  3:156118/149‑1 (MQ=255)
                                                                                                            tAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCga                                     <  3:57661/149‑1 (MQ=255)
                                                                                                              aaTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCgaga                                   <  3:636676/149‑1 (MQ=255)
                                                                                                              aaTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCgaga                                   <  4:497116/149‑1 (MQ=255)
                                                                                                              aaTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCgag                                    <  2:851890/148‑1 (MQ=255)
                                                                                                               aTCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCCTAACCGCCCGagag                                  <  3:74256/149‑1 (MQ=255)
                                                                                                                  ggTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGAt                               <  4:342675/149‑1 (MQ=255)
                                                                                                                   gTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATg                              <  4:894983/149‑1 (MQ=255)
                                                                                                                      aGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACt                           <  3:230204/149‑1 (MQ=255)
                                                                                                                      aGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACt                           <  3:539002/149‑1 (MQ=255)
                                                                                                                       gTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACGACAGGCAGCGTCGAggcgg                                                                                                            <  1:832739/67‑1 (MQ=255)
                                                                                                                       gTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACGACAGGCAGCGGCGAggcgg                                                                                                            >  2:832739/1‑67 (MQ=255)
                                                                                                                         tGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTAt                         <  1:966705/148‑1 (MQ=255)
                                                                                                                            tGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATccc                                                 >  2:634645/1‑121 (MQ=255)
                                                                                                                            tGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATccc                                                 <  1:634645/121‑1 (MQ=255)
                                                                                                                             gggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCg                                                                       <  3:292486/98‑1 (MQ=255)
                                                                                                                             gggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCg                                                                       >  4:292486/1‑98 (MQ=255)
                                                                                                                              ggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCaa                   >  3:567002/1‑149 (MQ=255)
                                                                                                                              ggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCa                    <  2:697207/148‑1 (MQ=255)
                                                                                                                              ggTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCa                    <  2:634900/148‑1 (MQ=255)
                                                                                                                               gTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGAc                            >  4:950452/1‑139 (MQ=255)
                                                                                                                               gTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGAc                            <  3:950452/139‑1 (MQ=255)
                                                                                                                                    tCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGGATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCa              >  2:162088/1‑148 (MQ=255)
                                                                                                                                     aacaaagcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTAAGTTCAAGCTCaa             >  1:44660/7‑148 (MQ=255)
                                                                                                                                      cACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtg           >  3:125771/1‑149 (MQ=255)
                                                                                                                                       aCGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtg           >  2:414546/1‑148 (MQ=255)
                                                                                                                                         aagcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtggt         >  2:414145/3‑148 (MQ=255)
                                                                                                                                          agcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATg                        >  1:884111/2‑132 (MQ=255)
                                                                                                                                          agcgcAGCCGTCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATg                        <  2:884111/131‑1 (MQ=255)
                                                                                                                                           gcgcAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAAtggtgg       <  1:758966/148‑1 (MQ=255)
                                                                                                                                              taGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGa    >  1:456288/2‑148 (MQ=255)
                                                                                                                                                gACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGc  >  1:47168/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGGGCGTCAGGACATCGATATCCCATAACCGCCCGAGAGGATGACTATGTTCAAGCTCAATGGTGGCGAGC  >  NC_000913/647430‑647721

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