Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,322,861 A→C 5.6% I332S (ATC→AGC)  PP_1155 ← YaiC/YhcK/ArdA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,322,8610AC 5.6% 228.5 / 3.4 71I332S (ATC→AGC) PP_1155YaiC/YhcK/ArdA family protein
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GGCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCA  >  NC_002947/1322716‑1323008
                                                                                                                                                  |                                                                                                                                                   
ggCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGAt                                                                                                                                                    <  4:62336/147‑1 (MQ=255)
   ggTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTa                                                                                                                                                 >  1:53901/1‑147 (MQ=255)
     tCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAgc                                                                                                                                             >  3:11825/1‑149 (MQ=255)
       cTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAgcgc                                                                                                                                           >  7:254757/1‑149 (MQ=255)
       cTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACATGATCGCGGCCAGCGTGCTTGGCGATGTACAgcgc                                                                                                                                           >  2:127568/1‑149 (MQ=255)
         ggCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGcc                                                                                                                                          <  1:276720/148‑1 (MQ=255)
         ggCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGcc                                                                                                                                          >  2:276720/1‑148 (MQ=255)
          gCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGcc                                                                                                                                          <  7:193169/147‑1 (MQ=255)
            cATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTg                                                                                                                                       >  8:113022/1‑148 (MQ=255)
             aTGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCtt                                                                                                                                        >  5:16523/1‑146 (MQ=255)
             aTGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCtt                                                                                                                                        <  6:16523/146‑1 (MQ=255)
             aTGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTc                                                                                                                                     >  2:200709/1‑149 (MQ=255)
              tGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCt                                                                                                                                    <  2:317086/149‑1 (MQ=255)
               gcgcGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaagcacg                                                                                                                                   >  2:120980/1‑142 (MQ=255)
                ctcgAGCGTCTCCTTCAGGCAACCTCCCGCCAGTCATTCAGGAAGGAGGTTCAATGAGGATAGGCCAGCCTGAGCGTGGCAGCTTCTCCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGc                                                                                                                                  <  6:80575/147‑1 (MQ=255)
                cgcgAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaagcacg                                                                                                                                   <  1:120980/148‑8 (MQ=255)
                  cgAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCt                                                                                                                                <  4:142352/149‑1 (MQ=255)
                       ggCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGa                                                                                                                                                     >  6:119090/1‑123 (MQ=255)
                       ggCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGa                                                                                                                                                     <  5:119090/123‑1 (MQ=255)
                            cTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGc                                                                                                                        <  7:113022/147‑1 (MQ=255)
                             ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAg                                                                                                                     >  4:362995/1‑149 (MQ=255)
                                    aaCCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCgg                                                                                                              <  7:53697/149‑1 (MQ=255)
                                     aCCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             <  1:200709/149‑1 (MQ=255)
                                     aCCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             >  4:365612/1‑149 (MQ=255)
                                        tCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGa                                                                                                          >  2:69278/1‑149 (MQ=255)
                                        tCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTACTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGa                                                                                                          >  7:355654/1‑149 (MQ=255)
                                                gtcagtcaGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             <  7:87395/138‑1 (MQ=255)
                                                gtcagtcaGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             >  8:87395/1‑138 (MQ=255)
                                                 tcagtcagGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGGGCCAGTGCTCGGGGGACTCGAGGCt                                                                                                 >  6:255046/1‑149 (MQ=255)
                                                     tcagGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGcc                                                                                             >  6:93545/1‑149 (MQ=255)
                                                       aggaaggaGGTTCAATCAGGATAGGCCATCCTGAGCGTGGCAGCTTCACCGCGCGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCAGCCTGTTCGAGCCAGTGCTCGGCGGACTCGAGGCTGGCCTg                                                                                           >  8:291770/1‑149 (MQ=255)
                                                       aggaaggaGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTACTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCt                                                                                            <  8:355654/148‑1 (MQ=255)
                                                            ggaggTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAgg                                                                                      >  8:230303/1‑149 (MQ=255)
                                                             gaggTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             <  2:330765/125‑1 (MQ=255)
                                                              aggTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCTGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGc                                                                                    <  6:210316/149‑1 (MQ=255)
                                                              aggTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGc                                                                                    <  1:314666/149‑1 (MQ=255)
                                                                gTTCAATCAGGAGAGGCCAGCCTGAGCGTGGCAGCGTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGa                                                                                  <  8:148130/149‑1 (MQ=255)
                                                                     aTCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGc                                                                             >  1:315089/1‑149 (MQ=255)
                                                                      tCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGTCCTAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGc                                                                             <  1:127568/149‑1 (MQ=255)
                                                                          gATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTc                                                                                                       >  7:143742/1‑118 (MQ=255)
                                                                          gATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTc                                                                                                       <  8:143742/118‑1 (MQ=255)
                                                                            tAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACt                                                                      <  5:14808/149‑1 (MQ=255)
                                                                               gCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAg                                                                   <  8:254757/149‑1 (MQ=255)
                                                                               gCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAg                                                                   <  3:362995/149‑1 (MQ=255)
                                                                               gCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAg                                                                   >  6:220444/1‑149 (MQ=255)
                                                                                ccAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGGGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAgg                                                                  <  5:255046/149‑1 (MQ=255)
                                                                                ccAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAgg                                                                  <  4:11825/149‑1 (MQ=255)
                                                                                 cAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGc                                                                 >  5:283111/1‑149 (MQ=255)
                                                                                             ggCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGt                                                     >  7:3994/1‑149 (MQ=255)
                                                                                             ggCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGt                                                     <  5:220444/149‑1 (MQ=255)
                                                                                              gCAGCTTCACCGCTGGAAAAATTGACCTGATGGCGGACAGAGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGACAGTGCTCGGGGGACTCGAGGCTTGCCTGGAAGGGAGACAGGCCGATACTCAGGCTGATGCTCAGGTg                                                    <  6:283111/149‑1 (MQ=255)
                                                                                                  cTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCgg                                                 >  2:373400/1‑148 (MQ=255)
                                                                                                   ttCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             <  1:189369/87‑1 (MQ=255)
                                                                                                   ttCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCggg                                                                                                             >  2:189369/1‑87 (MQ=255)
                                                                                                     cACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAAc                                             >  5:112432/1‑149 (MQ=255)
                                                                                                        cgcgGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACtgt                                          <  1:287510/149‑1 (MQ=255)
                                                                                                            ggCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGt                                      >  8:11806/1‑149 (MQ=255)
                                                                                                              gCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGtt                                     >  8:175575/1‑149 (MQ=255)
                                                                                                                    ttGTCCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTTGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTa                               <  1:169431/149‑1 (MQ=255)
                                                                                                                    ttGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCgg                                 >  5:32319/1‑147 (MQ=255)
                                                                                                                       aCCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCt                            >  2:45319/1‑149 (MQ=255)
                                                                                                                         ctcaTCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTACGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTCATGAGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTg                           >  4:74988/4‑148 (MQ=255)
                                                                                                                         cTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTg                           >  5:78733/1‑148 (MQ=255)
                                                                                                                                   ccAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGTGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTc                <  1:373400/149‑1 (MQ=255)
                                                                                                                                      gCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTcgcg             <  6:32319/149‑1 (MQ=255)
                                                                                                                                       cGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTc                <  4:113176/145‑1 (MQ=255)
                                                                                                                                       cGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTc                >  3:113176/1‑145 (MQ=255)
                                                                                                                                        gTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTa           <  6:78733/149‑1 (MQ=255)
                                                                                                                                        gTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTa           <  3:365612/149‑1 (MQ=255)
                                                                                                                                           cTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGa         <  7:291770/148‑1 (MQ=255)
                                                                                                                                           cTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGAc        <  2:103636/149‑1 (MQ=255)
                                                                                                                                                 gATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGt                                                     <  6:332619/98‑1 (MQ=255)
                                                                                                                                                 gATGTACAGCGCCTTGTCAGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATCCGCAGGTGCGGTAACTGTGGATTGAGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCa  >  6:89274/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCA  >  NC_002947/1322716‑1323008

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: