Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,554,459 G→C 7.5% L25L (CTG→CTC)  PP_3139 → group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,554,4590GC 7.5% 201.7 / 6.9 67L25L (CTG→CTC) PP_3139group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (2/3);  total (32/35)
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.40e-01

GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGG  >  NC_002947/3554311‑3554604
                                                                                                                                                    |                                                                                                                                                 
gCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTg                                                                                                                                                   <  2:61054/149‑1 (MQ=255)
 cccGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTgg                                                                                                                                                  >  7:57390/1‑149 (MQ=255)
  ccGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTg                                                                                                                                                   >  1:287948/1‑147 (MQ=255)
       gTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAgcgc                                                                                                                                             <  8:178660/148‑1 (MQ=255)
        tCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGcc                                                                                                                                            >  3:84638/1‑148 (MQ=255)
          gACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAg                                                                                                                                         <  2:287948/149‑1 (MQ=255)
          gACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCag                                                                                                                                          <  5:27519/148‑2 (MQ=255)
           aCAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAgg                                                                                                                                        <  2:113814/149‑1 (MQ=255)
            cAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAggg                                                                                                                                       <  8:101118/149‑1 (MQ=255)
            cAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAggg                                                                                                                                       <  3:203400/149‑1 (MQ=255)
                   aagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCTGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAgg                                                                                                                                <  3:3326/149‑1 (MQ=255)
                   aagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg                                                                                                                                 <  1:38035/148‑1 (MQ=255)
                    agttaagGTTTCCAGGGAAGGTTTGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGt                                                                                                                               <  8:57390/149‑1 (MQ=255)
                     gtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTg                                                                                                                              <  6:213330/149‑1 (MQ=255)
                     gtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTg                                                                                                                              >  7:23014/1‑149 (MQ=255)
                       gaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGc                                                                                                                             >  8:156803/1‑148 (MQ=255)
                          ggTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACg                                                                                                                         <  7:156803/149‑1 (MQ=255)
                              gCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGta                                                                                                                       >  7:91320/1‑147 (MQ=255)
                                         ttGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgg                                                                                                          <  8:141509/149‑1 (MQ=255)
                                            ggcatcaAGCGGGTTTGGATCTGGATTCGTATGCTCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCCCAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                       <  8:308804/147‑1 (MQ=255)
                                             tcatcaAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtgg                                                                                                      <  3:131708/149‑1 (MQ=255)
                                               atcaAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGa                                                                                                    >  2:360390/1‑149 (MQ=255)
                                                    gCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTAGCGAAATACTGGCGCTGGAGCGCCAAGGGCTTGAGGTGCAACGTATAGACATGCGTGGGTGGGAAGcc                                                                                                >  4:244318/1‑148 (MQ=255)
                                                       ggTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAg                                                                                            <  3:162704/149‑1 (MQ=255)
                                                            ggATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCa                                                                                        <  1:284732/148‑1 (MQ=255)
                                                             gATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAg                                                                                       <  4:166753/148‑1 (MQ=255)
                                                                  ggATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCTACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGcc                                                                                 >  3:344217/1‑149 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGcc                                                                                                                  <  6:176740/115‑1 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGcc                                                                                                                  >  5:176740/1‑115 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGAGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAAc                                                                                                                          <  8:166970/107‑1 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTCATCAAACAGTACCCGAAAGTCAGCCACACTATCATCCGTCGCGAAATACTGGCGCAGGAGCGCCATGGCGGTGAGGTGCTGCGTATAGCCCTGCGTGGGTGGGACGCCTCAGTGCAGGATGcc                                                                                 >  2:311353/1‑148 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTAATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGTGGTTGAGGTGCAAc                                                                                                                          >  7:166970/1‑107 (MQ=255)
                                                                           tGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGc                                                                         >  3:76819/1‑148 (MQ=255)
                                                                             cgcATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAg                                                                      <  2:298986/149‑1 (MQ=255)
                                                                                aTTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCg                                                                    >  8:147197/1‑148 (MQ=255)
                                                                                    cTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACg                                                                <  1:360390/148‑1 (MQ=255)
                                                                                        cTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGc                                                            >  2:110833/1‑148 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCGGAGGACGCACGCGACCGGGGCat                                                          >  4:310582/1‑148 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                             <  4:141570/146‑1 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                             >  3:141570/1‑146 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGc                                                            >  2:64548/1‑147 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGc                                                            >  3:93785/1‑147 (MQ=255)
                                                                                           catcaATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGa                                                        >  5:68021/1‑149 (MQ=255)
                                                                                             tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAcc                                                      >  1:275341/1‑149 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCTACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGc                                                   <  4:344217/149‑1 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCg                                                                                               >  2:138323/1‑105 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCg                                                                                               <  1:138323/105‑1 (MQ=255)
                                                                                                    gTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGACCGGGGCAAGACCCGCTACg                                               >  4:82884/1‑149 (MQ=255)
                                                                                                     tACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGt                                              <  3:244318/149‑1 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGc                                                                                  <  8:226461/110‑1 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGc                                                                                  >  7:226461/1‑110 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAtt                                           >  6:158753/1‑149 (MQ=255)
                                                                                                         cGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTg                                          <  2:275341/149‑1 (MQ=255)
                                                                                                                  gCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTAc                                                <  7:129880/134‑1 (MQ=255)
                                                                                                                  gCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTAc                                                >  8:129880/1‑134 (MQ=255)
                                                                                                                   ccACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                                >  2:357082/1‑149 (MQ=255)
                                                                                                                    cacaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGAAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCt                               >  1:287013/1‑149 (MQ=255)
                                                                                                                            cATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAgggg                        <  4:76819/148‑1 (MQ=255)
                                                                                                                                cGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgtt                   <  6:68021/149‑1 (MQ=255)
                                                                                                                                 gTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGttgttg                  <  7:147197/149‑1 (MQ=255)
                                                                                                                                    gcgATATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGTTGGGCCGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCa               <  1:147210/149‑1 (MQ=255)
                                                                                                                                    gcgAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgg                                                                                                          <  7:291722/58‑1 (MQ=255)
                                                                                                                                    gcgAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgg                                                                                                          >  8:291722/1‑58 (MQ=255)
                                                                                                                                          ccttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACg         >  5:292699/4‑149 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                                <  8:11611/123‑1 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                                >  7:11611/3‑125 (MQ=255)
                                                                                                                                           aCTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGt        >  7:280870/1‑149 (MQ=255)
                                                                                                                                           aCTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGt        >  7:194667/1‑149 (MQ=255)
                                                                                                                                            cTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTg       <  4:84638/149‑1 (MQ=255)
                                                                                                                                                 gcTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAgg  >  1:166798/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGG  >  NC_002947/3554311‑3554604

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