Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,355 C→A 8.2% E51* (GAG→TAG) ‡ eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,3550CA 8.2% 179.7 / 5.7 61E51* (GAG→TAG) ‡ecoecotin
Reads supporting (aligned to +/- strand):  ref base C (22/34);  new base A (3/2);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

GGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGC  >  NC_002947/3456209‑3456503
                                                                                                                                                  |                                                                                                                                                    
ggccgggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACtcgtc                                                                                                                                                   >  1:190741/1‑150 (MQ=255)
    gggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGcc                                                                                                                                               >  2:1264058/1‑150 (MQ=255)
     ggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCt                                                                                                                                              <  2:287426/150‑1 (MQ=255)
         gCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGAGGTGACGATTGCAGTCGACCTGCAGCCACTTGCTCGAAACGAGTTGCGGTTTATAGGCCGACTCGTCGGCCTGGGc                                                                                                                                          >  1:672946/1‑150 (MQ=255)
               gTTGTCCAGGCGGTAGTAATTGTAGCCGCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCGTAGCATCTTACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAg                                                                                                                                    <  1:1157267/150‑1 (MQ=255)
               gTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAg                                                                                                                                    >  2:562560/1‑150 (MQ=255)
                  gTCCAGGCGGTAGTAATTGTAGCCCCAGCCTGCCAGGGTGTGTGCTTCAAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTg                                                                                                                                 <  1:1020968/150‑1 (MQ=255)
                     cAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGAt                                                                                                                              <  1:258435/150‑1 (MQ=255)
                      aGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATa                                                                                                                             >  2:817518/1‑150 (MQ=255)
                           gtagtaATTGTAGCCCCAGCCTTCCAGGGTGAGTACCTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGTCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTg                                                                                                                        >  1:690493/1‑150 (MQ=255)
                                 aTTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGt                                                                                                                  <  1:106163/150‑1 (MQ=255)
                                           cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGt                                                                                                                                  <  2:149373/124‑1 (MQ=255)
                                           cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGt                                                                                                                                  >  1:149373/1‑124 (MQ=255)
                                             gCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTc                                                                                                      <  2:190741/150‑1 (MQ=255)
                                             gCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTc                                                                                                      <  1:562560/150‑1 (MQ=255)
                                              ccTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTcc                                                                                                     <  1:700359/150‑1 (MQ=255)
                                              ccTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTcc                                                                                                     <  1:1264058/150‑1 (MQ=255)
                                               cTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCg                                                                                                    >  2:645560/1‑150 (MQ=255)
                                                ttCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGc                                                                                                   >  2:255728/1‑150 (MQ=255)
                                                ttCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAATAGGACTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGc                                                                                                   <  1:1387087/150‑1 (MQ=255)
                                                   gaGGGTGTGGGGTACCCGGTTGCCGATCGGGAGACGTCAATTTCATTCGGCCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTc                                                                                                <  2:373768/149‑1 (MQ=255)
                                                     gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTg                                                                                                                                 >  1:1118576/1‑115 (MQ=255)
                                                     gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTg                                                                                                                                 <  2:1118576/115‑1 (MQ=255)
                                                     gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgg                                                                                              <  2:976295/150‑1 (MQ=255)
                                                     gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgg                                                                                              <  2:360626/150‑1 (MQ=255)
                                                      ggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCggg                                                                                             >  2:184369/1‑150 (MQ=255)
                                                          tgtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGACTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAg                                                                                         >  1:831265/1‑150 (MQ=255)
                                                            tgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAggg                                                                                       >  2:854510/1‑150 (MQ=255)
                                                               ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgcgc                                                                                    <  1:645560/150‑1 (MQ=255)
                                                               ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTCCCGGGTGAAGCCTTTTTCCGCCTCGGGGGAGGgcgc                                                                                    >  1:936896/1‑150 (MQ=255)
                                                                 ccAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGTGAGGCTTTCCAGTTTATAGGCCGGCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                                  <  2:1367203/150‑1 (MQ=255)
                                                                 ccAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCAAGCGCGGGTAGCCAGTTTATCGTCGTACGCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                                  >  1:884151/1‑150 (MQ=255)
                                                                        gccgccCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCtt                                                                           <  1:975517/150‑1 (MQ=255)
                                                                             ccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCGTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTGGAGGc                                                                      >  1:494059/1‑150 (MQ=255)
                                                                               aGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGCAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTg                                                                    <  2:714750/150‑1 (MQ=255)
                                                                                   gcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggcgg                                                                <  1:565225/150‑1 (MQ=255)
                                                                                   gcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTGTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggcgg                                                                >  2:945675/1‑150 (MQ=255)
                                                                                     tGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGaa                                                                                                               >  1:1396231/1‑101 (MQ=255)
                                                                                     tGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGaa                                                                                                               <  2:1396231/101‑1 (MQ=255)
                                                                                            tGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAg                                                       <  2:522710/150‑1 (MQ=255)
                                                                                            tGCAGGCGACCTGCAGCTTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAg                                                       <  1:21043/150‑1 (MQ=255)
                                                                                               aGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggcgg                                                                >  1:90465/1‑138 (MQ=255)
                                                                                               aGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggcgg                                                                <  2:90465/138‑1 (MQ=255)
                                                                                                gTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCTTGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgc                                                   <  1:1045499/150‑1 (MQ=255)
                                                                                                        gCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGcggcggcgg                                           <  1:817518/150‑1 (MQ=255)
                                                                                                          aGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcg                                               >  2:701960/1‑144 (MQ=255)
                                                                                                          aGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcg                                               <  1:701960/144‑1 (MQ=255)
                                                                                                          aGCGTCTTCCCAGCGAGTATTTCCAGTTTATAAGCTGACTCGTCTTCCTTGGCCGGCAGTTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGcc                                         <  2:515477/150‑1 (MQ=255)
                                                                                                                    cAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAggg                               <  2:1359013/150‑1 (MQ=255)
                                                                                                                         aGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGa                          >  2:991012/1‑150 (MQ=255)
                                                                                                                                    tctagaGGCCGACTCCTCGGCCTGGGCCGCCAGGAGGATAACCTGCCGGGTGAAGCCGTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCCGGGCCAGAATGGCGGTCAt               <  2:690493/145‑1 (MQ=255)
                                                                                                                                      tataGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg             <  1:534294/150‑1 (MQ=255)
                                                                                                                                       ccaGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                                  <  1:551791/78‑1 (MQ=255)
                                                                                                                                       ccaGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                                  >  2:551791/3‑80 (MQ=255)
                                                                                                                                       ccaGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgg            >  2:670462/3‑150 (MQ=255)
                                                                                                                                        caGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcggc                                             >  1:365663/2‑116 (MQ=255)
                                                                                                                                        caGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcggc                                             <  2:365663/115‑1 (MQ=255)
                                                                                                                                           gCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCtt                                                                           <  2:1352702/83‑1 (MQ=255)
                                                                                                                                           gCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCtt                                                                           >  1:1352702/1‑83 (MQ=255)
                                                                                                                                           gCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggc                                                                  <  2:534402/92‑1 (MQ=255)
                                                                                                                                           gCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggc                                                                  >  1:534402/1‑92 (MQ=255)
                                                                                                                                               aCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGgc    <  1:991012/150‑1 (MQ=255)
                                                                                                                                                 tcgtcgGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGgcgc  <  1:184369/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGC  >  NC_002947/3456209‑3456503

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