Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,405 G→C 10.2% G94A (GGC→GCC)  PP_1181 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,4050GC10.2% 199.5 / 11.5 69G94A (GGC→GCC) PP_1181two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (3/4);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAC  >  NC_002947/1356257‑1356550
                                                                                                                                                     |                                                                                                                                                 
atcatcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAggg                                                                                                                                                   >  5:266514/1‑149 (MQ=255)
    tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTg                                                                                                                                               >  7:96907/1‑149 (MQ=255)
    tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTg                                                                                                                                               <  6:325210/149‑1 (MQ=255)
       aCAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGCCCTGGACCTggtc                                                                                                                                            <  3:214967/149‑3 (MQ=255)
        cagccagcCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCGTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAcc                                                                                                                                           <  2:223463/149‑1 (MQ=255)
           ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTgg                                                                                                                                        <  3:58455/149‑1 (MQ=255)
            cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTggt                                                                                                                                       <  6:329062/149‑2 (MQ=255)
              gccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGTccagttca                                                                                                                                     <  5:328430/149‑7 (MQ=255)
                 cGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGcc                                                                                                                                   >  4:38354/1‑148 (MQ=255)
                  gACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcga                                                                                                                                 >  5:143172/1‑149 (MQ=255)
                              cTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCt                                                                                                                      <  3:306621/148‑1 (MQ=255)
                               tCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGCAGACTACGTCTg                                                                                                                     >  8:108407/1‑148 (MQ=255)
                               tCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTg                                                                                                                     >  5:298431/1‑148 (MQ=255)
                               tCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGAccg                                                                                                                                          >  5:317014/1‑126 (MQ=255)
                               tCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGAccg                                                                                                                                          <  6:317014/127‑2 (MQ=255)
                                 gACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGc                                                                                                                    <  2:115797/147‑1 (MQ=255)
                                  aCCTGATGCTCCCCGGCGAGGACGGCCTGACCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCaaa                                                                                                                 >  1:26880/1‑149 (MQ=255)
                                       aTGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGt                                                                                                            <  8:96907/149‑1 (MQ=255)
                                          cTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgc                                                                                                          <  2:26880/148‑1 (MQ=255)
                                           tCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgcgc                                                                                                        >  4:287268/1‑149 (MQ=255)
                                                         ggCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTgctgct                                                                                          >  5:262013/1‑149 (MQ=255)
                                                           ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTgg                                                                                        <  6:287022/149‑1 (MQ=255)
                                                           ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTgg                                                                                        >  6:80261/1‑149 (MQ=255)
                                                           ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTgg                                                                                        <  1:43544/149‑1 (MQ=255)
                                                            cTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGc                                                                                       >  3:167311/1‑149 (MQ=255)
                                                            cTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGc                                                                                       <  8:45937/149‑1 (MQ=255)
                                                             tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCa                                                                                      >  3:66641/1‑149 (MQ=255)
                                                             tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTgg                                                                                                                                        >  2:164893/1‑99 (MQ=255)
                                                             tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTgg                                                                                                                                        <  1:164893/99‑1 (MQ=255)
                                                              tcgcttCTGCCGCCGGGGGCGCAGCAAGTACCAGGGGCGGATCCTCATGCTCACCGCGCGCAGTTATGAACTGCACCAGGCCCCGGGCCTGCACGTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCAc                                                                                     <  2:148043/144‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCAc                                                                                                    <  5:231639/129‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCAc                                                                                                    >  6:231639/1‑129 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATc                                                                                >  3:272200/1‑149 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATc                                                                                >  2:4179/1‑149 (MQ=255)
                                                                       ccgccgGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATcc                                                                               <  6:143172/146‑1 (MQ=255)
                                                                       ccgccgGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATc                                                                                <  3:38354/145‑1 (MQ=255)
                                                                                       ccgaaCCAGGGGCCGATCCTCATGCACACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                            <  4:191971/145‑1 (MQ=255)
                                                                                         gTACCAGGGGCCGATCCTCATGCTCACCGCGCGCCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                            >  1:204437/1‑148 (MQ=255)
                                                                                         gTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                            <  2:204437/147‑1 (MQ=255)
                                                                                              aGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAcc                                                     >  2:228985/1‑149 (MQ=255)
                                                                                                gggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                                   <  6:298431/149‑1 (MQ=255)
                                                                                                   ccGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGa                                                >  4:72807/1‑149 (MQ=255)
                                                                                                    cGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGa                                                >  2:92326/1‑148 (MQ=255)
                                                                                                         cTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCaa                                          <  4:66641/149‑1 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                            >  7:82780/1‑130 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                            <  8:82780/130‑1 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAg                                         >  8:170021/1‑149 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAg                                         >  7:203971/1‑149 (MQ=255)
                                                                                                            aTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgc                                        <  6:262013/148‑1 (MQ=255)
                                                                                                            aTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAg                                         <  3:287268/147‑1 (MQ=255)
                                                                                                             tGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGTCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcgc                                      >  5:229005/1‑149 (MQ=255)
                                                                                                                        cgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGc                           >  5:31996/1‑149 (MQ=255)
                                                                                                                              tgatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTTCGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                      >  6:323743/1‑149 (MQ=255)
                                                                                                                              tgatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                      <  6:123863/149‑1 (MQ=255)
                                                                                                                               gatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCAAGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCgg                     >  3:112753/1‑149 (MQ=255)
                                                                                                                                  gaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                            <  3:199283/139‑1 (MQ=255)
                                                                                                                                  gaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                            >  4:199283/1‑139 (MQ=255)
                                                                                                                                   aaCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgcg                 >  8:82077/1‑149 (MQ=255)
                                                                                                                                    aCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgcgc                >  2:234707/1‑149 (MQ=255)
                                                                                                                                      tGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTg              >  3:217148/1‑149 (MQ=255)
                                                                                                                                        gcccAGGTCCAGGCCCTGGACCTGGTCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCg            >  8:57785/3‑149 (MQ=255)
                                                                                                                                        gACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc             <  1:92326/148‑1 (MQ=255)
                                                                                                                                        gACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc             <  1:28326/148‑1 (MQ=255)
                                                                                                                                        gACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCg            >  2:164894/1‑149 (MQ=255)
                                                                                                                                        gACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCCAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc             <  1:228985/148‑1 (MQ=255)
                                                                                                                                             ggTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc        >  3:161483/1‑148 (MQ=255)
                                                                                                                                              gTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGa      >  7:126590/1‑149 (MQ=255)
                                                                                                                                                ccAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                            <  8:122616/125‑1 (MQ=255)
                                                                                                                                                ccAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                            >  7:122616/1‑125 (MQ=255)
                                                                                                                                                  aGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAACCCCGTGCGCCCCCGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAc  >  5:182033/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAC  >  NC_002947/1356257‑1356550

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: