Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,561,243 G→C 9.2% I102M (ATC→ATG)  PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,561,2430GC 9.2% 186.0 / 7.6 63I102M (ATC→ATG) PP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (2/4);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGC  >  NC_002947/4561100‑4561387
                                                                                                                                               |                                                                                                                                                
gCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCgg                                                                                                                                             <  5:216132/149‑1 (MQ=255)
  aaCCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGgcg                                                                                                                                           >  8:3906/1‑149 (MQ=255)
     cACTGTCTGCATGATGACCTGACTTTCACGTATAATCAGCAAGCACAAACCACGCAAGATCGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgct                                                                                                                                         <  7:330603/148‑2 (MQ=255)
     cACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGGACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctg                                                                                                                                        <  5:139609/149‑3 (MQ=255)
      aCTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCAATCGATCCCACGGGAGTTCCCACCGCCCCCACCACAACGCCGATCAACGgcgctg                                                                                                                                        >  8:330603/1‑146 (MQ=255)
       cTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgtg                                                                                                                                      <  4:141061/149‑5 (MQ=255)
             gCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAggg                                                                                                                                  <  1:61284/147‑1 (MQ=255)
             gCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAggg                                                                                                                                  >  2:61284/1‑147 (MQ=255)
                ggATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatc                                                                                                                             <  5:366021/149‑1 (MQ=255)
                     aCCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggc                                                                                                                         <  7:3906/148‑1 (MQ=255)
                      ccTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAgg                                                                                                                                   >  5:333622/1‑137 (MQ=255)
                      ccTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAgg                                                                                                                                   <  6:333622/137‑1 (MQ=255)
                      ccTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                       <  7:250083/149‑1 (MQ=255)
                         gACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGcc                                                                                                                    >  4:96348/1‑149 (MQ=255)
                            tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccac                                                                                                                  >  8:119051/1‑148 (MQ=255)
                                  gTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                     >  4:171515/1‑139 (MQ=255)
                                  gTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                     <  3:171515/139‑1 (MQ=255)
                                     tCAACAGCAAGCACTAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGCAGTGCCCACCGCCACCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACAcc                                                                                                        >  4:46438/1‑149 (MQ=255)
                                        aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCa                                                                                                      >  7:92162/1‑148 (MQ=255)
                                                cacaAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCa                                                                                             >  7:195730/1‑149 (MQ=255)
                                                 acaAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCAt                                                                                            >  5:180599/1‑149 (MQ=255)
                                                  caAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCAtt                                                                                           >  2:30141/1‑149 (MQ=255)
                                                       cACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                      <  8:287397/149‑1 (MQ=255)
                                                        aCGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggcg                                                                                     >  2:358877/1‑149 (MQ=255)
                                                                  gctgctgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGt                                                                           >  5:84561/1‑149 (MQ=255)
                                                                         ctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTgg                                                                    <  6:366856/149‑1 (MQ=255)
                                                                            ccTGCCCAACCAGTGCGCCTGCCATCGAGCCCACTGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACCCTTGGCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                  <  1:214698/148‑1 (MQ=255)
                                                                             cTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTg                                                                 <  8:92162/148‑1 (MQ=255)
                                                                             cTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGc                                                                >  1:66703/1‑149 (MQ=255)
                                                                                cccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGt                                                             <  3:96348/149‑1 (MQ=255)
                                                                                     ccAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTccgcc                                                        >  2:245053/1‑149 (MQ=255)
                                                                                        gTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                  >  4:44717/1‑136 (MQ=255)
                                                                                        gTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                  <  3:44717/136‑1 (MQ=255)
                                                                                         tGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCcgcg                                                     >  6:164923/1‑148 (MQ=255)
                                                                                           cgcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                   <  6:180599/148‑1 (MQ=255)
                                                                                             ccTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAgcg                                                >  6:21386/1‑149 (MQ=255)
                                                                                                    tCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGcc                                          >  8:230202/1‑148 (MQ=255)
                                                                                                           cACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCa                                  <  5:164923/149‑1 (MQ=255)
                                                                                                             cGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAgg                                <  2:66703/149‑1 (MQ=255)
                                                                                                                gAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGAGCAGGca                              >  1:262536/1‑148 (MQ=255)
                                                                                                                 ggTGACCACCGCCCCAACCACAGCGCCGATGATCGGAGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGAAATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCAGCTTAGCGCTTGCCCAGGCGCCGGcacc                            >  8:88142/2‑148 (MQ=255)
                                                                                                                 aGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGcaca                            >  5:189906/1‑149 (MQ=255)
                                                                                                                  gTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                           >  8:152921/1‑117 (MQ=255)
                                                                                                                  gTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                           <  7:152921/117‑1 (MQ=255)
                                                                                                                   tGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACAcc                          >  2:58063/1‑149 (MQ=255)
                                                                                                                     cccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGc                        <  8:195730/149‑1 (MQ=255)
                                                                                                                      ccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCa                       <  7:119051/149‑1 (MQ=255)
                                                                                                                      ccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACCCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCg                         >  1:329726/1‑147 (MQ=255)
                                                                                                                       cACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAg                      >  3:57564/1‑149 (MQ=255)
                                                                                                                          cGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                     <  6:329223/99‑1 (MQ=255)
                                                                                                                          cGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                     >  5:329223/1‑99 (MQ=255)
                                                                                                                           gCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCATGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                  <  3:46438/149‑1 (MQ=255)
                                                                                                                             ccccACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTCGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGt                >  1:302070/1‑149 (MQ=255)
                                                                                                                                 ctgatcgcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACg            >  5:83430/7‑149 (MQ=255)
                                                                                                                                  ccacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACgg           <  5:21386/149‑1 (MQ=255)
                                                                                                                                    acaGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGc         <  1:30141/149‑1 (MQ=255)
                                                                                                                                       gcgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACt      <  6:84561/149‑1 (MQ=255)
                                                                                                                                       gcgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACt      <  2:262536/149‑1 (MQ=255)
                                                                                                                                       gcgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACt      <  1:358877/149‑1 (MQ=255)
                                                                                                                                           cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAg   >  7:138778/1‑148 (MQ=255)
                                                                                                                                           cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCAGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACgg           <  4:295946/140‑1 (MQ=255)
                                                                                                                                           cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCAGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACgg           >  3:295946/1‑140 (MQ=255)
                                                                                                                                           cGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGc  <  6:83430/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGC  >  NC_002947/4561100‑4561387

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: ‑