Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,266 C→A 7.3% 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 7.3% 622.2 / 28.1 176G207C (GGC→TGC) eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base C (75/87);  new base A (7/6);  total (83/93)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 182 total aligned reads displayed.

GGCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCG  >  NC_000913/2566118‑2566414
                                                                                                                                                    |                                                                                                                                                    
ggCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGcc                                                                                                                                                      <  2:873443/149‑1 (MQ=255)
 gagTAACACCTGCGGGCAGCAATCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  2:593101/147‑1 (MQ=255)
 gCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  2:1436760/149‑1 (MQ=255)
 gCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  4:1875566/149‑1 (MQ=255)
 gCGTAACACCTGCAGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  2:1691868/149‑1 (MQ=255)
      acacCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCAcc                                                                                                                                                >  2:769955/1‑149 (MQ=255)
          cTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAg                                                                                                                                            >  3:1775794/1‑149 (MQ=255)
            gCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCa                                                                                                                                          <  2:1836370/149‑1 (MQ=255)
             cGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCAc                                                                                                                                         <  1:745715/149‑1 (MQ=255)
                    caAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCTGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCca                                                                                                                                  <  4:1775794/149‑1 (MQ=255)
                     aaaCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                 >  2:1196039/1‑149 (MQ=255)
                      aaCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGaa                                                                                                                                             <  2:454964/136‑1 (MQ=255)
                      aaCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGaa                                                                                                                                             >  1:454964/1‑136 (MQ=255)
                       aCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACaa                                                                                                                               <  4:465782/149‑1 (MQ=255)
                        ccGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAat                                                                                                                              >  3:397862/1‑149 (MQ=255)
                          ggTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatca                                                                                                                            >  3:737065/1‑149 (MQ=255)
                           gTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcat                                                                                                                           <  3:1025256/149‑1 (MQ=255)
                            tttGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcatc                                                                                                                          <  2:1242125/149‑1 (MQ=255)
                                 aTCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCt                                                                                                                                                     <  1:1294188/117‑1 (MQ=255)
                                 aTCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGccc                                                                                                                     >  3:892880/1‑149 (MQ=255)
                                 aTCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGac                                                                                                                                             >  1:1002141/1‑124 (MQ=255)
                                  tCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGcccc                                                                                                                    >  1:929960/1‑149 (MQ=255)
                                  tCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGcccc                                                                                                                    >  4:156183/1‑149 (MQ=255)
                                    aaTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCgg                                                                                                                  >  2:1736048/1‑149 (MQ=255)
                                     aTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCa                                                                                                                                      >  2:264214/1‑128 (MQ=255)
                                     aTGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                 >  4:1213527/1‑149 (MQ=255)
                                       gCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGttt                                                                                                               >  1:835008/1‑149 (MQ=255)
                                         cTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAt                                                                                                             >  1:247397/1‑149 (MQ=255)
                                             gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         <  3:559868/149‑1 (MQ=255)
                                             gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         <  3:1409416/149‑1 (MQ=255)
                                             gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         <  4:1077454/149‑1 (MQ=255)
                                               gAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGt                                                                                                       <  3:923777/149‑1 (MQ=255)
                                                     tGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         >  2:1659740/1‑141 (MQ=255)
                                                     tGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                         <  1:1659740/141‑1 (MQ=255)
                                                      gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                >  2:1803963/1‑149 (MQ=255)
                                                                cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                >  3:667781/1‑139 (MQ=255)
                                                                cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                <  4:667781/139‑1 (MQ=255)
                                                                cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                      >  1:1428320/1‑149 (MQ=255)
                                                                 gTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGa                                                                                                           >  3:1163855/1‑127 (MQ=255)
                                                                  tCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCt                                                                                    <  4:1845206/149‑1 (MQ=255)
                                                                      tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcgc                                                                                              >  2:1124028/1‑135 (MQ=255)
                                                                      tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcgc                                                                                              <  1:1124028/135‑1 (MQ=255)
                                                                      tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCg                                                                                <  2:1645305/149‑1 (MQ=255)
                                                                        aCTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                     <  1:987247/142‑1 (MQ=255)
                                                                        aCTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                     >  2:987247/1‑142 (MQ=255)
                                                                        aCTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGt                                                                              <  1:1197422/149‑1 (MQ=255)
                                                                          ttCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGt                                                                                                                 >  3:840925/1‑112 (MQ=255)
                                                                          ttCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                            <  4:737065/149‑1 (MQ=255)
                                                                          ttCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGt                                                                              <  3:1264114/147‑1 (MQ=255)
                                                                            cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCAc                                                                                                                                                 >  4:1360562/1‑78 (MQ=255)
                                                                            cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCAc                                                                                                                                                 <  3:1360562/78‑1 (MQ=255)
                                                                            cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                          <  3:479037/149‑1 (MQ=255)
                                                                            cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                          <  2:247397/149‑1 (MQ=255)
                                                                                  tCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAgg                                                                    <  2:647926/149‑1 (MQ=255)
                                                                                  tCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAgg                                                                    >  3:38366/1‑149 (MQ=255)
                                                                                  tCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAgg                                                                    >  3:1801399/1‑149 (MQ=255)
                                                                                      ttCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGac                                                                >  1:1843599/1‑149 (MQ=255)
                                                                                      ttCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCATGGCTGTCGGTTTCCAGCAGGCGac                                                                >  1:1005898/1‑149 (MQ=255)
                                                                                         tGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGaccac                                                             <  1:665327/149‑1 (MQ=255)
                                                                                           ccATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                 >  2:1599055/1‑79 (MQ=255)
                                                                                           ccATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                 >  2:161583/1‑79 (MQ=255)
                                                                                           ccATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcac                                                                                                                                 <  1:161583/79‑1 (MQ=255)
                                                                                             aTCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGaa                                                                                                                                             >  2:1407269/1‑65 (MQ=255)
                                                                                                  ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                    <  4:1304436/149‑1 (MQ=255)
                                                                                                  ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                    <  2:1843599/149‑1 (MQ=255)
                                                                                                  ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                    <  4:883845/149‑1 (MQ=255)
                                                                                                  ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCATGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                    <  2:1005898/149‑1 (MQ=255)
                                                                                                      gTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggca                                                <  1:1467360/149‑1 (MQ=255)
                                                                                                       tCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcag                                               >  4:27857/1‑149 (MQ=255)
                                                                                                        cACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                          <  4:912631/121‑1 (MQ=255)
                                                                                                           cTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                      >  4:1094203/1‑106 (MQ=255)
                                                                                                           cTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                      <  3:1094203/106‑1 (MQ=255)
                                                                                                           cTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                           >  3:1054651/1‑149 (MQ=255)
                                                                                                                ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      >  3:202638/1‑149 (MQ=255)
                                                                                                                ccAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                      <  1:938907/149‑1 (MQ=255)
                                                                                                                 cAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGc                                     <  4:1054651/149‑1 (MQ=255)
                                                                                                                  aGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCt                                    <  4:323403/149‑1 (MQ=255)
                                                                                                                     tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                <  4:1530925/86‑1 (MQ=255)
                                                                                                                     tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                                >  3:1530925/1‑86 (MQ=255)
                                                                                                                     tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAt                                 <  2:929960/149‑1 (MQ=255)
                                                                                                                          gcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTc                            <  1:1745490/149‑1 (MQ=255)
                                                                                                                              gTCAGCGAACGGGCGTCGGTGCCTGCACCAAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggca                                                <  3:163768/125‑1 (MQ=255)
                                                                                                                                 aGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                                   <  3:13587/135‑1 (MQ=255)
                                                                                                                                 aGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                                   >  4:1484091/1‑135 (MQ=255)
                                                                                                                                  gCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAtt                                <  3:407462/137‑1 (MQ=255)
                                                                                                                                  gCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAtt                                >  4:407462/1‑137 (MQ=255)
                                                                                                                                     aaCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGAc                                                         <  2:523893/109‑1 (MQ=255)
                                                                                                                                     aaCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGAc                                                         >  1:523893/1‑109 (MQ=255)
                                                                                                                                      aCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACa                >  3:1823144/1‑149 (MQ=255)
                                                                                                                                      aCGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCAGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACa                >  1:251562/1‑149 (MQ=255)
                                                                                                                                        gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTg                                         <  4:661667/122‑1 (MQ=255)
                                                                                                                                        gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg              <  1:1803963/149‑1 (MQ=255)
                                                                                                                                        gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg              <  3:27857/149‑1 (MQ=255)
                                                                                                                                           cttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCt                                    <  1:1695663/122‑1 (MQ=255)
                                                                                                                                            ttCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTaaa                                                                                                    <  4:1799579/58‑1 (MQ=255)
                                                                                                                                            ttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt          >  3:839494/2‑149 (MQ=255)
                                                                                                                                            ttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt          >  1:1413586/2‑149 (MQ=255)
                                                                                                                                             tCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTa         <  4:839494/149‑1 (MQ=255)
                                                                                                                                                 tGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                                   >  1:406430/1‑119 (MQ=255)
                                                                                                                                                   ccTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTcgc   <  4:202638/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCG  >  NC_000913/2566118‑2566414

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