Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,522,499 C→G 5.5% A116A (GCC→GCG idh → isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,522,4990CG 5.5% 238.1 / 3.4 73A116A (GCC→GCGidhisocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base G (2/2);  total (39/34)
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.03e-01

GGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCT  >  NC_002947/4522355‑4522647
                                                                                                                                                |                                                                                                                                                     
ggAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGccgccg                                                                                                                                                   >  1:399175/1‑149 (MQ=255)
 gAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACcgc                                                                                                                                                  <  1:248718/149‑1 (MQ=255)
         ccAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaagaa                                                                                                                                           >  4:229633/1‑148 (MQ=255)
            gCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCGCCGCGGAAGAAAAAg                                                                                                                                       <  2:399175/149‑1 (MQ=255)
              tGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGaa                                                                                                                                     >  4:167043/1‑149 (MQ=255)
              tGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGaa                                                                                                                                     >  1:167120/1‑149 (MQ=255)
                gCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATc                                                                                                                                   <  2:167120/149‑1 (MQ=255)
                       ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        <  7:118254/137‑1 (MQ=255)
                       ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        >  8:118254/1‑137 (MQ=255)
                       ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGcc                                                                                                                            <  5:387345/149‑1 (MQ=255)
                          cTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGc                                                                                                                         >  5:84785/1‑149 (MQ=255)
                               gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaa                                                                                                                                              >  1:200640/1‑123 (MQ=255)
                               gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaa                                                                                                                                              <  2:200640/123‑1 (MQ=255)
                               gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGa                                                                                                                    >  1:153795/1‑149 (MQ=255)
                                ccAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGAc                                                                                                                   >  7:165819/1‑149 (MQ=255)
                                   aCATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        <  7:83194/125‑1 (MQ=255)
                                   aCATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        >  8:83194/1‑125 (MQ=255)
                                       catcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATc                                                                                                             >  2:42081/1‑148 (MQ=255)
                                        atcaAGCTGCCGAACATCAGGGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATc                                                                                                             >  8:279122/1‑147 (MQ=255)
                                         tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          <  8:165819/149‑1 (MQ=255)
                                         tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:195841/1‑149 (MQ=255)
                                         tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:203725/1‑149 (MQ=255)
                                           aaGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAggg                                                                                                        <  1:356844/149‑1 (MQ=255)
                                             gCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCt                                                                                                      >  8:127325/1‑149 (MQ=255)
                                                 ccGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTccgc                                                                                                  >  7:60852/1‑149 (MQ=255)
                                                 ccGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTccgc                                                                                                  >  3:278906/1‑149 (MQ=255)
                                                    aaCATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGt                                                                                               >  8:157300/1‑149 (MQ=255)
                                                       aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:281027/1‑149 (MQ=255)
                                                       aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:237685/1‑149 (MQ=255)
                                                       aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACAACCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:328822/1‑149 (MQ=255)
                                                            cgcCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCgg                                                                                       <  1:162204/149‑1 (MQ=255)
                                                                tCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCTGTGAACCCGGtgct                                                                                   >  8:146833/1‑149 (MQ=255)
                                                                  ggTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:284357/1‑124 (MQ=255)
                                                                  ggTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          <  4:284357/124‑1 (MQ=255)
                                                                      ccGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAAcc                                                                                          <  6:302555/136‑1 (MQ=255)
                                                                      ccGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAAcc                                                                                          >  5:302555/1‑136 (MQ=255)
                                                                       cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCTGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     <  7:110435/140‑1 (MQ=255)
                                                                       cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCTGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     >  8:110435/1‑140 (MQ=255)
                                                                       cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGaa                                                                            <  3:167043/149‑1 (MQ=255)
                                                                        gCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAg                                                                           <  1:42081/149‑1 (MQ=255)
                                                                         cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtgct                                                                                   <  3:372270/140‑1 (MQ=255)
                                                                         cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtgct                                                                                   >  4:372270/1‑140 (MQ=255)
                                                                         cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAgg                                                                          >  4:116808/1‑149 (MQ=255)
                                                                               aaaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTc                                                                   >  2:285219/1‑149 (MQ=255)
                                                                                aaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTTCGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTcc                                                                  >  2:200771/1‑149 (MQ=255)
                                                                                aaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTcc                                                                  <  3:116808/149‑1 (MQ=255)
                                                                                 aGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     >  8:234274/1‑130 (MQ=255)
                                                                                   ccGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGAc                                                               <  7:279122/149‑1 (MQ=255)
                                                                                             ggAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGc                                                                        >  4:2461/1‑130 (MQ=255)
                                                                                             ggAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGc                                                                        <  3:2461/130‑1 (MQ=255)
                                                                                                 cTGCAAGGCAAGGGCTTCAACATCCATGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACAGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGAAACTCCGACAGACGCGCGCAGCt                                                 >  7:151968/1‑149 (MQ=255)
                                                                                                  tGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTg                                                >  4:344697/1‑149 (MQ=255)
                                                                                                    caaggcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTc                                              >  5:246881/1‑149 (MQ=255)
                                                                                                        gcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          <  2:153795/149‑1 (MQ=255)
                                                                                                        gcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          <  6:246881/149‑1 (MQ=255)
                                                                                                          aaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          >  3:179850/1‑147 (MQ=255)
                                                                                                          aaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCaa                                        <  4:203725/149‑1 (MQ=255)
                                                                                                           aggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAg                                       <  2:328822/149‑1 (MQ=255)
                                                                                                            gggCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACAGCGGCGCGCGGCGGTCggagttgc                                      <  8:60852/149‑7 (MQ=255)
                                                                                                            gggCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACAGCGGCGCGCGGCGGTCggagttgc                                      <  4:278906/149‑7 (MQ=255)
                                                                                                                  cAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACg                                >  4:234157/1‑149 (MQ=255)
                                                                                                                      aTCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACTccc                             <  1:285219/149‑1 (MQ=255)
                                                                                                                         ccTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCaa                         <  3:229633/149‑1 (MQ=255)
                                                                                                                          cTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAg                        <  7:127325/149‑1 (MQ=255)
                                                                                                                                 cGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCAccc                   <  2:237685/147‑1 (MQ=255)
                                                                                                                                  gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCa                      <  4:81005/143‑1 (MQ=255)
                                                                                                                                  gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCa                      >  3:81005/1‑143 (MQ=255)
                                                                                                                                  gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcaca                >  3:397112/1‑149 (MQ=255)
                                                                                                                                  gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcac                 >  3:390835/1‑148 (MQ=255)
                                                                                                                                     gacgaGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACaa             <  7:157300/149‑1 (MQ=255)
                                                                                                                                      acgaGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcacac               >  3:280720/1‑146 (MQ=255)
                                                                                                                                          gCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATgg        <  6:84785/149‑1 (MQ=255)
                                                                                                                                              gCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGgtgc    <  4:397112/149‑1 (MQ=255)
                                                                                                                                              gCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGgtgc    <  4:390835/149‑1 (MQ=255)
                                                                                                                                                cACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCt  >  5:261575/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCT  >  NC_002947/4522355‑4522647

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: