Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 112,103 C→T 6.3% P175L (CCG→CTG)  PP_0105 → cytochrome c oxidase assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947112,1030CT 6.3% 183.7 / 2.0 63P175L (CCG→CTG) PP_0105cytochrome c oxidase assembly protein
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base T (2/2);  total (33/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 = 5.76e-01

GAGTGCTTCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACGGTATTCG  >  NC_002947/111955‑112247
                                                                                                                                                    |                                                                                                                                                
gAGTGCTTCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACccgcc                                                                                                                                                  >  2:339752/1‑149 (MQ=255)
     cTTCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTc                                                                                                                                              <  7:130575/148‑1 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  7:281682/1‑149 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  8:38617/1‑149 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  3:331255/1‑149 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  3:123821/1‑149 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  4:248135/1‑149 (MQ=255)
      ttCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTAGACATCACCGCTCGCCACCCGCCGGTcgc                                                                                                                                            >  2:301631/1‑149 (MQ=255)
           cTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCAc                                                                                                                                        <  5:118079/148‑1 (MQ=255)
              tACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGt                                                                                                                                      <  4:93834/147‑1 (MQ=255)
                 ccAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCTCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTgcgg                                                                                                                                 >  7:92964/1‑149 (MQ=255)
                 ccAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTgcgg                                                                                                                                 >  5:253500/1‑149 (MQ=255)
                 ccAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTgcgg                                                                                                                                 >  3:21232/1‑149 (MQ=255)
                       ggTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAg                                                                                                                            >  2:140318/1‑148 (MQ=255)
                       ggTGCAGCAGCCCGGCGAACGCACCGAGATGCCGGTGCGCTTCACCCTCGACCGCGACCGGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGATCGACATCACCGCCCGCCACCCGCCGCTCGCGCACGCCGCGGCGCAgg                                                                                                                           <  4:123821/149‑1 (MQ=255)
                           ctgcAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGt                                                                                                                       >  5:122373/1‑149 (MQ=255)
                             gcAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCTGGTCGCGCACGTTGCGGCGCAGGACGTc                                                                                                                      >  4:156045/1‑148 (MQ=255)
                               aGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAg                                                                                                                   >  7:327303/1‑149 (MQ=255)
                                 cccGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCa                                                                                                                    <  1:334591/146‑1 (MQ=255)
                                 cccGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCAGCCCCTCGCGCACGTTGCGACGCAGGACGTCCa                                                                                                                    >  2:334591/1‑146 (MQ=255)
                                           gCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGtt                                                                                                         >  7:147743/1‑147 (MQ=255)
                                           gCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTaa                                                                                                       <  4:331255/149‑1 (MQ=255)
                                            cATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAg                                                                                                      <  3:248135/149‑1 (MQ=255)
                                            cATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAg                                                                                                      <  4:346948/149‑1 (MQ=255)
                                                gagaTGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGaa                                                                                                  <  8:327303/149‑1 (MQ=255)
                                                       cGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAAc                                                                                           <  8:147743/149‑1 (MQ=255)
                                                         gTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAg                                                                                         <  1:140318/149‑1 (MQ=255)
                                                             gcTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCTGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAAt                                                                                     <  3:156045/149‑1 (MQ=255)
                                                              cTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATg                                                                                    >  3:282720/1‑149 (MQ=255)
                                                              cTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATg                                                                                    <  1:339752/149‑1 (MQ=255)
                                                               ttCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATgg                                                                                   >  5:146233/1‑149 (MQ=255)
                                                                 cATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATgg                                                                                   <  6:122373/147‑1 (MQ=255)
                                                                  aTCGTCGACCTCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCaa                                                                                <  8:92964/149‑1 (MQ=255)
                                                                  aTCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCaa                                                                                <  4:282720/149‑1 (MQ=255)
                                                                     gtcgACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGt                                                                              >  4:307760/1‑148 (MQ=255)
                                                                           cgcgACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAg                                                                        >  6:254738/1‑148 (MQ=255)
                                                                             cgACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCTGGTCGCGCACGTTGCGGCGCAGGACGTcc                                                                                                                     <  4:281598/101‑1 (MQ=255)
                                                                             cgACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCTGGTCGCGCACGTTGCGGCGCAGGACGTcc                                                                                                                     >  3:281598/1‑101 (MQ=255)
                                                                             cgACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCac                                                                     >  4:291932/1‑149 (MQ=255)
                                                                                 cTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCactac                                                                  >  8:262774/1‑148 (MQ=255)
                                                                                         cAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCg                                                         <  3:307760/149‑1 (MQ=255)
                                                                                            cGTGAAGCACATGACACTGGACTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGTCGACCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGAACTACTATGATCCgtc                                                       >  3:249066/1‑146 (MQ=255)
                                                                                                       tGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATg                                           >  5:278899/1‑149 (MQ=255)
                                                                                                         acacTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGATGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGc                                         <  7:38617/149‑1 (MQ=255)
                                                                                                         acacTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGc                                         <  8:281682/149‑1 (MQ=255)
                                                                                                         acacTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGc                                         <  7:262774/149‑1 (MQ=255)
                                                                                                               gCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGatc                                    <  6:253500/148‑1 (MQ=255)
                                                                                                               gCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGAGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGatc                                    <  1:301631/148‑1 (MQ=255)
                                                                                                                 cTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGatcatc                                 >  2:117053/1‑149 (MQ=255)
                                                                                                                 cTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGa                                      <  8:356908/144‑1 (MQ=255)
                                                                                                                 cTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGa                                      >  7:356908/1‑144 (MQ=255)
                                                                                                                           ttCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCg                       >  7:367197/1‑149 (MQ=255)
                                                                                                                                 aTCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGt                 >  8:4315/1‑149 (MQ=255)
                                                                                                                                 aTCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGt                 <  4:21232/149‑1 (MQ=255)
                                                                                                                                   cACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGtt                >  6:307438/1‑148 (MQ=255)
                                                                                                                                   cACCGCTCGCCACCCGCCCGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGtt                >  2:248960/1‑148 (MQ=255)
                                                                                                                                    aCCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTtca              >  1:252428/1‑149 (MQ=255)
                                                                                                                                     ccGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTtcat             >  7:167670/1‑149 (MQ=255)
                                                                                                                                      cGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTtcatc            <  3:291932/149‑1 (MQ=255)
                                                                                                                                        cTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGat                                     <  8:112353/122‑1 (MQ=255)
                                                                                                                                        cTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGat                                     >  7:112353/1‑122 (MQ=255)
                                                                                                                                           gCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACg         <  5:307438/147‑1 (MQ=255)
                                                                                                                                           gCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACg         <  8:367197/147‑1 (MQ=255)
                                                                                                                                              aCCCGCCCGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACGGTAtt    <  1:248960/149‑1 (MQ=255)
                                                                                                                                                ccgccgGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACGGTATTCg  <  6:146233/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAGTGCTTCTGCTTTACCCAGCAGGTGCTGCAGCCCGGCGAACGCATCGAGATGCCGGTGCGCTTCATCGTCGACCGCGACCTGCCGGCCAGCGTGAAGCACCTGACACTGGCCTACACCTTGTTCGACATCACCGCTCGCCACCCGCCGGTCGCGCACGTTGCGGCGCAGGACGTCCAGGGCGCCCGTTAAGGGAAGGAGAACAGCAATGGCAAGTCATGAGCACTACTATGTTCCGGCGCAGAGCAAATGGCCGATCATCGCCACCATCGGCATGTTCATCACGGTATTCG  >  NC_002947/111955‑112247

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