Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,546,790 C→T 5.6% P616P (CCC→CCT PP_2238 → metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,546,7900CT 5.6% 202.4 / 2.4 71P616P (CCC→CCTPP_2238metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base T (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

AAGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCC  >  NC_002947/2546644‑2546937
                                                                                                                                                  |                                                                                                                                                   
aaGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTg                                                                                                                                                  <  1:56774/150‑1 (MQ=255)
      gCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCaa                                                                                                                                            <  1:147563/150‑1 (MQ=255)
       ccTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                                                                                                           <  1:1231436/150‑1 (MQ=255)
       ccTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                                                                                                           <  1:1231675/150‑1 (MQ=255)
       ccTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                                                                                                           <  1:1458613/150‑1 (MQ=255)
       ccTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTAAATGAGCGCCTGGCGCCCTTGCTCCAAg                                                                                                                                           <  2:223041/150‑1 (MQ=255)
         tttGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAggg                                                                                                                                         >  2:1736122/1‑150 (MQ=255)
                aCTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTTGTGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCTTTGCTCCAAGGGGATgcgc                                                                                                                                  >  1:1095225/1‑150 (MQ=255)
                aCTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATTGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGGGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCTTTGCTCCAAGGGGATgcgc                                                                                                                                  >  1:1095432/1‑150 (MQ=255)
                  ttCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTa                                                                                                                                >  2:1775995/1‑150 (MQ=255)
                    ccTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGc                                                                                                                              <  2:1130231/150‑1 (MQ=255)
                     cTGCAGTCCGAAACCTAGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGACTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCGAGGGGATGCGCTAGcc                                                                                                                             <  2:975152/150‑1 (MQ=255)
                       gCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctg                                                                                                                           >  1:1797755/1‑150 (MQ=255)
                              gAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGga                                                                                                                     >  2:33698/1‑149 (MQ=255)
                              gAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGga                                                                                                                     <  1:33698/149‑1 (MQ=255)
                               aaaCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGgaga                                                                                                                   >  2:1642360/1‑150 (MQ=255)
                               aaaCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGgaga                                                                                                                   >  2:287726/1‑150 (MQ=255)
                                aaCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGAt                                                                                                                  <  2:1249008/150‑1 (MQ=255)
                                 aCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATg                                                                                                                 >  2:1623587/1‑150 (MQ=255)
                                  ccTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGc                                                                                                                <  2:645526/150‑1 (MQ=255)
                                   taagACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCg                                                                                                               >  1:544199/4‑150 (MQ=255)
                                    ttGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctgg                                                                                                                          <  1:1201933/138‑1 (MQ=255)
                                    ttGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctgg                                                                                                                          >  2:1201933/1‑138 (MQ=255)
                                    ttGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGt                                                                                                              >  1:1628522/1‑150 (MQ=255)
                                     tGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTa                                                                                                             <  1:1736122/150‑1 (MQ=255)
                                         gCTGTGCCCTATCGACACCCGCTGTCTGTTGTCCGAGTGTCGGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCATGCGTACATGAGCGCCTGGCGCCTTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGc                                                                                                         <  2:1095432/150‑1 (MQ=255)
                                         gCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGGGTGGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCTTGCGTACATGAGCGCCTGGCGCCTTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGc                                                                                                         <  2:1095225/150‑1 (MQ=255)
                                                      gACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGc                                                                                            >  1:1634541/1‑150 (MQ=255)
                                                         acCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTaa                                                                                               >  2:153166/1‑144 (MQ=255)
                                                         acCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTaa                                                                                               <  1:153166/144‑1 (MQ=255)
                                                                gtctgtTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGaaaaa                                                                                  >  1:1530177/1‑150 (MQ=255)
                                                                   tgtTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAggg                                                                               <  2:1770811/150‑1 (MQ=255)
                                                                        cccGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGtt                                                                          >  1:1369430/1‑150 (MQ=255)
                                                                         ccGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAggg                                                                                                                                         <  2:605971/86‑1 (MQ=255)
                                                                         ccGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAggg                                                                                                                                         >  1:605971/1‑86 (MQ=255)
                                                                           gAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                       >  2:563583/1‑150 (MQ=255)
                                                                                   gACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc                                                               >  1:1735761/1‑150 (MQ=255)
                                                                                        aggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtttt                                                          <  2:1369430/150‑1 (MQ=255)
                                                                                               ggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa                                                     <  1:660194/148‑1 (MQ=255)
                                                                                               ggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa                                                     >  2:660194/1‑148 (MQ=255)
                                                                                                gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAttt                                                  <  1:287726/150‑1 (MQ=255)
                                                                                                gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGCGCAGTTTGCGAACTGCCCTTTTCTCTAttt                                                  <  1:1642360/150‑1 (MQ=255)
                                                                                                  ggAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACAGTCATTTTCTCTATTTAc                                                >  2:1295645/1‑150 (MQ=255)
                                                                                                   gAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTACACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACa                                               <  2:544199/150‑1 (MQ=255)
                                                                                                     gTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAga                                             <  2:1634541/150‑1 (MQ=255)
                                                                                                          tGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatg                                        >  1:249340/1‑150 (MQ=255)
                                                                                                          tGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatg                                        >  1:877762/1‑150 (MQ=255)
                                                                                                                gCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggct                                                                                                                        >  1:1430270/1‑64 (MQ=255)
                                                                                                                gCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggct                                                                                                                        <  2:1430270/64‑1 (MQ=255)
                                                                                                                gCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAggg                                                                               >  1:1407535/1‑105 (MQ=255)
                                                                                                                gCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAggg                                                                               <  2:1407535/105‑1 (MQ=255)
                                                                                                                 cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAtt                                                   <  1:378274/132‑1 (MQ=255)
                                                                                                                 cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAtt                                                   >  2:378274/1‑132 (MQ=255)
                                                                                                                  tACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTCTTTACAGATGATGATCATGCt                                >  2:280277/1‑150 (MQ=255)
                                                                                                                          cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                       <  2:1236393/103‑1 (MQ=255)
                                                                                                                          cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                       >  1:1236393/1‑103 (MQ=255)
                                                                                                                          cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTgg                        >  2:795324/1‑150 (MQ=255)
                                                                                                                          cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTgg                        >  2:672981/1‑150 (MQ=255)
                                                                                                                           tGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                     >  2:504584/1‑136 (MQ=255)
                                                                                                                           tGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                     <  1:504584/136‑1 (MQ=255)
                                                                                                                              gTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggc                                                                                                                         <  1:1141358/49‑1 (MQ=255)
                                                                                                                              gTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggc                                                                                                                         >  2:1141358/1‑49 (MQ=255)
                                                                                                                              gTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc                                                               <  1:133306/107‑1 (MQ=255)
                                                                                                                              gTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc                                                               >  2:133306/1‑107 (MQ=255)
                                                                                                                              gTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCATACTGCACCGCTGTGACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCGGgtgtc                    >  2:524938/1‑149 (MQ=255)
                                                                                                                              gTACATGAACGAAAGGAGCCCTATCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTa                    >  2:556389/1‑150 (MQ=255)
                                                                                                                               tACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGGTCGTTTTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGATGTGCTCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAg                   <  1:347514/150‑1 (MQ=255)
                                                                                                                                      gcgcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCccgccg            >  1:1063768/1‑150 (MQ=255)
                                                                                                                                       cgcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCgg           <  2:1530177/150‑1 (MQ=255)
                                                                                                                                         ccTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGt         <  1:563583/150‑1 (MQ=255)
                                                                                                                                               gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc   >  1:511130/1‑150 (MQ=255)
                                                                                                                                                cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGcc  <  2:1063768/150‑1 (MQ=255)
                                                                                                                                                 ccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtttt                                                          >  1:1617299/1‑93 (MQ=255)
                                                                                                                                                 ccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtttt                                                          <  2:1617299/93‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCC  >  NC_002947/2546644‑2546937

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