Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,599,859 C→T 5.6% R442H (CGC→CAC)  bglX ← beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,599,8590CT 5.6% 217.1 / 2.2 71R442H (CGC→CAC) bglXbeta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATCAGCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGTTTTCCAG  >  NC_002947/1599712‑1600007
                                                                                                                                                   |                                                                                                                                                    
tCGATCAGCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCAcggcgcggcg                                                                                                                                                     >  4:138666/1‑149 (MQ=255)
tCGATCAGCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCAcggcgcggcg                                                                                                                                                     >  1:151002/1‑148 (MQ=255)
   aTCAGCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACggcgcggcgc                                                                                                                                                    >  6:263671/1‑147 (MQ=255)
       gCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCt                                                                                                                                              >  5:257329/1‑149 (MQ=255)
       gCACCGCAGGCGGGCGTGGGTCGTCGACGAGTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAGGATCGCCTTGTCGGCCGTGAGGTTGGAGCCTGTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGGAGGCCt                                                                                                                                              >  4:336751/1‑149 (MQ=255)
        cACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCtt                                                                                                                                             >  4:257123/1‑149 (MQ=255)
         aCCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGAAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTc                                                                                                                                            <  5:266868/149‑1 (MQ=255)
           cGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTcgc                                                                                                                                          <  6:4543/149‑1 (MQ=255)
             cAGGCTGGCTTTGGTATTCGACGATTTACGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATAGCCTTGTCGACCGTGACGTTGTAGCCTTTTTAGTAGACCAGCTTGGCATTTCCTTCCACTGCGCGTCTCAAGCCTTcgcg                                                                                                                                         >  2:313443/1‑148 (MQ=255)
             cAGGCGGGCGTGGGTCGTCGACAATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGACGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCgcgc                                                                                                                                        <  8:67819/149‑1 (MQ=255)
               ggcgggcgTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCAc                                                                                                                                      <  3:121995/149‑1 (MQ=255)
                gcgggcgTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACg                                                                                                                                     <  5:87621/149‑1 (MQ=255)
                gcgggcgTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACg                                                                                                                                     <  5:93810/149‑1 (MQ=255)
                gcgggcgTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACg                                                                                                                                     <  1:45666/149‑1 (MQ=255)
                         gggcgtcgACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGa                                                                                                                             <  4:126017/145‑1 (MQ=255)
                          gtcgtcGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGt                                                                                                                           <  7:138055/149‑1 (MQ=255)
                          gtcgtcGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGt                                                                                                                           <  4:149218/149‑1 (MQ=255)
                           tcgtcgACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTg                                                                                                                          <  1:260620/149‑1 (MQ=255)
                                gacgaTTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACgg                                                                                                                     <  4:318718/149‑1 (MQ=255)
                                gacgaTGTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCATGTCGCCCGGGACGTTGGAGCCTTTGGCGTAGACCAGCATGGCGATGCCTGCCACGGCGCGCCGCACGCCGTCGCGCACGCTCACCGCGAGGCCgg                                                                                                                     >  8:301836/1‑149 (MQ=255)
                                     tttCCGGTGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTa                                                                                                                <  1:313443/149‑1 (MQ=255)
                                      ttCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATTGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTa                                                                                                                <  1:281097/148‑1 (MQ=255)
                                      ttCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACg                                                                                                                                     >  1:343197/1‑127 (MQ=255)
                                      ttCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACg                                                                                                                                     <  2:343197/127‑1 (MQ=255)
                                          ggggCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTc                                                                                                           <  5:18289/149‑1 (MQ=255)
                                           gggCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTcg                                                                                                          >  8:98078/1‑149 (MQ=255)
                                                         aGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGc                                                                                             <  5:166858/148‑1 (MQ=255)
                                                          ggAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGccc                                                                                           <  3:336751/149‑1 (MQ=255)
                                                           gAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGccc                                                                                           >  4:42126/1‑148 (MQ=255)
                                                           gAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCca                                                                                          <  2:239275/149‑1 (MQ=255)
                                                              gTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCcatc                                                                                        >  7:296280/1‑148 (MQ=255)
                                                              gTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCGTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCcatc                                                                                        <  8:296280/148‑1 (MQ=255)
                                                                 cAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATg                                                                                     <  6:257329/148‑1 (MQ=255)
                                                                  aGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGt                                                                                    <  3:42126/148‑1 (MQ=255)
                                                                  aGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGt                                                                                    <  2:8716/148‑1 (MQ=255)
                                                                  aGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTc                                                                                   >  2:217499/1‑149 (MQ=255)
                                                                     tAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGt                                                                                                                   <  6:10263/114‑1 (MQ=255)
                                                                     tAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGt                                                                                                                   >  5:10263/1‑114 (MQ=255)
                                                                          gAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGa                                                                            <  2:228321/148‑1 (MQ=255)
                                                                           aaGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGTGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCCGCTGGCGATCATGTGGATCGGGGc                                                                          >  4:191727/1‑149 (MQ=255)
                                                                           aaGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAg                                                                           >  2:255953/1‑148 (MQ=255)
                                                                           aaGATCGCCTTGTCGCCCGGGACGGTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGTGCAAGCCGTCGCGCACGGGCCCCGAGTGGACGGGGTTCCCGTCGGCGGCCCCGCGGCCCATCATGTCGCTCggggg                                                                          >  4:196672/1‑146 (MQ=255)
                                                                              aTCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATc                                                                       >  6:327507/1‑149 (MQ=255)
                                                                                      gTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAg                                                                 <  7:98078/147‑1 (MQ=255)
                                                                                       tCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGTGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCg                                                               <  3:196672/148‑1 (MQ=255)
                                                                                       tCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGTGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCg                                                               <  3:191727/148‑1 (MQ=255)
                                                                                         gCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGTCGGCCCAGCACCCCATCATGTCGAACCGAGCAAAGGCCAGCCGGc                                                            >  2:263356/1‑149 (MQ=255)
                                                                                            cGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGa                                                         >  3:94758/1‑149 (MQ=255)
                                                                                             gTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCCGCTGGCCGat                                                        <  2:151002/149‑2 (MQ=255)
                                                                                                   ttGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAgg                                                   >  7:325772/1‑148 (MQ=255)
                                                                                                   ttGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAgg                                                   >  8:149304/1‑148 (MQ=255)
                                                                                                       aGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATGGGCGCATAGGCCAGCGGGCCGACCAAGGCAATg                                              >  4:44264/1‑149 (MQ=255)
                                                                                                       aGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATg                                              <  3:138666/149‑1 (MQ=255)
                                                                                                           tttGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGc                                          >  5:340475/1‑149 (MQ=255)
                                                                                                              ggCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGGCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCg                                        <  5:263671/148‑1 (MQ=255)
                                                                                                                 gTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGTGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCt                                    >  1:108406/1‑149 (MQ=255)
                                                                                                                 gTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCCTCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCt                                    >  7:287398/1‑149 (MQ=255)
                                                                                                                 gTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCt                                    >  2:275114/1‑149 (MQ=255)
                                                                                                                 gTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATACGCCAGCGGGCCGACCAACGCAATGGTGCCgggcg                                    >  8:305772/1‑146 (MQ=255)
                                                                                                                    gACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCtttt                                 >  3:319113/1‑149 (MQ=255)
                                                                                                                    gACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCtttt                                 <  3:91553/149‑1 (MQ=255)
                                                                                                                     aCCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGAGCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCttttt                                <  2:219087/149‑1 (MQ=255)
                                                                                                                     aCCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCCGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCttt                                  <  3:158080/147‑1 (MQ=255)
                                                                                                                     aCCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCCGCACAAAATGTCGATCGGAGCAACCGCCAGCCGGCCGACCAACGCAAAGGTGGCCGCCgtt                                  >  4:158080/1‑147 (MQ=255)
                                                                                                                       cAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCCGCATGTCGAGCGGAGCATCGGCCCGCGGGCCGCCCCCGGCAATGGTGCCGGCCGTTTTcc                              >  6:293812/1‑148 (MQ=255)
                                                                                                                          cTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTGAGCg                           >  5:212003/1‑149 (MQ=255)
                                                                                                                          cTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCg                           >  3:202065/1‑149 (MQ=255)
                                                                                                                                   gccttCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCt                  >  8:6369/1‑149 (MQ=255)
                                                                                                                                    ccttCCACGGCGCGGCGCATGCCTTCGGGCACGGTGACGGCGTGGACGGGTTTACGGTCGGCGGCCCAGCTGCCCTTGATGTCGATCGGAGCAGCGGCCTGCGGGGCGTCCATGGCTTTGGTGCGGGCCTTTTTCAGCGGCAGGGTCTg                 <  7:301836/149‑1 (MQ=255)
                                                                                                                                    ccttCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGCGCATCGGCCAGCGGGCCGCCCAAGGCAATGGTGCCGGCCGTTTTCCGCGGCAGGGTcgg                 >  8:194919/1‑147 (MQ=255)
                                                                                                                                     cttCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTgg                <  1:255953/149‑1 (MQ=255)
                                                                                                                                      ttCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGt               >  3:29942/1‑149 (MQ=255)
                                                                                                                                      ttCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCAGGCCATCATGTCGATCGGAGCATCCGCCCGCGGGCCGACCAAGGCAATGGTGCCGGCCGTTTTAAGCGGCAGGGTGTggg               >  8:126290/1‑148 (MQ=255)
                                                                                                                                         cACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGCTCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGc            >  4:250556/1‑149 (MQ=255)
                                                                                                                                             gcgcggcgcATGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCa                              <  7:355441/127‑1 (MQ=255)
                                                                                                                                             gcgcggcgcAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCTGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGtt        <  3:129275/149‑1 (MQ=255)
                                                                                                                                             gcgcggcgcAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCGGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGtt        >  4:129275/1‑149 (MQ=255)
                                                                                                                                                cggcgcAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGTTTTc     <  1:217499/149‑1 (MQ=255)
                                                                                                                                                  gcgcAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGGTTTc     >  7:348146/1‑147 (MQ=255)
                                                                                                                                                   cgcAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGTTTTCcag  <  6:340475/149‑1 (MQ=255)
                                                                                                                                                   cgcAAGCCTTCGCGCACGGGCACCGAGTGGACGGGATTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGTTTTCcac  <  7:179143/149‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGATCAGCACCGCAGGCGGGCGTGGGTCGTCGACGATTTCCGGGGCATCGAAGTTGAGGAAGTTCAGGTAATCGAAGATCGCCTTGTCGCCCGTGACGTTGGAGCCTTTGGCGTAGACCAGCTTGGCCTTGCCTTCCACGGCGCGGCGCAAGCCTTCGCGCACGGTCACCGAGTGGACGGGTTTACCGTCGGCGGCCCAGCTGCCCATCATGTCGATCGGAGCATCGGCCAGCGGGCCGACCAAGGCAATGGTGCCGGCCTTTTTCAGCGGCAGGGTCTGGTTGCGGTTTTCCAG  >  NC_002947/1599712‑1600007

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