Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,485,292 T→C 7.8% W61R (TGG→CGG) ‡ yhdY → amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,485,2920TC 7.8% 215.1 / 6.6 77W61R (TGG→CGG) ‡yhdYamino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (39/32);  new base C (3/3);  total (42/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 = 7.97e-01

GCCACCGGTGAAAACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGC  >  NC_002947/1485144‑1485440
                                                                                                                                                    |                                                                                                                                                    
gccACCGCGGAACACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTCTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGATACTGCCGTGGCCGTTGATCGATGCCAACt                                                                                                                                                      <  2:391332/149‑1 (MQ=255)
 ccACCGGTGAAAACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTg                                                                                                                                                     <  5:82784/149‑1 (MQ=255)
   aCCGGTGAAAACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTCAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTggg                                                                                                                                                   <  6:19700/149‑1 (MQ=255)
        tGAAAACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGc                                                                                                                                              <  8:257650/149‑1 (MQ=255)
          aaaaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcac                                                                                                                                            >  3:313260/1‑149 (MQ=255)
          aaaaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGca                                                                                                                                             <  4:209485/148‑1 (MQ=255)
           aaaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcacc                                                                                                                                           <  4:209834/149‑1 (MQ=255)
           aaaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcac                                                                                                                                            >  4:302199/1‑148 (MQ=255)
            aaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTTCCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcacca                                                                                                                                          <  8:89438/149‑1 (MQ=255)
            aaCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcacca                                                                                                                                          >  3:240656/1‑149 (MQ=255)
             aCCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGcaccac                                                                                                                                         >  6:361744/1‑149 (MQ=255)
               cGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCAcgc                                                                                                                                       >  1:13242/1‑149 (MQ=255)
               cGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACCCAGTTGGcatcgatc                                                                                                                                       <  7:5932/149‑7 (MQ=255)
                   ggCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTTCTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACgcgcgc                                                                                                                                   >  1:354101/1‑149 (MQ=255)
                    gCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGcc                                                                                                                                  <  2:13242/149‑1 (MQ=255)
                          tCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTg                                                                                                                             >  6:28155/1‑148 (MQ=255)
                          tCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACCCAGTTGGcat                                                                                                                                            <  1:21392/133‑2 (MQ=255)
                          tCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACCCAGTTGGcat                                                                                                                                            >  2:21392/1‑132 (MQ=255)
                            gCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTTCTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCa                                                                                                                           <  2:354101/148‑1 (MQ=255)
                                ggATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCaag                                                                                                                      >  3:133326/1‑149 (MQ=255)
                                 gATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCaagg                                                                                                                     <  8:366587/149‑1 (MQ=255)
                                   tGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCaaggaa                                                                                                                   <  5:361744/149‑1 (MQ=255)
                                     cGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCaaggaagg                                                                                                                 <  3:186763/149‑1 (MQ=255)
                                      gTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGgc                                                                                                                >  2:57138/1‑149 (MQ=255)
                                       tGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACATGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGgcg                                                                                                               >  1:347561/1‑149 (MQ=255)
                                        gCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCGCCAAGGAAGgcgc                                                                                                              >  4:58745/1‑149 (MQ=255)
                                           aaTCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCctg                                                                                                           >  1:327667/1‑149 (MQ=255)
                                            aTCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCctgc                                                                                                          >  5:95842/1‑149 (MQ=255)
                                                gTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTggg                                                                                                      >  1:219754/1‑149 (MQ=255)
                                                 ttCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCctgct                                                                                                         <  6:11453/145‑1 (MQ=255)
                                                 ttCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACCCAGTTGGcatc                                                                                                                                           <  8:4876/111‑3 (MQ=255)
                                                 ttCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACCCAGTTGGcatc                                                                                                                                           >  7:4876/1‑109 (MQ=255)
                                                   ctcCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGgtg                                                                                                    <  5:28155/148‑1 (MQ=255)
                                                   ctcCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGGTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGgtgt                                                                                                   >  2:229134/1‑149 (MQ=255)
                                                      cAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCg                                                                                                >  3:291295/1‑149 (MQ=255)
                                                       aGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGt                                                                                               <  4:261523/149‑1 (MQ=255)
                                                            gCTTAACACGCTGTAGACCCTCTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTAACAGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGAGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTCTTCGAgcag                                                                                           >  2:54671/1‑148 (MQ=255)
                                                              ttAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTgg                                                                                                                                               >  3:217618/1‑94 (MQ=255)
                                                              ttAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTgg                                                                                                                                               <  4:217618/94‑1 (MQ=255)
                                                              ttAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGcagcag                                                                                        >  7:28400/1‑149 (MQ=255)
                                                              ttAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGcagcag                                                                                        >  6:346518/1‑149 (MQ=255)
                                                              ttAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGcagcag                                                                                        >  6:264592/1‑149 (MQ=255)
                                                                  cacGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCt                                                                                    <  6:95842/149‑1 (MQ=255)
                                                                  cacGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCt                                                                                    >  6:139684/1‑149 (MQ=255)
                                                                  cacGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCt                                                                                    >  8:135440/1‑149 (MQ=255)
                                                                          tGACCATGTTCGCCATCTTCCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTTGGCGTTGATCGATGCCAACTGGGTTGGCACCACGAGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCa                                                                             <  2:347561/148‑1 (MQ=255)
                                                                           gACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGt                                                                                               >  8:306306/1‑129 (MQ=255)
                                                                           gACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGt                                                                                               <  7:306306/129‑1 (MQ=255)
                                                                             cccTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTc                                                                         <  2:327667/149‑1 (MQ=255)
                                                                              ccTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTAATGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCa                                                                        >  7:196401/1‑149 (MQ=255)
                                                                                    cGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTAc                                                                   <  4:291295/148‑1 (MQ=255)
                                                                                            aCCTGGTCTGGCTGATCGTGCCGCAGTTACTGCAGTGGGCGTTGATCGATGCCAAATTGGTTGGAACCACGCGCGCCTATTGCACCAAGCAACCCGCCTGCTGGGTGTTCGTGCAGAACAGCTTTGGCCAGTTCAAGTACgg                                                                 >  5:313845/1‑142 (MQ=255)
                                                                                             ccTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTAcc                                                         <  5:346518/149‑1 (MQ=255)
                                                                                             ccTGGTGTGGCTGAGCGTGCCGCCGTTACTGCAGGGGGGGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTAc                                                          <  8:360713/148‑1 (MQ=255)
                                                                                                ggTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATCTACGGTTATTACCCGa                                                      >  6:7876/1‑149 (MQ=255)
                                                                                                 gtgtGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGAc                                                     >  5:292502/1‑149 (MQ=255)
                                                                                                 gtgtGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGAc                                                     >  2:202691/1‑149 (MQ=255)
                                                                                                  tgtgGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACt                                                    <  5:264592/149‑1 (MQ=255)
                                                                                                      gCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGaa                                                 <  7:204848/148‑1 (MQ=255)
                                                                                                      gCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGaa                                                 <  4:133326/148‑1 (MQ=255)
                                                                                                      gCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGaa                                                 <  4:347637/148‑1 (MQ=255)
                                                                                                        tGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTg                                              >  2:98527/1‑149 (MQ=255)
                                                                                                          aTCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACt                                                    <  6:125295/141‑1 (MQ=255)
                                                                                                          aTCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACt                                                    >  5:125295/1‑141 (MQ=255)
                                                                                                                  gccgTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCgtg                                    <  1:98527/149‑1 (MQ=255)
                                                                                                                     gTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTTGCCAGTTCTTGTTCGGTTATTACCCGACTGAACTGCGCTGGCGTGTgg                                 <  1:54671/149‑1 (MQ=255)
                                                                                                                                 ggCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCgtgt                     <  4:313260/149‑1 (MQ=255)
                                                                                                                                   cGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTgg                   <  3:302199/149‑1 (MQ=255)
                                                                                                                                    gTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTAccc                                                        <  4:50405/111‑1 (MQ=255)
                                                                                                                                    gTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTAccc                                                        >  3:50405/1‑111 (MQ=255)
                                                                                                                                      tGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCtt                >  5:217364/1‑149 (MQ=255)
                                                                                                                                        aTCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGc              >  2:386627/1‑149 (MQ=255)
                                                                                                                                            gtGCCAACCCAGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTg          >  7:23291/2‑149 (MQ=255)
                                                                                                                                                 aaCTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGTAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCgg     >  2:312231/1‑149 (MQ=255)
                                                                                                                                                 aaCTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGTAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCgg     >  2:307213/1‑149 (MQ=255)
                                                                                                                                                  aCTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGgc    >  7:333700/1‑149 (MQ=255)
                                                                                                                                                   cTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGgcg   >  7:258402/1‑149 (MQ=255)
                                                                                                                                                   cTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCATGCTGGGTGTTCGTGCAGCAGCGCTTTGGACAGTTCATGTACGGTTATTACACGACTGCACTGCGCTCGCGTGTGGACATGACCGTGTGGCTTGCCGTGCTAGgcg   >  1:141186/1‑149 (MQ=255)
                                                                                                                                                    tGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGgcgc  >  1:193681/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCACCGGTGAAAACCGTCGGCGTGCTCGCATGGATGCGTGCCAATCTGTTCTCCAGCTGGCTTAACACGCTGTTGACCCTGTTCGCCATCTACCTGGTGTGGCTGATCGTGCCGCCGTTACTGCAGTGGGCGTTGATCGATGCCAACTGGGTTGGCACCACGCGCGCCGATTGCACCAAGGAAGGCGCCTGCTGGGTGTTCGTGCAGCAGCGCTTTGGCCAGTTCATGTACGGTTATTACCCGACTGAACTGCGCTGGCGTGTGGACCTGACCGTGTGGCTTGCCGTGCTCGGCGC  >  NC_002947/1485144‑1485440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: