Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,459,946 A→T 9.8% L139Q (CTG→CAG)  PP_3075 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,459,9460AT 9.8% 172.7 / 7.3 61L139Q (CTG→CAG) PP_3075transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base T (3/3);  total (31/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 = 9.68e-01

CGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCGGGGA  >  NC_002947/3459812‑3460093
                                                                                                                                      |                                                                                                                                                   
cgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATg                                                                                                                                       <  7:160592/149‑1 (MQ=255)
        ccGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGa                                                                                                                                <  2:315700/148‑1 (MQ=255)
           acgaACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAg                                                                                                                            >  5:147963/1‑149 (MQ=255)
                cTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAgcgcgc                                                                                                                       >  8:222042/1‑149 (MQ=255)
                cTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAgcgcg                                                                                                                        >  6:292497/1‑148 (MQ=255)
                cTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAgcgcg                                                                                                                        <  5:292497/148‑1 (MQ=255)
                   ggCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCt                                                                                                                    <  1:230744/149‑1 (MQ=255)
                     cGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTgg                                                                                                                  >  2:306358/1‑149 (MQ=255)
                           tGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCg                                                                                                            >  7:140250/1‑149 (MQ=255)
                            gTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                           <  6:147963/149‑1 (MQ=255)
                            gTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGCCCATcca                                                                                                           >  7:325772/1‑147 (MQ=255)
                             tATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                           <  2:98308/148‑1 (MQ=255)
                               tttGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGa                                                                                                        <  1:306358/149‑1 (MQ=255)
                                     cTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATggg                                                                                                    >  7:289434/1‑147 (MQ=255)
                                          gggcgCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGgcg                                                                                              >  6:1997/1‑148 (MQ=255)
                                              gcgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGa                                                                                         >  6:336713/1‑149 (MQ=255)
                                                  aGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGaccac                                                                                     <  3:31060/149‑1 (MQ=255)
                                                       cGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACg                                                                                >  5:22341/1‑149 (MQ=255)
                                                           cGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTc                                                                            <  8:289434/149‑1 (MQ=255)
                                                             ggTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGt                                                                          >  4:159317/1‑149 (MQ=255)
                                                                  aGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGc                                                                     <  6:22341/149‑1 (MQ=255)
                                                                  aGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGc                                                                     <  5:336713/149‑1 (MQ=255)
                                                                     ccAGTTCCTGCTGCAGGCGGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgtg                                                                  <  8:286925/149‑1 (MQ=255)
                                                                              gctgcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATg                                                         >  7:210645/1‑149 (MQ=255)
                                                                                 gcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTg                                                                                                                   <  6:232775/88‑1 (MQ=255)
                                                                                 gcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTg                                                                                                                   >  5:232775/1‑88 (MQ=255)
                                                                                    tGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAg                                                                                                 >  8:146615/1‑103 (MQ=255)
                                                                                    tGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAg                                                                                                 <  7:146615/103‑1 (MQ=255)
                                                                                    tGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcag                                                    >  1:339379/1‑148 (MQ=255)
                                                                                      cTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcagcag                                                 >  5:246769/1‑149 (MQ=255)
                                                                                      cTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcagcag                                                 >  2:57982/1‑149 (MQ=255)
                                                                                       tGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcagcag                                                 >  8:122475/1‑148 (MQ=255)
                                                                                              gcgcTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCATGATGggcc                                         <  3:159317/149‑1 (MQ=255)
                                                                                               cgcTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGggccg                                        <  6:246769/149‑1 (MQ=255)
                                                                                                      gcagcGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTg                                 <  7:222042/149‑1 (MQ=255)
                                                                                                        agcGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCt                               >  6:119724/1‑149 (MQ=255)
                                                                                                        agcGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCt                               <  6:25554/149‑1 (MQ=255)
                                                                                                         gcggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGggccgg                                       >  8:155220/1‑140 (MQ=255)
                                                                                                         gcggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGggccgg                                       <  7:155220/140‑1 (MQ=255)
                                                                                                         gcggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCt                               >  4:12947/1‑148 (MQ=255)
                                                                                                          cggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTca                             <  5:119724/149‑1 (MQ=255)
                                                                                                              ttggctttGCGCTGTTCGACGAACTGCGCAGAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTcacca                          >  2:171507/8‑148 (MQ=255)
                                                                                                              gACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGc                                                            >  6:47217/1‑114 (MQ=255)
                                                                                                              gACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGc                                                            <  5:47217/114‑1 (MQ=255)
                                                                                                              gACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTcacca                          >  1:265408/1‑148 (MQ=255)
                                                                                                               aCAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACcaca                        <  2:339379/149‑1 (MQ=255)
                                                                                                                ggctttGCGCTGTTCGACGAACTGCGC‑‑AGCC‑‑TGTCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgc                    <  2:339510/143‑1 (MQ=37)
                                                                                                                   tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgc                    >  3:130481/3‑149 (MQ=255)
                                                                                                                   tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgc                    >  7:99961/3‑149 (MQ=255)
                                                                                                                     tGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgcg                   >  3:74902/1‑148 (MQ=255)
                                                                                                                      gcgcAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCCCCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCgcgc                  >  3:3899/1‑148 (MQ=255)
                                                                                                                       cgcAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCAt                >  5:364589/1‑149 (MQ=255)
                                                                                                                        gctgTTCGACGAACTGCGCAGAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATc               <  1:171507/146‑1 (MQ=255)
                                                                                                                                tcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTg                                 <  6:77786/123‑1 (MQ=255)
                                                                                                                                tcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTg                                 >  5:77786/1‑123 (MQ=255)
                                                                                                                                acgAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGcc       <  8:99961/148‑1 (MQ=255)
                                                                                                                                 cgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCt                               <  2:107705/124‑1 (MQ=255)
                                                                                                                                 cgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGACGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCt                               >  1:107705/1‑124 (MQ=255)
                                                                                                                                  gAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCg      <  8:210645/148‑1 (MQ=255)
                                                                                                                                   aaCAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCggg    <  4:74902/149‑1 (MQ=255)
                                                                                                                                     cAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCGGGGa  <  7:122475/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                   
CGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCGGGGA  >  NC_002947/3459812‑3460093

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: