Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,430,496 G→A 6.1% G89D (GGC→GAC)  PP_3925 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,430,4960GA 6.1% 198.7 / 2.1 66G89D (GGC→GAC) PP_3925hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (2/2);  total (32/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 = 9.82e-01

TCCTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGACGCCAGCCGGACCG  >  NC_002947/4430353‑4430639
                                                                                                                                               |                                                                                                                                               
tCCTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTc                                                                                                                                             <  2:115256/148‑1 (MQ=255)
 ccTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTcc                                                                                                                                            <  4:71806/148‑1 (MQ=255)
  cTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTcccc                                                                                                                                          >  1:228241/1‑149 (MQ=255)
  cTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTcccc                                                                                                                                          >  3:306155/1‑149 (MQ=255)
        gcgcCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGc                                                                                                                                    <  5:217491/149‑1 (MQ=255)
             tCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGcccc                                                                                                                                 >  8:58262/1‑147 (MQ=255)
               cAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAgg                                                                                                                             <  3:231858/149‑1 (MQ=255)
               cAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAg                                                                                                                              <  2:228241/148‑1 (MQ=255)
                      acctacctGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGTCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTAc                                                                                                                        <  6:211072/147‑1 (MQ=255)
                        ctacctGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGca                                                                                                                     >  1:198885/1‑148 (MQ=255)
                          acctGGTCAGCTACCAGGGCCTCGACTTCGTGCTGCGTCGACCAACCTTCGGTCAAAAGTCCAAATCGGCCAACGACATGGGTCGCCACTAACGCCTCCTCAACCACCTGCACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCacac                                                                                                                  >  5:327682/1‑149 (MQ=255)
                               gTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCAc                                                                                                             <  4:306155/149‑1 (MQ=255)
                                tCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCAcc                                                                                                            <  2:198885/149‑1 (MQ=255)
                                    ctacaggggCCGCGACTTCGTTCTGGGTCGACCACGCTTCGGTCAAAAGGCCAAATAGGCCCACGTCATGGGTCGCGAGTACCGCATCGTGAACCAACTGATCATAGTCTTCCCCTACTGCCCCAAGGCCTAGGCACACTGCACCg                                                                                                           <  6:327682/140‑1 (MQ=255)
                                        ccGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCgg                                                                                                    >  4:321997/1‑149 (MQ=255)
                                         cGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGc                                                                                                   <  6:155759/149‑1 (MQ=255)
                                          gggCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAg                                                                                                                              >  2:280853/1‑121 (MQ=255)
                                          gggCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAg                                                                                                                              <  1:280853/121‑1 (MQ=255)
                                             ccGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGcccc                                                                                                                                 <  1:261797/115‑1 (MQ=255)
                                             ccGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGcccc                                                                                                                                 >  2:261797/1‑115 (MQ=255)
                                                   cTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAATGCCTACGCACACTGCACCGACACCGGCCTGATcggcg                                                                                         >  4:181387/1‑149 (MQ=255)
                                                    ttCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATcggcgg                                                                                        >  8:38255/1‑149 (MQ=255)
                                                      cGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCggcggc                                                                                       <  3:321997/148‑1 (MQ=255)
                                                      cGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCggcggc                                                                                       <  7:38255/148‑1 (MQ=255)
                                                         gctgcGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCTGcc                                                                                                  >  1:332987/1‑134 (MQ=255)
                                                         gctgcGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCTGcc                                                                                                  <  2:332987/134‑1 (MQ=255)
                                                              gTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAt                                                                              <  7:58262/149‑1 (MQ=255)
                                                                     accCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATg                                                                        >  8:333954/1‑148 (MQ=255)
                                                                      cccTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCACCCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATgg                                                                       >  5:254024/1‑148 (MQ=255)
                                                                       ccTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGa                                                                                     <  6:322929/133‑1 (MQ=255)
                                                                       ccTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGa                                                                                     >  5:322929/1‑133 (MQ=255)
                                                                       ccTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGaa                                                                     >  3:228808/1‑149 (MQ=255)
                                                                           cGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACgg                                                                  <  7:197038/148‑1 (MQ=255)
                                                                             gTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGt                                                                >  1:175329/1‑148 (MQ=255)
                                                                                aaaaGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAtg                                                                             <  2:269278/132‑1 (MQ=255)
                                                                                    ggCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAATGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGca                                                        <  3:181387/149‑1 (MQ=255)
                                                                                    ggCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGca                                                        <  6:254024/149‑1 (MQ=255)
                                                                                      ccAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGt                                                                <  4:35733/139‑1 (MQ=255)
                                                                                      ccAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGt                                                                >  3:35733/1‑139 (MQ=255)
                                                                                      ccAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCGCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGGGATGGACCGGGTCAAGGGcag                                                       >  6:358218/1‑147 (MQ=255)
                                                                                           tCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCt                                                 <  4:228808/149‑1 (MQ=255)
                                                                                              gCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGt                                                                <  6:41678/131‑1 (MQ=255)
                                                                                              gCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGt                                                                >  5:41678/1‑131 (MQ=255)
                                                                                                ccACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGaa                                                                     <  7:88603/124‑1 (MQ=255)
                                                                                                ccACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGaa                                                                     >  8:88603/1‑124 (MQ=255)
                                                                                                ccACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGACTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCt                                            >  3:121362/1‑149 (MQ=255)
                                                                                                        tGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATg                                                                        <  7:273730/113‑1 (MQ=255)
                                                                                                        tGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATg                                                                        >  8:273730/1‑113 (MQ=255)
                                                                                                               cgAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGcc                              >  6:283411/1‑148 (MQ=255)
                                                                                                                gAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAtgt                                                                            >  5:43536/1‑101 (MQ=255)
                                                                                                                gAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAtgt                                                                            <  6:43536/101‑1 (MQ=255)
                                                                                                                gAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCt                                                 <  6:126954/128‑1 (MQ=255)
                                                                                                                gAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCt                                                 >  5:126954/1‑128 (MQ=255)
                                                                                                                   tACCGCATCCTCAACCAACTGAACAGCGACTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACt                         <  4:121362/149‑1 (MQ=255)
                                                                                                                        cATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAAcc                    >  1:346222/1‑149 (MQ=255)
                                                                                                                        cATCCTCAACCAACTGAACAGCGACTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGcatca                                                     >  6:251637/1‑116 (MQ=255)
                                                                                                                        cATCCTCAACCAACTGAACAGCGACTTCCCCTACTGCACCAAGGCATACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGcatca                                                     <  5:251637/116‑1 (MQ=255)
                                                                                                                            cTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAtgtg                                                                           <  1:63631/90‑1 (MQ=255)
                                                                                                                            cTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTAtgtg                                                                           >  2:63631/1‑90 (MQ=255)
                                                                                                                             tCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGc                               <  8:151053/133‑1 (MQ=255)
                                                                                                                             tCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGc                               >  7:151053/1‑133 (MQ=255)
                                                                                                                             tCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGAc               >  1:314450/1‑149 (MQ=255)
                                                                                                                               aaccaacTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGACGc             >  8:34351/1‑149 (MQ=255)
                                                                                                                                          cAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGACGCCAGCCGGACCg  <  5:283411/149‑1 (MQ=255)
                                                                                                                                           aGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGAc               >  7:286157/1‑135 (MQ=255)
                                                                                                                                           aGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGAc               <  8:286157/135‑1 (MQ=255)
                                                                                                                                               gCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCg                                          >  4:187398/1‑104 (MQ=255)
                                                                                                                                               gCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCg                                          <  3:187398/104‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCCTGGCGGCGCCTCCAACCTCACCTACCTGGTCAGCTACCCGGGCCGCGACTTCGTGCTGCGTCGACCACCCTTCGGTCAAAAGGCCAAATCGGCCCACGACATGGGTCGCGAGTACCGCATCCTCAACCAACTGAACAGCGGCTTCCCCTACTGCCCCAAGGCCTACGCACACTGCACCGACACCGGCCTGATCGGCGGCGAGTTCTATGTGATGGAACGGGTCAAGGGCATCATCCTGCGCTCGGACATCCCAGCCGAACTGAACCTGGACGCCAGCCGGACCG  >  NC_002947/4430353‑4430639

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