Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,490 C→T 8.1% G132D (GGT→GAT)  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,4900CT 8.1% 231.7 / 10.0 74G132D (GGT→GAT) eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (4/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

GCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTT  >  NC_000913/2566344‑2566633
                                                                                                                                                  |                                                                                                                                               
gcagGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGAcc                                                                                                                                                <  4:1017897/148‑1 (MQ=255)
   ggCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGc                                                                                                                                             >  3:151844/1‑148 (MQ=255)
       accaccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATc                                                                                                                                         >  3:131220/1‑148 (MQ=255)
           ccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  4:151844/148‑1 (MQ=255)
              aCGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACAGCCGATGTCGATATTCAGTACCAGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGAtt                                                                                                                                  >  4:581790/1‑148 (MQ=255)
                gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                <  3:581790/148‑1 (MQ=255)
                                  gCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              <  4:131220/148‑1 (MQ=255)
                                   cTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     >  4:507171/1‑124 (MQ=255)
                                   cTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  3:507171/124‑1 (MQ=255)
                                        ttttCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  4:361546/1‑148 (MQ=255)
                                           tCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        <  2:367104/145‑1 (MQ=255)
                                           tCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  1:367104/1‑145 (MQ=255)
                                                  tCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                                                  <  4:966028/144‑1 (MQ=255)
                                                  tCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                                                  >  3:966028/1‑144 (MQ=255)
                                                     aaCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCg                                                                                                                    >  4:441624/1‑123 (MQ=255)
                                                     aaCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCg                                                                                                                    <  3:441624/123‑1 (MQ=255)
                                                        aGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACgg                                                                                                                                     <  2:824277/103‑1 (MQ=255)
                                                        aGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACgg                                                                                                                                     >  1:824277/1‑103 (MQ=255)
                                                           gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                          >  1:228520/1‑95 (MQ=255)
                                                           gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                          <  2:228520/95‑1 (MQ=255)
                                                           gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACt                                                                                     >  3:1064957/1‑148 (MQ=255)
                                                            cGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    <  4:1064957/148‑1 (MQ=255)
                                                              tAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGAc                                                                                                                                                 <  3:94639/85‑1 (MQ=255)
                                                              tAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGAc                                                                                                                                                 >  4:94639/1‑85 (MQ=255)
                                                                gTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  2:176026/95‑1 (MQ=255)
                                                                gTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     >  1:176026/1‑95 (MQ=255)
                                                                     cgGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgcgc                                                                           >  4:45131/1‑148 (MQ=255)
                                                                            aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAAt                                                                                                 <  3:479205/119‑1 (MQ=255)
                                                                            aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAAt                                                                                                 >  4:479205/1‑119 (MQ=255)
                                                                                      tCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAg                                                          >  1:143683/1‑148 (MQ=255)
                                                                                       cGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                         >  2:519674/1‑148 (MQ=255)
                                                                                         atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgc                                                                             <  1:841069/126‑1 (MQ=255)
                                                                                         atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgc                                                                             >  2:841069/1‑126 (MQ=255)
                                                                                         atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAt                                                       >  1:387403/1‑148 (MQ=255)
                                                                                          tatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAtt                                                      >  4:372405/1‑148 (MQ=255)
                                                                                            ttCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGc                                                    <  1:519674/148‑1 (MQ=255)
                                                                                             tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                <  4:232720/71‑1 (MQ=255)
                                                                                             tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                >  3:232720/1‑71 (MQ=255)
                                                                                                  aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTgg                                                                                                             <  2:679243/85‑1 (MQ=255)
                                                                                                  aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTgg                                                                                                             >  1:679243/1‑85 (MQ=255)
                                                                                                  aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGttt                                              <  2:387403/148‑1 (MQ=255)
                                                                                                   cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGtttt                                             <  3:361546/148‑1 (MQ=255)
                                                                                                   cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGtttt                                             <  2:143683/148‑1 (MQ=255)
                                                                                                    ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCggg                                                              <  1:1054212/130‑1 (MQ=255)
                                                                                                    ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCggg                                                              >  2:1054212/1‑130 (MQ=255)
                                                                                                       acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    >  3:1029730/1‑105 (MQ=255)
                                                                                                       acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    <  4:1029730/105‑1 (MQ=255)
                                                                                                             ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAgg                                                                                                                             >  1:1017598/1‑58 (MQ=255)
                                                                                                             ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAgg                                                                                                                             <  2:1017598/58‑1 (MQ=255)
                                                                                                             ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              <  2:821735/73‑1 (MQ=255)
                                                                                                             ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              >  1:821735/1‑73 (MQ=255)
                                                                                                               gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCTGGCg                                                            <  2:662367/121‑1 (MQ=255)
                                                                                                               gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCg                                                            >  1:662367/1‑121 (MQ=255)
                                                                                                                ttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagag                                                                                <  1:197332/100‑1 (MQ=255)
                                                                                                                ttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagag                                                                                >  2:197332/1‑100 (MQ=255)
                                                                                                                     cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                                           >  4:592104/1‑132 (MQ=255)
                                                                                                                     cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                                           <  3:592104/132‑1 (MQ=255)
                                                                                                                     cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAGCGGtgat                           >  2:77911/1‑148 (MQ=255)
                                                                                                                       gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgccg                                                                >  3:151594/1‑109 (MQ=255)
                                                                                                                       gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgccg                                                                <  4:151594/109‑1 (MQ=255)
                                                                                                                       gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatga                         <  3:372405/148‑1 (MQ=255)
                                                                                                                        agaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTg                                                     <  1:731652/117‑1 (MQ=255)
                                                                                                                        aaaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTg                                                     >  2:731652/1‑119 (MQ=255)
                                                                                                                         aaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        <  2:709790/67‑1 (MQ=255)
                                                                                                                         aaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  1:709790/1‑67 (MQ=255)
                                                                                                                              tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGc                  <  3:45131/148‑1 (MQ=255)
                                                                                                                              tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGCGCGAGCATTGCGGGTTTTCGCGCTTTCa                                   >  2:294420/1‑131 (MQ=255)
                                                                                                                              tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTAGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCa                                   <  1:294420/131‑1 (MQ=255)
                                                                                                                                       cGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATc         >  4:993567/1‑148 (MQ=255)
                                                                                                                                       cGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATc         >  1:49221/1‑148 (MQ=255)
                                                                                                                                         gCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCaa       <  1:77911/148‑1 (MQ=255)
                                                                                                                                           tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgatga                      >  4:522667/1‑131 (MQ=255)
                                                                                                                                           tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgatga                      <  3:522667/131‑1 (MQ=255)
                                                                                                                                           tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACg     >  4:992314/1‑148 (MQ=255)
                                                                                                                                            ccGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGc    >  3:803444/1‑148 (MQ=255)
                                                                                                                                            ccGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGc    <  3:993567/148‑1 (MQ=255)
                                                                                                                                              gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCtt  >  4:277872/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTT  >  NC_000913/2566344‑2566633

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: