Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,977,992 C→T 6.8% S121S (TCG→TCA)  PP_2604 ← major facilitator family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,977,9920CT 6.8% 173.7 / 2.1 59S121S (TCG→TCA) PP_2604major facilitator family transporter
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base T (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

GCCAACATTGGTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAACCAGAAATTGATGTTG  >  NC_002947/2977845‑2978139
                                                                                                                                                   |                                                                                                                                                   
gcCAACATTGGTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGcc                                                                                                                                                     >  1:296330/1‑148 (MQ=255)
  cAACATTGGTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAg                                                                                                                                                  <  8:67056/149‑1 (MQ=255)
     cATTGGTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTg                                                                                                                                               <  3:36482/149‑1 (MQ=255)
          gTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCa                                                                                                                                           <  8:267021/148‑1 (MQ=255)
          gTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAg                                                                                                                                          <  5:138139/149‑1 (MQ=255)
                gCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTg                                                                                                                                    <  4:7548/149‑1 (MQ=255)
                gCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTg                                                                                                                                    <  2:99454/149‑1 (MQ=255)
                gCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTg                                                                                                                                    >  2:97783/1‑149 (MQ=255)
                      aacaacCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCAcc                                                                                                                              >  3:228236/1‑149 (MQ=255)
                         aacCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGt                                                                                                                           >  7:237851/1‑149 (MQ=255)
                            cccGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGc                                                                                                                        >  6:72373/1‑149 (MQ=255)
                            cccGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGc                                                                                                                        >  5:183376/1‑149 (MQ=255)
                               gTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAg                                                                                                                     >  1:221907/1‑149 (MQ=255)
                               gTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAg                                                                                                                     >  3:204536/1‑149 (MQ=255)
                               gTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAg                                                                                                                     >  7:336276/1‑149 (MQ=255)
                                      tCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTg                                                                                                                          <  2:340561/137‑1 (MQ=255)
                                      tCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTg                                                                                                                          >  1:340561/1‑137 (MQ=255)
                                        gCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCAcag                                                                                                            >  1:246074/1‑149 (MQ=255)
                                         cccGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCAcagc                                                                                                           <  5:72373/149‑1 (MQ=255)
                                         cccGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCAcag                                                                                                            >  8:35925/1‑148 (MQ=255)
                                                aCCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAg                                                                                                    <  2:296330/149‑1 (MQ=255)
                                                aCCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAg                                                                                                    >  4:94904/1‑149 (MQ=255)
                                                 cccGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGc                                                                                                   <  1:97783/149‑1 (MQ=255)
                                                 cccGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCAGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGTCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGc                                                                                                   >  5:125445/1‑149 (MQ=255)
                                                             cacaGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGcc                                                                                                                       >  6:182379/1‑117 (MQ=255)
                                                             cacaGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGcc                                                                                                                       <  5:182379/117‑1 (MQ=255)
                                                              acaGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGcc                                                                                      <  4:204536/149‑1 (MQ=255)
                                                                    gcACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCTGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGttt                                                                                 >  6:88031/1‑148 (MQ=255)
                                                                    gcACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCTGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGttt                                                                                 >  6:93344/1‑148 (MQ=255)
                                                                    gcACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGttt                                                                                 >  4:157907/1‑148 (MQ=255)
                                                                      aCGGTCTTGACGCAGGCCGGGTAGGTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAATGCCCGTTTGAc                                                                              <  6:183376/149‑1 (MQ=255)
                                                                        ggTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGAGGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTc                                                                            <  8:249271/149‑1 (MQ=255)
                                                                          tCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCg                                                                          <  3:330592/149‑1 (MQ=255)
                                                                          tCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCg                                                                          <  3:117016/149‑1 (MQ=255)
                                                                           cTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGa                                                                         >  1:60648/1‑149 (MQ=255)
                                                                              gACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGAcc                                                                       <  8:237851/148‑1 (MQ=255)
                                                                                  cAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTc                                                                  <  4:228236/149‑1 (MQ=255)
                                                                                     gCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGAcc                                                                       >  2:160524/1‑141 (MQ=255)
                                                                                     gCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGAcc                                                                       <  1:160524/141‑1 (MQ=255)
                                                                                      ccGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATc                                                              <  7:35925/149‑1 (MQ=255)
                                                                                              ttGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTg                                                                                                                                    >  6:91040/1‑71 (MQ=255)
                                                                                                   ggcCGCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAACAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTgca                                                 <  6:128514/149‑1 (MQ=255)
                                                                                                           tCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTgca                                                 <  4:108780/141‑1 (MQ=255)
                                                                                                           tCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTgca                                                 >  3:108780/1‑141 (MQ=255)
                                                                                                               gCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCTGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGa                                     <  5:88031/149‑1 (MQ=255)
                                                                                                               gCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCTGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGa                                     <  5:93344/149‑1 (MQ=255)
                                                                                                                ccGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGTGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGAt                                    <  7:4771/149‑1 (MQ=255)
                                                                                                                  gAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGATGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATgg                                  >  1:59406/1‑149 (MQ=255)
                                                                                                                     gATGAAGCGGCGGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGGAAATAGCCGATGGCa                                >  6:123370/1‑148 (MQ=255)
                                                                                                                      aTGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATa                              <  2:60648/149‑1 (MQ=255)
                                                                                                                        gAAGCGGCAGACCTTGTAACCTGCAGCCGAGGTGGCGAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAg                             <  8:336276/148‑1 (MQ=255)
                                                                                                                          aGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCg                          <  3:94904/149‑1 (MQ=255)
                                                                                                                              gCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCt                      >  6:222275/1‑149 (MQ=255)
                                                                                                                               cAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTg                     <  2:153959/149‑1 (MQ=255)
                                                                                                                                 gACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTgg                    >  5:149072/1‑148 (MQ=255)
                                                                                                                                       gAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAAc                 <  4:307131/145‑1 (MQ=255)
                                                                                                                                       gAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAAc                 >  3:307131/1‑145 (MQ=255)
                                                                                                                                             tGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAACCAGAAATTGa       >  1:207654/1‑149 (MQ=255)
                                                                                                                                               cAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAACCAGAAATCGAt      <  3:157907/148‑1 (MQ=255)
                                                                                                                                               cAGCCGAGGTGGCCAGGCCGGGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAACAACAACGCCCCTTTGACTCCGACCTTCTCGATCAGCCGGCCTTGCAGCAAAAACCCGATGGCATAGCCGACCTGGAACCAGAAATTGa       <  2:246074/147‑1 (MQ=255)
                                                                                                                                                  ccGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAACCAGAAATTGATGTTg  <  2:221907/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCAACATTGGTACCGGCGTTGAACAACCCCGTTGCAATCGCCCGCTCACCCGCCGGGAACCACAGGCGCACGGTCTTGACGCAGGCCGGGTAGTTGGCGGCCTCGGTCAGGCCGAGGATGAAGCGGCAGACCATGAAACCTGCAGCCGAGGTGGCCAGGCCGTGGGCACCCGTGGCCAGGCTCCACAGCAGTACAGCGAAGAAGAAGGCCCGTTTGACTCCGACCTTGTCGATCAGCCGGCCTTGCAGCAAAAAGCCGATGGCATAGCCGACCTGGAACCAGAAATTGATGTTG  >  NC_002947/2977845‑2978139

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