Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,580,599 C→A 6.5% Q943H (CAG→CAT) ‡ treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,580,5990CA 6.5% 239.0 / 6.1 77Q943H (CAG→CAT) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base A (2/3);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

TAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTG  >  NC_002947/4580456‑4580744
                                                                                                                                               |                                                                                                                                                 
tAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgctt                                                                                                                                              >  5:245898/1‑149 (MQ=255)
 aaGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgcttg                                                                                                                                             >  6:4863/1‑149 (MQ=255)
  aGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgcttg                                                                                                                                             <  7:253965/148‑1 (MQ=255)
   gggCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgcttgct                                                                                                                                           <  5:4863/149‑1 (MQ=255)
     gCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCa                                                                                                                                         <  8:331639/149‑1 (MQ=255)
      cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAg                                                                                                                                        <  3:81670/149‑1 (MQ=255)
         tgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                       <  7:73261/147‑1 (MQ=255)
         ggtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                       <  4:230751/146‑1 (MQ=255)
           ttttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtgt                                                                                                                                   <  6:245898/149‑1 (MQ=255)
           ttttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtgg                                                                                                                                   <  3:246179/149‑2 (MQ=255)
             ttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTc                                                                                                                                  <  2:189166/148‑1 (MQ=255)
              tgagcGGCGCTCTTCACGTTGCCGGGAGGGTTCAACTTCAAACTCCTGGAGGTAGGCATCACCCTGCACCACGAGCCCCTCGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGCGCCTGCTTGCTCAGGGTGTCg                                                                                                                                 <  3:228611/144‑1 (MQ=255)
               gCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                               >  8:28536/1‑149 (MQ=255)
                cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACTCGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                               <  1:76888/148‑1 (MQ=255)
                cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtg                                                                                                                              >  2:285041/1‑149 (MQ=255)
                cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                               <  6:105412/148‑1 (MQ=255)
                 ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgt                                                                                                                             >  4:349132/1‑149 (MQ=255)
                  cGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                                                                            >  8:299390/1‑149 (MQ=255)
                   ggCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGc                                                                                                                           <  2:51317/149‑1 (MQ=255)
                       ctctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                       >  2:223730/1‑149 (MQ=255)
                       ctctTCAAGCGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCCTCACCCTTCACCACGAGCACCTGGCCCAGCTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCAGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                       >  5:242780/1‑149 (MQ=255)
                        tctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                      >  5:265437/1‑149 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                      >  3:214567/1‑147 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                      <  5:96804/147‑1 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAggg                                                                                                                                      <  4:102968/131‑1 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAggg                                                                                                                                      >  3:102968/1‑131 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTGCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                      >  5:305837/1‑147 (MQ=255)
                            cAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAgg                                                                                                                  >  6:159686/1‑149 (MQ=255)
                                        cgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGt                                                                                                      <  3:281276/149‑1 (MQ=255)
                                            ttCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgct                                                                                                  <  6:265437/149‑1 (MQ=255)
                                                 cTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAgg                                                                                             <  7:28536/149‑1 (MQ=255)
                                                  ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgt                                                                                            >  3:211244/1‑149 (MQ=255)
                                                      aaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCAT‑AGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtcgtc                                                                                        >  5:243246/1‑148 (MQ=255)
                                                       aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTc                                                                                                             <  6:131246/127‑1 (MQ=255)
                                                       aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTc                                                                                                             >  5:131246/1‑127 (MQ=255)
                                                        aTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCTGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                      <  1:223730/149‑1 (MQ=255)
                                                                ggTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGAc                                                                              >  5:18081/1‑149 (MQ=255)
                                                                 gTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACt                                                                             <  1:285041/149‑1 (MQ=255)
                                                                 gTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACt                                                                             <  3:200881/149‑1 (MQ=255)
                                                                   aGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgctt                                                                                                                                              <  8:30406/82‑1 (MQ=255)
                                                                   aGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgctt                                                                                                                                              >  7:30406/1‑82 (MQ=255)
                                                                    ggCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                           <  5:159686/148‑1 (MQ=255)
                                                                    ggCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                           <  4:214567/148‑1 (MQ=255)
                                                                                caccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCa                                                              <  4:211244/149‑1 (MQ=255)
                                                                                caccaaGAGCTCCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTc                                                                <  6:242780/147‑1 (MQ=255)
                                                                                  ccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                            >  2:65256/1‑149 (MQ=255)
                                                                                  ccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                            <  3:349132/149‑1 (MQ=255)
                                                                                      gAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtgtgt                                                        >  3:358250/1‑149 (MQ=255)
                                                                                            cTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGAt                                                  >  2:220882/1‑149 (MQ=255)
                                                                                               gCCCAGGTGCAGGTCGCCGTGTACACGCAT‑AGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                              <  6:243246/149‑1 (MQ=255)
                                                                                                   aGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGc                                           >  8:18340/1‑149 (MQ=255)
                                                                                                     gTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGc                                           >  4:116614/1‑147 (MQ=255)
                                                                                                           gTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGt                                   >  5:139931/1‑149 (MQ=255)
                                                                                                                    tACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                          <  7:299390/149‑1 (MQ=255)
                                                                                                                      cacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                        <  7:18340/149‑1 (MQ=255)
                                                                                                                       acGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                          >  4:50453/1‑146 (MQ=255)
                                                                                                                       acGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGGCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                          <  3:50453/146‑1 (MQ=255)
                                                                                                                         gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                 >  7:90862/1‑89 (MQ=255)
                                                                                                                         gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                 <  8:90862/89‑1 (MQ=255)
                                                                                                                         gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                <  6:296283/138‑1 (MQ=255)
                                                                                                                         gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                >  5:296283/1‑138 (MQ=255)
                                                                                                                          cATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                    <  1:65256/149‑1 (MQ=255)
                                                                                                                          cATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                    <  2:121568/149‑1 (MQ=255)
                                                                                                                           aTAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcag                   >  4:335168/1‑149 (MQ=255)
                                                                                                                             aaGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGTTGGGCCAGCAGGTCGAGGGCGTGGGtcagtt                 >  5:8214/1‑149 (MQ=255)
                                                                                                                              aGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTTGCGATGCTGGGCCAGCAGGTCGAGGGCTTGGGtcattt                 <  4:73536/148‑1 (MQ=255)
                                                                                                                              aGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCAATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagtt                 >  1:179711/1‑148 (MQ=255)
                                                                                                                                cAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcag              <  6:18081/149‑1 (MQ=255)
                                                                                                                                 aGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcagt             <  7:24583/149‑1 (MQ=255)
                                                                                                                                   gCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGc           >  5:315165/1‑149 (MQ=255)
                                                                                                                                   gCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGc           >  7:206118/1‑149 (MQ=255)
                                                                                                                                    ccACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt          >  3:298613/1‑149 (MQ=255)
                                                                                                                                     cACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCAACCAAGGACTgg                                                                           >  1:326008/1‑83 (MQ=255)
                                                                                                                                     cACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCAACCAAGGACTgg                                                                           <  2:326008/83‑1 (MQ=255)
                                                                                                                                        ccTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgc                                                                                                   >  7:60740/1‑56 (MQ=25)
                                                                                                                                        ccTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgc                                                                                                   <  8:60740/56‑1 (MQ=25)
                                                                                                                                            aGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCtgtg  <  8:177096/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTG  >  NC_002947/4580456‑4580744

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: ‑