Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,452,860 T→C 6.6% L341P (CTG→CCG)  PP_1271 → putative Multidrug efflux MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,452,8600TC 6.6% 221.0 / 4.0 75L341P (CTG→CCG) PP_1271putative Multidrug efflux MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base C (2/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTACCTGCT  >  NC_002947/1452715‑1453004
                                                                                                                                                  |                                                                                                                                                
gTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTCTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGAc                                                                                                                                               >  1:242838/1‑149 (MQ=255)
     ccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCCGTcagcgg                                                                                                                                          >  2:43419/1‑149 (MQ=255)
     ccTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCCGTcagcgg                                                                                                                                          >  3:169509/1‑149 (MQ=255)
       tGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGCTCGCCCTGCCGCTGACGGCGGCa                                                                                                                                        >  7:167947/1‑149 (MQ=255)
         gCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCa                                                                                                                                        <  1:7091/147‑1 (MQ=255)
         gCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCCGTCAGCGGCag                                                                                                                                       <  4:169509/148‑2 (MQ=255)
         gCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCCCAACTGATCGCCCTGCCGCCGTCAGCGGCagg                                                                                                                                      <  1:43419/149‑3 (MQ=255)
           cggcggGCTGGGCGTGCTTTCGGGTGGGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCCGTCAGCGGCagggc                                                                                                                                    <  4:237484/149‑5 (MQ=255)
            ggCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACtgtg                                                                                                                                   >  3:210010/1‑149 (MQ=255)
              cGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTGACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACtgtg                                                                                                                                   >  1:289509/1‑148 (MQ=255)
               gTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCaa                                                                                                                                >  2:311115/1‑149 (MQ=255)
                  cTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCagggc                                                                                                                                    >  4:23976/1‑138 (MQ=255)
                  cTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACAt                                                                                                                             <  8:276885/149‑1 (MQ=255)
                            tttCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTg                                                                                                                   >  8:165063/1‑149 (MQ=255)
                              tCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGa                                                                                                                 >  4:216214/1‑149 (MQ=255)
                                gggTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACt                                                                                                               <  4:210010/149‑1 (MQ=255)
                                      gggCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCt                                                                                                         <  1:139667/149‑1 (MQ=255)
                                        gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTg                                                                                                              <  7:251474/142‑1 (MQ=255)
                                        gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTg                                                                                                              >  8:251474/1‑142 (MQ=255)
                                        gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTgg                                                                                                       <  5:315564/149‑1 (MQ=255)
                                        gcagcaTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTgg                                                                                                       <  7:165063/149‑1 (MQ=255)
                                            caTCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGc                                                                                                   >  3:86395/1‑149 (MQ=255)
                                              tCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGa                                                                                                                                                <  5:121943/102‑1 (MQ=255)
                                              tCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGa                                                                                                                                                >  6:121943/1‑102 (MQ=255)
                                              tCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTg                                                                                                 >  1:56729/1‑149 (MQ=255)
                                               ccccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACtgtg                                                                                                                                   >  4:241400/1‑114 (MQ=255)
                                               ccccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACtgtg                                                                                                                                   <  3:241400/114‑1 (MQ=255)
                                                cccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGt                                                                                                     >  1:162445/1‑143 (MQ=255)
                                                cccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGt                                                                                                     <  2:162445/143‑1 (MQ=255)
                                                 ccTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCAGCGCAGTGGACTGCCGCTGGGTGCTGGGc                                                                                              >  3:227944/1‑149 (MQ=255)
                                                      gCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTgg                                                                                         >  8:16598/1‑149 (MQ=255)
                                                      gCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTgg                                                                                         >  6:299101/1‑149 (MQ=255)
                                                       ccTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGc                                                                                        <  5:223387/149‑1 (MQ=255)
                                                         tACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCTTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCt                                                                                      >  4:19069/1‑149 (MQ=255)
                                                           ccTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTgg                                                                                    <  8:167947/149‑1 (MQ=255)
                                                           ccTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTgg                                                                                    <  7:16598/149‑1 (MQ=255)
                                                            cTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCCGCACAGTGCAACATCCGCGACGTGGACTGCCGCTGGGTGATGGGAATTGGCCTGGc                                                                                   >  4:62335/1‑149 (MQ=255)
                                                             ttGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGt                                                                                                     <  8:58975/130‑1 (MQ=255)
                                                             ttGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGt                                                                                                     >  7:58975/1‑130 (MQ=255)
                                                                 gcAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGc                                                                              >  8:79024/1‑149 (MQ=255)
                                                                       ccAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTACCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGAt                                                                                <  3:186851/141‑1 (MQ=255)
                                                                       ccAGGGCTACCGCCCGGCGCAGACAAGCCCGTTTATGATGCTGGAGGCCATGCCGCAACTGATCGCCCTACAGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGATGGGTGCTGGCCATTGGCCTCGCGAt                                                                                >  4:186851/1‑141 (MQ=255)
                                                                        cAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCgg                                                                        <  3:216214/148‑1 (MQ=255)
                                                                        cAGGGCTACCGCCCAGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCTCGCAGTTGACTGCCGCTGGGTGCTGGGCATTGGCCTGTCGATGCTGGCGgt                                                                       <  3:354318/149‑1 (MQ=255)
                                                                          gggCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCTACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGgtgt                                                                     >  1:276972/1‑149 (MQ=255)
                                                                                 ccGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCtgt                                                                <  1:311115/148‑1 (MQ=255)
                                                                                  cGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCAGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTgtgt                                                              <  4:227944/149‑1 (MQ=255)
                                                                                   gCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGtt                                                             <  2:276972/149‑1 (MQ=255)
                                                                                      cGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGgc                                                          <  2:289509/149‑1 (MQ=255)
                                                                                       ggCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTgg                                                           >  2:45148/1‑147 (MQ=255)
                                                                                             gACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGgcagc                                                       >  5:259027/1‑145 (MQ=255)
                                                                                             gACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGgcagc                                                       <  6:259027/145‑1 (MQ=255)
                                                                                              aCCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCTTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCct                                                  <  3:19069/149‑1 (MQ=255)
                                                                                                             atgCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGt                                    <  5:15262/148‑1 (MQ=255)
                                                                                                              tgCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTgg                                  <  1:45148/149‑1 (MQ=255)
                                                                                                               gCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGa                                 >  2:348588/1‑149 (MQ=255)
                                                                                                                 tggtggCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATc                               >  3:269826/1‑149 (MQ=255)
                                                                                                                       ccATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGc                         >  5:346607/1‑149 (MQ=255)
                                                                                                                        cATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCAGCTGGGTGCTGGGCATTGGACTTGCGATGCTGGCGGTGTCCTGTGTTGGCATCAGACTGCTGACACGAGAGTGGAGACGTgtc                         >  6:250773/1‑146 (MQ=255)
                                                                                                                            ccGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACt                     <  2:56729/148‑1 (MQ=255)
                                                                                                                             cGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGc                                                                                        <  8:57700/80‑1 (MQ=255)
                                                                                                                             cGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGc                                                                                        >  7:57700/1‑80 (MQ=255)
                                                                                                                              gCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGAt                                <  1:227013/135‑1 (MQ=255)
                                                                                                                              gCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGAt                                >  2:227013/1‑135 (MQ=255)
                                                                                                                               cAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTa                 <  4:269826/149‑1 (MQ=255)
                                                                                                                               cAACTGATCGCCCTGCCGATGACGGCGGCACTGTGCACCATCCGCGCCGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGGTGGCAGCAGCCTTCTGACCAGCGAGTGGATCCGTGGCGACTTcaa                 <  5:250773/149‑3 (MQ=255)
                                                                                                                                    gATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGCCTTCTACCCGt            >  4:198849/1‑149 (MQ=255)
                                                                                                                                        gCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGgcagca                                                      >  6:55661/1‑103 (MQ=255)
                                                                                                                                        gCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGgcagca                                                      <  5:55661/103‑1 (MQ=255)
                                                                                                                                        gCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTa        >  5:218070/1‑149 (MQ=255)
                                                                                                                                         cccTGCCGCTGACGGTGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTAc       <  1:330887/149‑1 (MQ=255)
                                                                                                                                          ccTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTAcc      <  1:348588/149‑1 (MQ=255)
                                                                                                                                          ccTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTAcc      >  2:178600/1‑149 (MQ=255)
                                                                                                                                              ccGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTACctgct  <  6:346607/149‑1 (MQ=255)
                                                                                                                                              ccGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTACctgct  >  4:324555/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTGACCCTGGCCGGCGTGCTGGTCGTGCTTTCGGGTGTGGGCAGCATCCCCTCGGCCTACCTTGCGCAGATCCAGGGCTACCGCCCGGCGCAGACCAGCCCGTTGATGATGCTGGTGGCCATGCCGCAACTGATCGCCCTGCCGCTGACGGCGGCACTGTGCAACATCCGCGCAGTGGACTGCCGCTGGGTGCTGGGCATTGGCCTGGCGATGCTGGCGGTGTCCTGTGTTGGCAGCAGCCTGCTGACCAGCGAGTGGATCCGTGGCGACTTCTACCCGTTCTACCTGCT  >  NC_002947/1452715‑1453004

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: