Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,122,074 G→T 5.9% P26T (CCC→ACC)  PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,122,0740GT 5.9% 255.9 / 4.4 84P26T (CCC→ACC) PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (32/47);  new base T (2/3);  total (34/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

GCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTAGAGCGCCTC  >  NC_002947/6121926‑6122221
                                                                                                                                                     |                                                                                                                                                   
gCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGTGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggg                                                                                                                                                     <  3:264769/149‑1 (MQ=255)
 cGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCgggg                                                                                                                                                    <  4:360392/149‑1 (MQ=255)
   gCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggggg                                                                                                                                                   >  1:194477/1‑148 (MQ=255)
   gCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggggg                                                                                                                                                   >  2:93808/1‑148 (MQ=255)
   gCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGc                                                                                                                                                  <  1:155765/149‑1 (MQ=255)
      aCCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACTGGGGGCAAt                                                                                                                                               <  5:345375/149‑1 (MQ=255)
        cAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGc                                                                                                                                             >  8:73069/1‑149 (MQ=255)
             gTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgtc                                                                                                                                         >  8:35352/1‑148 (MQ=255)
             gTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgtc                                                                                                                                         <  7:73069/148‑1 (MQ=255)
             gTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgtc                                                                                                                                         <  7:35352/148‑1 (MQ=255)
                 cgcCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCgg                                                                                                                                    <  5:173975/149‑1 (MQ=255)
                   ccTGGTTAACCGGCGGGGTAGGCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGTCCTTGCTCTCTGGGCATAGGCAGTCTTCGGTCACGTGCTCCCAGCGGGTACCGGGGGCAATGCCGTGGTCGGGc                                                                                                                                  <  8:284474/149‑1 (MQ=255)
                   ccTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGTGGCAATGCCGTCGTCGGGc                                                                                                                                  >  7:336769/1‑149 (MQ=255)
                   ccTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGTGGCAATGCCGTCGTCGGGc                                                                                                                                  >  7:331831/1‑149 (MQ=255)
                      ggATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCaa                                                                                                                               >  4:349069/1‑149 (MQ=255)
                           aacGGCGCGGTAGCCACTGCTGCGATCGGGGTGCCTAGTTCGCTCATTTCAAAGGCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                                                                                          <  6:223123/147‑1 (MQ=255)
                           aCCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                                                                                          <  4:321752/149‑1 (MQ=255)
                            ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCTGGGGCAATGCCGTCGTCGGGCCAACCTTTg                                                                                                                         <  5:22061/149‑1 (MQ=255)
                            ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTg                                                                                                                         <  4:338595/149‑1 (MQ=255)
                             cGGCGCGGTAGCCAATGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAGCTTGCACTCAGGGCATAGCTAGTCTTCTTGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTgg                                                                                                                        <  7:244611/149‑1 (MQ=255)
                                cgcgGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCtt                                                                                                                     <  3:339224/149‑1 (MQ=255)
                                 gcgGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtc                                                                                                                    <  8:208187/149‑1 (MQ=255)
                                       gCCACTTCTGCGATCTGGTTGCTGAGTTCTATCATTTCAAAGTCGGCTTAGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCTGTGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtcg                                                                                                                   <  6:307512/144‑1 (MQ=255)
                                       gCCAATGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCAGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGTGCAATGCCGTAGTCTGTCCAACCTTTGGCTtcg                                                                                                                   >  5:307512/1‑144 (MQ=255)
                                        ccACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCTGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCTTCGTCGGGCCAACCTTTGGCTTCGTCGTAg                                                                                                             <  8:169737/149‑1 (MQ=255)
                                        ccACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAg                                                                                                             <  5:327981/149‑1 (MQ=255)
                                        ccACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAg                                                                                                             <  1:303487/149‑1 (MQ=255)
                                         cACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCACCCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAg                                                                                                             <  1:93808/149‑1 (MQ=255)
                                          aCTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGTACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtcgtcg                                                                                                                >  7:245398/1‑144 (MQ=255)
                                              ctgcGATCGTGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGGCTTCGGGCACGTCCTCCCAGCGTGTACCGTGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGa                                                                                                       <  5:236337/149‑1 (MQ=255)
                                                  gATCGGGGGGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACcac                                                                                                   <  2:259621/149‑1 (MQ=255)
                                                       gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACTTCCTCCCAGCGGGTACAGGGGTCAATGACGTCGTCTGGACAACATTTGGCTTAGTAGTAGATCAGACCAAAgat                                                                                               >  8:99044/1‑148 (MQ=255)
                                                          tGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGatgatg                                                                                           >  5:151619/1‑149 (MQ=255)
                                                             tGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGTGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCag                                                                                        <  8:336769/149‑1 (MQ=255)
                                                             tGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGTGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCag                                                                                        <  8:331831/149‑1 (MQ=255)
                                                             tGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCt                                                                                                                      <  5:316963/119‑1 (MQ=255)
                                                             tGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCt                                                                                                                      >  6:316963/1‑119 (MQ=255)
                                                                 ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatgat                                                                                            >  4:242780/1‑141 (MQ=255)
                                                                 ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatgat                                                                                            <  3:242780/141‑1 (MQ=255)
                                                                 ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGcc                                                                                    <  2:106481/149‑1 (MQ=255)
                                                                 ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCCTCATCGGGCCAACCTTTGGCATCGTCGTAGATCAGACCACAgatcaagc                                                                                          >  8:280650/1‑138 (MQ=255)
                                                                 ttCGATCATTTCAAAGTCGGCTTAGCCGACCTTGCAGTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGc                                                                                          <  7:280650/143‑1 (MQ=255)
                                                                    gATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAgg                                                                                 >  1:141992/1‑149 (MQ=255)
                                                                    gATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAgg                                                                                 >  8:283362/1‑149 (MQ=255)
                                                                    gATCATTTCAAAGTCGGCTTTGCCGAACTTGCACTAAGGGCATAGAAAGTCTTCGGGAAAGTCATCCCAGAGGGTCCCGGGGGCAATGACGTCGTAGGGCCAACATTTGTCTTCTTAGTAGATAAGACCACAGATGATGAAGAGCCAgt                                                                                 >  2:284898/1‑148 (MQ=255)
                                                                        aTTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGtttt                                                                             <  4:321143/149‑1 (MQ=255)
                                                                         tttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTc                                                                            <  5:169512/149‑1 (MQ=255)
                                                                           tCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCa                                                                           <  7:203979/148‑1 (MQ=255)
                                                                           gcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCATCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCAt                                                                          >  3:153283/2‑149 (MQ=255)
                                                                             aaaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAga                                                                                                <  6:337458/125‑1 (MQ=255)
                                                                             aaaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAga                                                                                                >  5:337458/1‑125 (MQ=255)
                                                                                   ggCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgctgc                                                                  <  7:75327/149‑1 (MQ=255)
                                                                                    gCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgctgc                                                                  >  4:307848/1‑148 (MQ=255)
                                                                                     cTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcct                                                                <  2:194477/149‑1 (MQ=255)
                                                                                          ccGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCa                                                           >  8:288872/1‑149 (MQ=255)
                                                                                             aCCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGTTACCGCGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTc                                                        <  5:175607/149‑1 (MQ=255)
                                                                                                ttGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGcc                                                     >  7:169043/1‑149 (MQ=255)
                                                                                                  gCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGTGCAATGCCGTCGTCGGGCCAACCGTTGGCTTCGTCGGAGATCAGACCACAGATGATGCAGAGCCAGTTTTTCATTATGATGCCTCCTCACTAGCCCt                                                   >  1:308276/1‑149 (MQ=255)
                                                                                                  gCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGcc                                                     <  3:307848/147‑1 (MQ=255)
                                                                                                    aCTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCtt                                                  >  5:215247/1‑148 (MQ=255)
                                                                                                    aCTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCtt                                                  >  7:196889/1‑148 (MQ=255)
                                                                                                     cTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGTGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTa                                                                        >  8:26427/1‑125 (MQ=255)
                                                                                                     cTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTa                                                                        <  7:26427/125‑1 (MQ=255)
                                                                                                      tCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGa                                                >  3:13910/1‑148 (MQ=255)
                                                                                                           gCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCa                                          >  5:7833/1‑149 (MQ=255)
                                                                                                                  cAGTCTTAGTGCACTTACTCCCAGCGGGTACGGGGGGCAATGCCGTCGTCGGTCCAACCTTTGGCTTCGTCGTAGATCAGACAACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGTCCAg                                   <  7:99044/149‑1 (MQ=255)
                                                                                                                    gTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGtt                                >  4:329424/1‑149 (MQ=255)
                                                                                                                     tCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGcc                                                     >  3:127859/1‑128 (MQ=255)
                                                                                                                     tCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGcc                                                     <  4:127859/128‑1 (MQ=255)
                                                                                                                         caggCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTc                               <  3:244459/144‑1 (MQ=255)
                                                                                                                         cGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCtt                             <  7:283362/148‑1 (MQ=255)
                                                                                                                          gggCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATTAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcctcc                                                              <  1:78279/114‑1 (MQ=255)
                                                                                                                          gggCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTgg                           <  6:215247/149‑1 (MQ=255)
                                                                                                                          gggCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTgg                           <  8:196889/149‑1 (MQ=255)
                                                                                                                              aCGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCt                       <  3:296713/149‑1 (MQ=255)
                                                                                                                              aCGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCt                       >  5:338544/1‑149 (MQ=255)
                                                                                                                                 tcctccCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGAc                    <  7:288872/149‑1 (MQ=255)
                                                                                                                                      ccAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAg                <  6:7833/148‑1 (MQ=255)
                                                                                                                                      ccAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGATCTTGGTCCTGAAGTAGc               >  5:94263/1‑149 (MQ=255)
                                                                                                                                       cAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCt              >  3:234459/1‑149 (MQ=255)
                                                                                                                                       cAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCt              >  3:111284/1‑149 (MQ=255)
                                                                                                                                           cGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTa           <  4:13910/149‑1 (MQ=255)
                                                                                                                                           cGGGTACCGGGGGCAATGCCGGCGGCGGGCCAACCTTTGGCGTCGGCCTAGATCAGACCACACCTGCAGCAGACCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTa           <  8:245398/149‑1 (MQ=255)
                                                                                                                                             ggTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTaga         <  3:349069/149‑1 (MQ=255)
                                                                                                                                              gTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTagag        >  7:140664/1‑149 (MQ=255)
                                                                                                                                              gTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTagag        <  6:151619/149‑1 (MQ=255)
                                                                                                                                               tACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTAGAgc       >  7:318301/1‑149 (MQ=255)
                                                                                                                                                  cGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTAGAGCGcc    <  8:140664/149‑1 (MQ=255)
                                                                                                                                                    ggggCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATTCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTAGAGCGcctc  <  6:338544/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTCACTCGCCCTTGATACGCACGGCCAGGTTCTTGGTCCTGACGTAGCTGTAGAGCGCCTC  >  NC_002947/6121926‑6122221

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: