Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 628,160 A→T 5.9% I212N (ATT→AAT)  eutC ← ethanolamine ammonia‑lyase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947628,1600AT 5.9% 214.6 / 2.9 68I212N (ATT→AAT) eutCethanolamine ammonia‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base T (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCTTCATTCTCCAGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGCCGACCTCATCGGC  >  NC_002947/628019‑628305
                                                                                                                                             |                                                                                                                                                 
tCTTCATTCTCCAGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGtt                                                                                                                                            <  5:3574/149‑1 (MQ=255)
         tcCAGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACACCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAgg                                                                                                                                    >  7:59130/1‑148 (MQ=255)
          ccAGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAATTGCAGTTGCGGTAGGCa                                                                                                                                  >  7:292512/1‑149 (MQ=255)
            aGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATc                                                                                                                                >  3:6152/1‑149 (MQ=255)
                tATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGt                                                                                                                             >  3:151424/1‑148 (MQ=255)
                tATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGt                                                                                                                             >  5:83372/1‑148 (MQ=255)
                  tGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTc                                                                                                                            <  3:368698/147‑1 (MQ=255)
                   ggACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAgg                                                                                                                         >  3:274956/1‑149 (MQ=255)
                    gACCTCGGCTACGTCCTTTAGATTCACCCAGCAAAGCTGGCGCCGGCAGGCTACGCGCATCAAGTCCAACACGCGGTGCGACGCCATGCCATAGCTCAGTAACTCCAGGCGGATAATCTAAATGAAGTTGCGGTAGGCATCGTTCAGGc                                                                                                                        >  6:4393/1‑149 (MQ=255)
                      ccTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGcc                                                                                                                       >  5:50860/1‑148 (MQ=255)
                         cGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGAcc                                                                                                                   >  4:24783/1‑149 (MQ=255)
                               cGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTgg                                                                                                               <  3:18103/147‑1 (MQ=255)
                                gTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAATTGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGt                                                                                                             <  8:292512/148‑1 (MQ=255)
                                gTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAAGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGt                                                                                                             <  4:193704/148‑1 (MQ=255)
                                   cTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCg                                                                                                          >  8:186262/1‑148 (MQ=255)
                                    tttAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGc                                                                                                                                          >  7:313535/1‑115 (MQ=255)
                                    tttAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGc                                                                                                                                          <  8:313535/115‑1 (MQ=255)
                                    tttAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGt                                                                                                         >  5:263160/1‑148 (MQ=255)
                                      tAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTaa                                                                                                       <  7:186262/148‑1 (MQ=255)
                                       aGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCg                                                                                                                                         >  6:331436/1‑113 (MQ=255)
                                       aGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCg                                                                                                                                         <  5:331436/113‑1 (MQ=255)
                                          ttCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGa                                                                                                   >  3:177717/1‑148 (MQ=255)
                                             aCCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAAt                                                                                                >  7:110048/1‑148 (MQ=255)
                                             aCCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAAt                                                                                                <  8:110048/148‑1 (MQ=255)
                                             aCCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATa                                                                                               >  4:201885/1‑149 (MQ=255)
                                               ccaggAAAGCTGGCGCCGGCAGAATTCCCGCATCAAGTACAACACGCGGTGCGCCGCCATCCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCATGGGTGAGTAAGTGAAATa                                                                                               >  2:155503/4‑147 (MQ=255)
                                                   gAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGc                                                                                          >  1:112521/1‑148 (MQ=255)
                                                     aaGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCaa                                                                                       <  8:59130/149‑1 (MQ=255)
                                                            cgcCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTc                                                                                <  4:177717/149‑1 (MQ=255)
                                                                ggCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTgg                                                                                                               <  6:358199/114‑1 (MQ=255)
                                                                ggCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTgg                                                                                                               >  5:358199/1‑114 (MQ=255)
                                                                         cgcgCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCa                                                                     <  3:24783/147‑1 (MQ=255)
                                                                              aTCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCgg                                                              >  6:324301/1‑149 (MQ=255)
                                                                              aTCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCACCgg                                                              >  6:262864/1‑149 (MQ=255)
                                                                                   gTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAgg                                                                                     <  1:22694/121‑1 (MQ=255)
                                                                                   gTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAgg                                                                                     >  2:22694/1‑121 (MQ=255)
                                                                                    tACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTc                                                        >  1:398630/1‑149 (MQ=255)
                                                                                      caacaGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAATTGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCgg                                                              <  3:127608/141‑1 (MQ=255)
                                                                                      caacaGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAATTGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCgg                                                              >  4:127608/1‑141 (MQ=255)
                                                                                       aacaGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGcc                                                     >  4:357026/1‑149 (MQ=255)
                                                                                        acaGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGcc                                                     >  6:126442/1‑148 (MQ=255)
                                                                                        acaGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGcc                                                     >  4:211172/1‑148 (MQ=255)
                                                                                          aGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGAt                                                  <  5:262864/149‑1 (MQ=255)
                                                                                            gCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGaaa                                                                                                 <  2:398493/100‑1 (MQ=255)
                                                                                            gCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGaaa                                                                                                 >  1:398493/1‑100 (MQ=255)
                                                                                                  cgccgcCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATc                                           <  6:83372/148‑1 (MQ=255)
                                                                                                  cgccgcCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCa                                          <  3:211172/149‑1 (MQ=255)
                                                                                                  cgccgcCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCa                                          <  5:126442/149‑1 (MQ=255)
                                                                                                             catAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATcc                                 >  8:4869/1‑147 (MQ=255)
                                                                                                             catAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGa                               <  4:187746/149‑1 (MQ=255)
                                                                                                             catAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGa                               <  5:324301/149‑1 (MQ=255)
                                                                                                             catAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAAGCCAAGGCTGTCGGGTGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGa                               <  6:263160/149‑1 (MQ=255)
                                                                                                              atAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCCTCACCGTCATCCgcg                              >  8:387253/1‑147 (MQ=255)
                                                                                                              atAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAg                              >  2:243591/1‑149 (MQ=255)
                                                                                                              atAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAg                              >  5:79604/1‑149 (MQ=255)
                                                                                                              atAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAg                              <  4:6152/149‑1 (MQ=255)
                                                                                                               tAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAgc                             >  2:129676/1‑149 (MQ=255)
                                                                                                                 gCTCAGTGCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACAGTCATCCGAgcgc                           >  8:142863/1‑149 (MQ=255)
                                                                                                                         ccTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGTATCACCGTCATCCGAGCGcc                          <  2:376612/142‑1 (MQ=255)
                                                                                                                         ccTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGTATCACCGTCATCCGAGCGcc                          >  1:376612/1‑142 (MQ=255)
                                                                                                                         ccTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCcagcag                    <  4:151424/148‑1 (MQ=255)
                                                                                                                         ccTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCagcagc                   <  2:398630/149‑1 (MQ=255)
                                                                                                                             cAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGc               <  1:243591/149‑1 (MQ=255)
                                                                                                                                gCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAg                                                                    <  2:95981/93‑1 (MQ=255)
                                                                                                                                gCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTAGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAg                                                                    >  1:95981/1‑93 (MQ=255)
                                                                                                                                       ttGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGCCGACCTAATc     >  2:290280/1‑149 (MQ=255)
                                                                                                                                          gAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGCCGACCTCATCGGc  <  3:5216/149‑1 (MQ=255)
                                                                                                                                          gAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGCCGACCTCATCGGc  >  4:5216/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TCTTCATTCTCCAGGCTATGGACCTCGGCTTCGTCCTTTAGATTCACCCCGGAAAGCTGGCGCCGGCAGGCTTCGCGCATCAAGTACAACAGGCGGTGCGCCGCCATGCCATAGCTCAGTCCCTCCAGGCGGATATTGGAAATGCAGTTGCGGTAGGCATCGGTCAGGCCGACCTTGGGTGCGTAAGTGAAATACAGGCCAAGGCTGTCGGGCGAGCTCAGCCCCGGGCGCTCGCCGATCAGCATCACCGTCATCCGAGCGCCCAGCAGCTCGCCGACCTCATCGGC  >  NC_002947/628019‑628305

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