Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,354 T→G 6.1% E51A (GAG→GCG)  eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,3540TG 6.1% 208.1 / 2.4 76E51A (GAG→GCG) ecoecotin
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base G (2/3);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGC  >  NC_002947/3456205‑3456503
                                                                                                                                                     |                                                                                                                                                     
tGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACt                                                                                                                                                       >  1:774260/1‑150 (MQ=255)
    ggccgggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACtcgtc                                                                                                                                                   <  1:360283/150‑1 (MQ=255)
       cgggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCACAGGGCCTGAAGAGTGCAGTCGACCTGCAGCGTCTCGACAGCGAGGACCAGCCGCTTGTAGGCCGACTCGTCGGc                                                                                                                                                >  1:337496/1‑150 (MQ=255)
         ggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCt                                                                                                                                              <  2:389617/150‑1 (MQ=255)
         ggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCt                                                                                                                                              <  1:666085/150‑1 (MQ=255)
           ccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTgg                                                                                                                                            <  1:284975/150‑1 (MQ=255)
             gCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGc                                                                                                                                          >  1:1738244/1‑150 (MQ=255)
              cTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGcc                                                                                                                                         >  1:569738/1‑150 (MQ=255)
                   gTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGTGCCGGCAg                                                                                                                                    >  2:1156061/1‑150 (MQ=255)
                    ttGTCGAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGGGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAgg                                                                                                                                   <  1:122787/150‑1 (MQ=255)
                      ggccAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCATCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTg                                                                                                                                 <  2:669466/148‑1 (MQ=255)
                        ccAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGTGCCGGCAGGTGGa                                                                                                                               <  1:1156061/150‑1 (MQ=255)
                         cAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGAt                                                                                                                              >  1:1515927/1‑150 (MQ=255)
                           ggCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATaa                                                                                                                            <  2:1521498/150‑1 (MQ=255)
                            gCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCCGGTTGCCGCCCAGGCGCTGACGATTGCCGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCAGACTCGCGGGCCGGGGCCGGCAGGGGGATCAc                                                                                                                           >  1:911724/1‑150 (MQ=255)
                            gCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAAc                                                                                                                           >  1:495909/1‑150 (MQ=255)
                             cGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACCCGCCGGCCTGGGCCGGCAGGTGGATAAcc                                                                                                                          <  1:1022705/150‑1 (MQ=255)
                               gtagtaATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTGTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACGAGTCGGCCtataaa                                                                                                                                         <  1:883906/133‑6 (MQ=255)
                               gtagtaATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACGAGTCGGCCtataaa                                                                                                                                         >  2:883906/1‑128 (MQ=255)
                                   taATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGTGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCgg                                                                                                                    >  1:1219740/1‑150 (MQ=255)
                                    aaTTGTAGCCCCAGCCTTCCAGGGTGTGCGCTTCCACGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCggg                                                                                                                   <  1:6545/150‑1 (MQ=255)
                                     aTTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGt                                                                                                                  <  2:1738244/150‑1 (MQ=255)
                                        gTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGaa                                                                                                               <  2:1642429/150‑1 (MQ=255)
                                            ccccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCTGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCt                                                                                                           <  1:1551858/150‑1 (MQ=255)
                                              ccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAACAttt                                                                                                         >  2:1458176/1‑150 (MQ=255)
                                               cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCtttt                                                                                                        <  2:266188/150‑1 (MQ=255)
                                                        aGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCg                                                                                               >  2:1122001/1‑150 (MQ=255)
                                                         gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgg                                                                                              >  2:1770186/1‑150 (MQ=255)
                                                         gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCGGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCATCgg                                                                                              <  2:254665/150‑1 (MQ=255)
                                                         gcgtgtgtCCTTCCAGGTGGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATCGGCCGACTCGTCGGCCTCGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgg                                                                                              <  2:911724/148‑1 (MQ=255)
                                                          ggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCggg                                                                                             <  2:1179711/150‑1 (MQ=255)
                                                           gtgtgtGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgggg                                                                                            <  1:187442/150‑1 (MQ=255)
                                                              tgtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACtcgt                                                                                                                                                    <  1:1128668/91‑1 (MQ=255)
                                                              tgtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACtcgt                                                                                                                                                    >  2:1128668/1‑91 (MQ=255)
                                                                   ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATaa                                                                                                                            <  2:1458244/110‑1 (MQ=255)
                                                                   ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATaa                                                                                                                            >  1:1458244/1‑110 (MQ=255)
                                                                   ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgcgc                                                                                    >  2:1645134/1‑150 (MQ=255)
                                                                   ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAAGATTTCCAAATTATAGGCCGACTCGTCGGCCTGGGCCTGCAGGTGGATAACCTGCCGGGTGAAACCTTTTTCCGCCTCGGGGTAGGgcgg                                                                                    >  2:655837/1‑149 (MQ=255)
                                                                    tCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTGTATATGCGGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGcgcg                                                                                   >  2:161957/1‑150 (MQ=255)
                                                                     ccAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                                  >  2:976304/1‑150 (MQ=255)
                                                                      cAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGAt                                                                                 <  1:533397/150‑1 (MQ=255)
                                                                       aGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGTGGTTTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATg                                                                                >  1:592909/1‑150 (MQ=255)
                                                                          ttGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTcc                                                                             <  2:495909/150‑1 (MQ=255)
                                                                          ttGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTcc                                                                             <  2:774260/150‑1 (MQ=255)
                                                                               gccCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAg                                                                        >  2:661342/1‑150 (MQ=255)
                                                                                cccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAgg                                                                       <  2:569738/150‑1 (MQ=255)
                                                                                cccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAgg                                                                       <  1:976304/150‑1 (MQ=255)
                                                                                 ccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGcc                                                                                                                                               >  2:1491670/1‑77 (MQ=255)
                                                                                 ccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGcc                                                                                                                                               <  1:1491670/77‑1 (MQ=255)
                                                                                 ccAGGCGATGACGATTGCAGACGACATGCGTGTTCTGGCCAGCGAGGATTACAACTTTATAGTATGCAGTTTCGTCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGc                                                                      <  1:527160/150‑1 (MQ=255)
                                                                                     gcgcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggc                                                                  <  1:1122001/150‑1 (MQ=255)
                                                                                      cgcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGTCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggga                                                                 <  1:1341601/150‑3 (MQ=255)
                                                                                       gcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTc                                                                                                      <  1:335788/112‑1 (MQ=255)
                                                                                       gcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTc                                                                                                      >  2:335788/1‑112 (MQ=255)
                                                                                                           tGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCa                                                                         >  2:1316319/1‑121 (MQ=255)
                                                                                                           tGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCa                                                                         <  1:1316319/121‑1 (MQ=255)
                                                                                                                cggattGCCAGCGGGGATTTTCAGTTTATAGGCCGACTTGGCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAg                                       <  2:338033/146‑1 (MQ=255)
                                                                                                                  tCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGc                                     >  2:1923190/1‑150 (MQ=255)
                                                                                                                       ccAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAgg                                >  1:63049/1‑150 (MQ=255)
                                                                                                                        cAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAggg                               >  2:1691152/1‑150 (MQ=255)
                                                                                                                          gCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGcc                             <  2:600285/150‑1 (MQ=255)
                                                                                                                           cGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgc                                                                                      <  1:1119470/92‑1 (MQ=255)
                                                                                                                           cGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGCCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgc                                                                                      >  2:1119470/1‑92 (MQ=255)
                                                                                                                               gATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGc                              >  1:1648060/1‑144 (MQ=255)
                                                                                                                               gATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGc                              <  2:1648060/144‑1 (MQ=255)
                                                                                                                                aTTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATc                       <  1:526168/150‑1 (MQ=255)
                                                                                                                                  ttCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgc                     <  2:1219740/150‑1 (MQ=255)
                                                                                                                                  ttCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgc                     <  2:1174911/150‑1 (MQ=255)
                                                                                                                                     cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGt                  >  2:1145428/1‑150 (MQ=255)
                                                                                                                                       gTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCACGGGGAAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCa                >  2:1673894/1‑150 (MQ=255)
                                                                                                                                           ataGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCa                <  1:1697906/146‑1 (MQ=255)
                                                                                                                                           ataGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgg            >  1:362044/1‑150 (MQ=255)
                                                                                                                                           ataGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGGGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCa                >  2:1697906/1‑146 (MQ=255)
                                                                                                                                            taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg           <  2:337496/150‑1 (MQ=255)
                                                                                                                                            caGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCgg                   >  2:547156/2‑142 (MQ=255)
                                                                                                                                            caGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCgg                   <  1:547156/141‑1 (MQ=255)
                                                                                                                                             aGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgggg          <  2:362044/150‑1 (MQ=255)
                                                                                                                                                  gACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAggg     <  1:161957/150‑1 (MQ=255)
                                                                                                                                                  gACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGCATGGAGCGTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGGCATGGGGGTAggg     >  1:1902856/1‑150 (MQ=255)
                                                                                                                                                   aCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGAAGGCCAAGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGgc    >  1:433526/1‑150 (MQ=255)
                                                                                                                                                     tcgtcgGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGgcgc  >  2:98825/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGC  >  NC_002947/3456205‑3456503

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