Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,266 C→A 8.0% G207C (GGC→TGC)  eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,566,2660CA 8.0% 299.6 / 13.7 88G207C (GGC→TGC) eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base C (37/44);  new base A (3/4);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

GCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCG  >  NC_000913/2566119‑2566414
                                                                                                                                                   |                                                                                                                                                    
gCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  4:603732/149‑1 (MQ=255)
        ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                              <  2:779041/148‑1 (MQ=255)
        ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                              <  2:495365/148‑1 (MQ=255)
        ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGac                                                                                                                                             <  4:526117/149‑2 (MQ=255)
          tGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGacg                                                                                                                                            <  2:451235/148‑3 (MQ=255)
           gCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCa                                                                                                                                          <  4:626498/149‑1 (MQ=255)
            cGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCAc                                                                                                                                         <  4:267798/149‑1 (MQ=255)
                cagcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCAt                                                                                                                                     <  3:318076/149‑1 (MQ=255)
                   caAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCca                                                                                                                                  <  4:812698/149‑1 (MQ=255)
                      aCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACaa                                                                                                                               >  3:715628/1‑149 (MQ=255)
                           tttGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcat                                                                                                                           <  2:594412/148‑1 (MQ=255)
                            ttGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcatc                                                                                                                          >  2:681695/1‑148 (MQ=255)
                            ttGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCt                                                                                                                         <  3:672588/149‑1 (MQ=255)
                                aTCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGcc                                                                                                                      <  1:340832/148‑1 (MQ=255)
                                  cAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGcccc                                                                                                                    >  1:132165/1‑148 (MQ=255)
                                   aaTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCgg                                                                                                                  >  3:431603/1‑149 (MQ=255)
                                    aTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                 >  3:718253/1‑149 (MQ=255)
                                     tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                 >  1:575610/1‑148 (MQ=255)
                                     tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                 <  2:262798/148‑1 (MQ=255)
                                           cgcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGACGcc                                                                                                                                          >  3:258167/1‑116 (MQ=255)
                                           cgcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGACGcc                                                                                                                                          <  4:258167/117‑2 (MQ=255)
                                            gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         <  4:196129/149‑1 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTTACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                <  3:769206/149‑1 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAt                                                                                                             <  3:18645/136‑1 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCa                                                                                                 <  2:872464/148‑1 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCAACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAt                                                                                                             >  4:18645/1‑136 (MQ=255)
                                                      agagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                <  2:376605/148‑1 (MQ=255)
                                                       gagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcg                                                                                               >  1:786192/1‑148 (MQ=255)
                                                            ttCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGccc                                                                                          <  2:427848/148‑1 (MQ=255)
                                                            ttCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGccc                                                                                          >  2:574613/1‑148 (MQ=255)
                                                            ttCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCg                                                                                         >  4:110376/1‑149 (MQ=255)
                                                             tCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGt                                                                                        >  3:133993/1‑149 (MQ=255)
                                                              ccGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGt                                                                                        <  2:630610/148‑1 (MQ=255)
                                                               cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTg                                                                                       <  2:527848/148‑1 (MQ=255)
                                                                gTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                     <  4:829606/149‑1 (MQ=255)
                                                                 tCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                     <  1:5251/148‑1 (MQ=255)
                                                                 tCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                     >  1:10371/1‑148 (MQ=255)
                                                                    aTCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTc                                                                                 <  4:138536/149‑1 (MQ=255)
                                                                     tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTc                                                                                 >  2:604269/1‑148 (MQ=255)
                                                                     tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTc                                                                                 >  2:160060/1‑148 (MQ=255)
                                                                      cACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCgg                                                                               >  3:142900/1‑149 (MQ=255)
                                                                       aCTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAg                                                                                                          >  3:483282/1‑121 (MQ=255)
                                                                       aCTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAg                                                                                                          <  4:483282/121‑1 (MQ=255)
                                                                        cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCa                                                                                                 <  2:909611/129‑1 (MQ=255)
                                                                        cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCa                                                                                                 >  1:909611/1‑129 (MQ=255)
                                                                        cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGtt                                                                             <  3:492009/149‑1 (MQ=255)
                                                                        cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGtt                                                                             <  3:33884/149‑1 (MQ=255)
                                                                         ttCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                            >  3:414035/1‑149 (MQ=255)
                                                                          tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                            >  1:249109/1‑148 (MQ=255)
                                                                          tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                            >  1:522465/1‑148 (MQ=255)
                                                                            gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                          <  2:854260/148‑1 (MQ=255)
                                                                            gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                          <  1:731701/148‑1 (MQ=255)
                                                                            gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCca                                                                         <  3:47625/149‑1 (MQ=255)
                                                                                aTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGATTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagca                                                                      <  2:805384/148‑1 (MQ=255)
                                                                                 tCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagcag                                                                     >  2:772233/1‑148 (MQ=255)
                                                                                  cAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAgg                                                                    >  1:624438/1‑148 (MQ=255)
                                                                                            aTCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGAc                                                         >  3:294487/1‑149 (MQ=255)
                                                                                                gccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTg                                                     >  3:891940/1‑149 (MQ=255)
                                                                                                 ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                    <  4:308324/149‑1 (MQ=255)
                                                                                                      tCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggca                                                >  1:659476/1‑148 (MQ=255)
                                                                                                       cACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcag                                               >  1:578454/1‑148 (MQ=255)
                                                                                                        aCCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                      >  4:440504/1‑108 (MQ=255)
                                                                                                        aCCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                      <  3:440504/108‑1 (MQ=255)
                                                                                                         ccTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAggcaggca                                            <  4:294487/149‑1 (MQ=255)
                                                                                                           tGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                           >  1:593110/1‑148 (MQ=255)
                                                                                                            gCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGt                                          >  1:684913/1‑148 (MQ=255)
                                                                                                            gCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTACCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTg                                         >  3:666927/1‑149 (MQ=255)
                                                                                                             cACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGc                                        >  4:169898/1‑149 (MQ=255)
                                                                                                               ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGcc                                       <  2:735852/148‑1 (MQ=255)
                                                                                                               ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      >  4:223459/1‑149 (MQ=255)
                                                                                                               ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      <  4:305360/149‑1 (MQ=255)
                                                                                                               ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGACTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      >  3:721305/1‑149 (MQ=255)
                                                                                                                cAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      <  2:593110/148‑1 (MQ=255)
                                                                                                                 aGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGc                                     >  2:756131/1‑148 (MQ=255)
                                                                                                                   cTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGa                                  <  4:702135/149‑1 (MQ=255)
                                                                                                                   cTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTACCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGa                                  <  4:666927/149‑1 (MQ=255)
                                                                                                                    tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAt                                 <  4:891833/149‑1 (MQ=255)
                                                                                                                         gcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAttttt                             <  1:604269/148‑1 (MQ=255)
                                                                                                                                aGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATc                     >  3:214408/1‑149 (MQ=255)
                                                                                                                                   gAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGaa                  >  3:309360/1‑149 (MQ=255)
                                                                                                                                     aCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAAc                 <  1:160060/148‑1 (MQ=255)
                                                                                                                                         gcttCGGTGCAGGCACCGACGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAggg             >  1:93762/4‑148 (MQ=255)
                                                                                                                                         gcttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACa                >  4:456655/4‑145 (MQ=255)
                                                                                                                                         gcttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACa                <  3:456655/142‑1 (MQ=255)
                                                                                                                                           gTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg              >  3:860122/1‑145 (MQ=255)
                                                                                                                                           gTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg              <  4:860122/145‑1 (MQ=255)
                                                                                                                                            tCGGTGCCTGCACCGAAGCACTCATCCACACTCAACTGCCCCGGTTTATGAGCGTAATCCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt          <  2:132165/148‑1 (MQ=255)
                                                                                                                                             cGGTGCCTCCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTa         <  2:786192/148‑1 (MQ=255)
                                                                                                                                                   cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTcgcg  <  4:431603/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCG  >  NC_000913/2566119‑2566414

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: ‑