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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,566,2660CA 6.4% 425.7 / 11.0 125G207C (GGC→TGC) eutAreactivating factor for ethanolamine ammonia‑lyase EutBC
Reads supporting (aligned to +/- strand):  ref base C (63/54);  new base A (3/5);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 128 total aligned reads displayed.

GGCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGG  >  NC_000913/2566118‑2566415
                                                                                                                                                    |                                                                                                                                                     
ggCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACATTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                      <  1:1053786/150‑1 (MQ=255)
   gTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCa                                                                                                                                                   <  2:1301687/150‑1 (MQ=255)
    tAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCAGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCAc                                                                                                                                                  >  1:238893/1‑150 (MQ=255)
       cacCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                               <  1:2312625/150‑1 (MQ=255)
        acCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATGAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCTTCGGTGCAGGCACCGac                                                                                                                                              <  2:322058/150‑2 (MQ=255)
         ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGacg                                                                                                                                             <  2:1981772/150‑3 (MQ=255)
         ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTGCTGCCATCCGCCGGGGCACCTGCACCAGCTGCGCGCCGGGCAGCGAACGGGCGTCGGTGCAGGCACCGacg                                                                                                                                             <  2:2424980/150‑3 (MQ=255)
          cTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGc                                                                                                                                            >  2:1907857/1‑150 (MQ=255)
           tGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCa                                                                                                                                           <  2:754930/150‑1 (MQ=255)
           tGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGACGcc                                                                                                                                           <  1:847597/150‑2 (MQ=255)
            cattgatgcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCAc                                                                                                                                          >  1:118279/8‑150 (MQ=255)
                gcagcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCAt                                                                                                                                      >  1:520339/1‑150 (MQ=255)
                    caAACCGGTTTGCATCAATGCTTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                  >  2:1533513/1‑150 (MQ=255)
                    caAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                  <  1:1625187/150‑1 (MQ=255)
                    caAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                  <  2:837310/150‑1 (MQ=255)
                    caAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                  >  2:1351615/1‑150 (MQ=255)
                      aaCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACaa                                                                                                                                <  2:1116725/150‑1 (MQ=255)
                        ccGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatc                                                                                                                              >  1:210677/1‑150 (MQ=255)
                              tGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGc                                                                                                                        >  1:972820/1‑150 (MQ=255)
                               gCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGcc                                                                                                                       <  1:712115/150‑1 (MQ=255)
                                cATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGccc                                                                                                                      >  2:721858/1‑150 (MQ=255)
                                    aaTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                  >  1:33130/1‑150 (MQ=255)
                                    aaTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCAGGt                                                                                                                  <  1:976809/150‑1 (MQ=255)
                                         cTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATg                                                                                                             <  2:1772764/150‑1 (MQ=255)
                                          tGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGa                                                                                                            >  2:588224/1‑150 (MQ=255)
                                          tGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGa                                                                                                            >  2:616274/1‑150 (MQ=255)
                                           gcgcGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAg                                                                                                           <  1:1907857/150‑1 (MQ=255)
                                              cgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGt                                                                                                        <  2:2224530/150‑1 (MQ=255)
                                              cgAGCGGTGAGAGCGTTCCGTCAATCAATTCGAAAATCAGTTCTGCCATCCGCAGGGTCACCTGCAACAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGCAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGt                                                                                                        >  2:1897922/1‑150 (MQ=255)
                                                   ggTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAAcc                                                                                                   >  2:1610693/1‑150 (MQ=255)
                                                       agagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcgc                                                                                               <  2:1373216/150‑1 (MQ=255)
                                                        gagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcgcg                                                                                              <  1:275358/150‑1 (MQ=255)
                                                        gagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTACGCTCATAAAccggggc                                                                                              <  1:2007951/150‑6 (MQ=255)
                                                         agCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACgcgc                                                                                             >  1:1777231/1‑150 (MQ=255)
                                                                cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                      <  2:1820559/150‑1 (MQ=255)
                                                                cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                      >  1:1225437/1‑150 (MQ=255)
                                                                  tCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTg                                                                                    <  1:1351615/150‑1 (MQ=255)
                                                                     aTCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCg                                                                                 <  2:409561/150‑1 (MQ=255)
                                                                      tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCgg                                                                                <  1:1571524/150‑1 (MQ=255)
                                                                         cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGACCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                             >  2:1312868/1‑150 (MQ=255)
                                                                         cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                             >  1:286149/1‑150 (MQ=255)
                                                                         cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                             >  1:2185276/1‑150 (MQ=255)
                                                                           tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                           >  2:1289993/1‑150 (MQ=255)
                                                                           tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                           <  2:1128376/150‑1 (MQ=255)
                                                                           tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                           >  2:350047/1‑150 (MQ=255)
                                                                            cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCca                                                                          <  1:2415996/150‑1 (MQ=255)
                                                                             gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcag                                                                         <  1:848558/150‑1 (MQ=255)
                                                                             caCAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcag                                                                         <  1:1263776/149‑1 (MQ=255)
                                                                              aCAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagc                                                                        <  2:972820/150‑1 (MQ=255)
                                                                              aCAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagc                                                                        <  2:448694/150‑1 (MQ=255)
                                                                               cAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagca                                                                       >  1:1814188/1‑150 (MQ=255)
                                                                               cAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagca                                                                       >  2:1455348/1‑150 (MQ=255)
                                                                               cAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagca                                                                       >  1:66695/1‑150 (MQ=255)
                                                                                aaTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagcag                                                                      >  1:810611/1‑150 (MQ=255)
                                                                                aaTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagcag                                                                      >  1:1079322/1‑150 (MQ=255)
                                                                                 aTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTGTCCAGCAgg                                                                     >  2:715422/1‑150 (MQ=255)
                                                                                  tCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGc                                                                    <  1:426116/150‑1 (MQ=255)
                                                                                     gTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGac                                                                 <  2:1929914/150‑1 (MQ=255)
                                                                                        cTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGaccac                                                              >  2:1192532/1‑150 (MQ=255)
                                                                                          gCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCCCCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCg                                                            >  1:1425515/1‑150 (MQ=255)
                                                                                              tCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGt                                                        <  1:1192532/150‑1 (MQ=255)
                                                                                               ccgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGtt                                                       <  2:1996510/150‑1 (MQ=255)
                                                                                                cgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGttt                                                      >  2:1450429/1‑149 (MQ=255)
                                                                                                cgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTg                                                      <  2:1310692/150‑1 (MQ=255)
                                                                                                  ccgTGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACTGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGag                                                    <  2:1299/150‑1 (MQ=255)
                                                                                                   cgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGagg                                                   <  1:404035/150‑1 (MQ=255)
                                                                                                      gTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcag                                                >  2:2080619/1‑150 (MQ=255)
                                                                                                        cACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTc                                                                            <  2:319876/120‑1 (MQ=255)
                                                                                                        cACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcaggc                                              <  1:616415/150‑1 (MQ=255)
                                                                                                         aCCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGAACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGAAGTTGAggcaggca                                             >  2:1652738/1‑150 (MQ=255)
                                                                                                          ccTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                            >  1:770646/1‑150 (MQ=255)
                                                                                                          ccTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                            <  2:1869695/150‑1 (MQ=255)
                                                                                                             gCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGCCCACCGACGTTGAGGCAGGCAGTGc                                         >  1:1910203/1‑150 (MQ=255)
                                                                                                             gCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGc                                         >  1:1096818/1‑150 (MQ=255)
                                                                                                             gCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGc                                         <  1:1925085/150‑1 (MQ=255)
                                                                                                                ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGc                                      >  1:1139103/1‑150 (MQ=255)
                                                                                                                  aGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                                    <  1:974523/150‑1 (MQ=255)
                                                                                                                   gCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGa                                   >  2:1426544/1‑150 (MQ=255)
                                                                                                                   gCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGa                                   <  1:291668/150‑1 (MQ=255)
                                                                                                                     tGCGCGCCGGTAAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAtt                                 <  1:1315344/150‑1 (MQ=255)
                                                                                                                          gcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTcc                            <  1:635372/150‑1 (MQ=255)
                                                                                                                           ccGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCg                           <  2:810142/150‑1 (MQ=255)
                                                                                                                            cGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCgg                          >  2:613239/1‑150 (MQ=255)
                                                                                                                                 aGCGAACGGGCGGCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCg                     <  2:2436886/150‑1 (MQ=255)
                                                                                                                                    gAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAAc                  >  1:110490/1‑150 (MQ=255)
                                                                                                                                      aCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAg                <  1:1906321/150‑1 (MQ=255)
                                                                                                                                       cGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg               >  2:485905/1‑150 (MQ=255)
                                                                                                                                       agtgcttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg               >  1:1695780/7‑150 (MQ=255)
                                                                                                                                        gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAggg              <  2:1366874/150‑1 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGc             >  2:2062763/1‑150 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGc             >  1:543439/1‑150 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGc             >  1:729724/1‑150 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCAGGAATCGAACAGGGc             >  2:1597695/1‑150 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCAGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGc             >  2:866714/1‑150 (MQ=255)
                                                                                                                                         ggCGTCGGTGCCTGCACCGAAGCACTCAACCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGCCGACCACCGACGTTGAGGCAGGCAGTGCCGCAGAATTTTCCGGCATCGAACACGGc             >  2:476341/1‑150 (MQ=255)
                                                                                                                                           cttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAt                                  >  2:2370152/3‑127 (MQ=255)
                                                                                                                                           cttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAt                                  <  1:2370152/125‑1 (MQ=255)
                                                                                                                                                   ccTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTcgcg   >  2:178792/1‑150 (MQ=255)
                                                                                                                                                    cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  <  1:1533513/150‑1 (MQ=255)
                                                                                                                                                    cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  >  1:24758/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGG  >  NC_000913/2566118‑2566415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: