Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,709,233 G→A 6.4% Q110* (CAG→TAG)  PP_4166 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,709,2330GA 6.4% 183.6 / 2.5 62Q110* (CAG→TAG) PP_4166hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/24);  new base A (2/2);  total (36/26)
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.23e-01

TGCGACATGGCGCCTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGC  >  NC_002947/4709086‑4709382
                                                                                                                                                   |                                                                                                                                                     
tgcGACATGGCGCCTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCa                                                                                                                                                     >  2:958799/1‑150 (MQ=255)
           gcCTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACt                                                                                                                                          >  2:701785/1‑150 (MQ=255)
           gcCTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACt                                                                                                                                          >  1:69480/1‑150 (MQ=255)
            ccTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCGCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACtt                                                                                                                                         <  2:696574/150‑1 (MQ=255)
                 tCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGctc                                                                                                                                    >  1:139821/1‑150 (MQ=255)
                   cggaccggaATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCctgcctgcgtcgg                                                                                                                                                 >  1:1113240/1‑130 (MQ=255)
                   cggaccggaATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCctgcctgcgtcgg                                                                                                                                                 <  2:1113240/135‑6 (MQ=255)
                         ggaATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAgcgc                                                                                                                            >  2:57212/1‑150 (MQ=255)
                          gaATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGcc                                                                                                                           >  1:840676/1‑150 (MQ=255)
                           aaTGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCt                                                                                                                          >  2:566047/1‑150 (MQ=255)
                            atgatgGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTc                                                                                                                         <  2:618226/150‑1 (MQ=255)
                               atgGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTcgcg                                                                                                                      <  1:701785/150‑1 (MQ=255)
                                  gCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGc                                                                                                                   >  1:796062/1‑150 (MQ=255)
                                   cccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCtctc                                                                                                                                  <  1:431913/134‑1 (MQ=255)
                                   cccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCtctc                                                                                                                                  >  2:431913/1‑134 (MQ=255)
                                    ccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCa                                                                                                                 >  1:980300/1‑150 (MQ=255)
                                     cATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAg                                                                                                                >  2:76727/1‑150 (MQ=255)
                                      aTTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGa                                                                                                               <  1:958799/150‑1 (MQ=255)
                                          tCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCt                                                                                                           <  2:980300/150‑1 (MQ=255)
                                            ggCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTgg                                                                                                         <  1:566047/150‑1 (MQ=255)
                                                    ggggCAGTTCGGCCATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTa                                                                                                 <  2:243724/150‑1 (MQ=255)
                                                     gggCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAg                                                                                                <  2:1054467/150‑1 (MQ=255)
                                                     gggCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAg                                                                                                >  1:845561/1‑150 (MQ=255)
                                                       gCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTc                                                                                              <  1:632385/150‑1 (MQ=255)
                                                          gTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAg                                                                                           >  2:878289/1‑150 (MQ=255)
                                                           ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTACATGCAGGcctgcc                                                                                                                                         <  1:1187265/103‑5 (MQ=255)
                                                           ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTACATGCAGGcctgcc                                                                                                                                         <  1:1187574/103‑5 (MQ=255)
                                                           ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTACATGCAGGcctgcc                                                                                                                                         >  2:1187265/1‑99 (MQ=255)
                                                           ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTACATGCAGGcctgcc                                                                                                                                         >  2:1187574/1‑99 (MQ=255)
                                                              ggCAATGTCCATAAAGGTCATCGGCGAGGGGGATAAGATCATGCGGAGAGTGCCGAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGc                                                                                       >  1:210728/1‑150 (MQ=255)
                                                                     tCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAggcggc                                                                                >  2:884962/1‑150 (MQ=255)
                                                                     tCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGACTCTCCAGCGCCTCGCGGGCGAGACACTGGTTCGTGTAGTCCAGTTggcggcggc                                                                                >  2:1103341/1‑150 (MQ=255)
                                                                      ccATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAgc                                                                                                                              >  2:789283/1‑103 (MQ=255)
                                                                      ccATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAgc                                                                                                                              <  1:789283/103‑1 (MQ=255)
                                                                       cATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGcc                                                                                                                  <  2:1193023/114‑1 (MQ=255)
                                                                       cATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGcc                                                                                                                  >  1:1193023/1‑114 (MQ=255)
                                                                        aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGgcgcc                                                                             >  1:110524/1‑149 (MQ=255)
                                                                             ggTCATCGGCGAGGGGGATCAGGTCCTGCCGCGCCGTCTTAGCCGCTCAGCGACCGACGCAGGCAGGCCTGCCTGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGcc                                                                        <  1:269800/150‑1 (MQ=255)
                                                                             ggTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGcc                                                                        >  1:94542/1‑150 (MQ=255)
                                                                                 aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAg                                                                                     <  2:659158/133‑1 (MQ=255)
                                                                                 aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAg                                                                                     >  1:659158/1‑133 (MQ=255)
                                                                                 aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGAc                                                                    >  2:625732/1‑150 (MQ=255)
                                                                                      cGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTg                                                               <  2:210728/150‑1 (MQ=255)
                                                                                      cGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTg                                                               <  1:878289/150‑1 (MQ=255)
                                                                                            ggATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTgc                                                         >  1:580953/1‑150 (MQ=255)
                                                                                                  ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAgg                                                   >  2:392663/1‑150 (MQ=255)
                                                                                                    tCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTg                                                 <  2:110524/150‑1 (MQ=255)
                                                                                                         ccGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCa                                                                       >  1:788737/1‑123 (MQ=255)
                                                                                                         ccGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCa                                                                       <  2:788737/123‑1 (MQ=255)
                                                                                                                    aaaCTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCgcaggcag                                 <  1:57212/150‑1 (MQ=255)
                                                                                                                     aaCTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGcaggcagg                                <  2:191519/150‑1 (MQ=255)
                                                                                                                           tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                   >  2:804673/1‑77 (MQ=255)
                                                                                                                           tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                   <  1:804673/77‑1 (MQ=255)
                                                                                                                           tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCgcaggcag                                 >  2:1018574/1‑143 (MQ=255)
                                                                                                                           tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCgcaggcag                                 <  1:1018574/143‑1 (MQ=255)
                                                                                                                             aGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGc                        >  1:809534/1‑150 (MQ=255)
                                                                                                                                 accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCg                    >  1:1162059/1‑150 (MQ=255)
                                                                                                                                    gacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTAc                 >  1:546477/1‑150 (MQ=255)
                                                                                                                                     acGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                     <  1:808317/97‑1 (MQ=255)
                                                                                                                                     acGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                     >  2:808317/1‑97 (MQ=255)
                                                                                                                                        caggcaggCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAg             <  2:796062/150‑1 (MQ=255)
                                                                                                                                           gcaggcCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTc          >  2:116583/1‑150 (MQ=255)
                                                                                                                                                 cTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTc    >  1:797268/1‑150 (MQ=255)
                                                                                                                                                 cTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGG‑TTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCg   >  2:353702/1‑150 (MQ=255)
                                                                                                                                                   gCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGc  >  2:683387/1‑150 (MQ=255)
                                                                                                                                                   gCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGc  <  1:76727/150‑1 (MQ=255)
                                                                                                                                                   gCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGc  >  2:983240/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGCGACATGGCGCCTGGTCCGGACCGGAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGC  >  NC_002947/4709086‑4709382

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