Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,155,470 A→T 8.0% I210F (ATC→TTC)  PP_1013 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,155,4700AT 8.0% 219.2 / 7.3 75I210F (ATC→TTC) PP_1013sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (3/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CCATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAGGCGG  >  NC_002947/1155325‑1155618
                                                                                                                                                 |                                                                                                                                                    
ccATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCg                                                                                                                                                   <  3:361390/149‑1 (MQ=255)
   ttattCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCt                                                                                                                                                >  7:297711/1‑149 (MQ=255)
     attCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCtt                                                                                                                                               <  8:297711/148‑1 (MQ=38)
      ttCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTc                                                                                                                                              <  1:36521/148‑1 (MQ=38)
       tCGCTGACGATGGATCCGCTGATACCCCCGGAACTGGTCACTCTGATCCGCTTTGGCGAGGTCGAATGGTTGTACCTCGCCTCTCTGTTGCCCGAACCTTTCACGGGGCTGGTAGCCGAACGCCTGCCGCGCCAGCTGATCGGCttga                                                                                                                                             <  6:373737/148‑3 (MQ=255)
         gCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCg                                                                                                                                          >  6:40697/1‑149 (MQ=255)
         gCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCg                                                                                                                                          >  4:372858/1‑149 (MQ=255)
             aCCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCt                                                                                                                                      >  6:314708/1‑149 (MQ=255)
                    aGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTGGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACg                                                                                                                               <  1:261361/149‑1 (MQ=255)
                      ccGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGc                                                                                                                             >  3:45979/1‑149 (MQ=255)
                        gCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCt                                                                                                                           >  8:56662/1‑149 (MQ=255)
                        gCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCt                                                                                                                           >  8:230817/1‑149 (MQ=255)
                          tGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGAt                                                                                                                                                     >  1:355335/1‑121 (MQ=255)
                          tGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGAt                                                                                                                                                     <  2:355335/121‑1 (MQ=255)
                               cccccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgtt                                                                                                                    >  3:224234/1‑149 (MQ=255)
                                ccccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                   <  8:83752/149‑1 (MQ=255)
                                ccccGGTACTGGGCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCGGTTGCCCGAACCGTACACGGGGCGGGAAGCCGAACGCCTGCCGCGCCAGCAGAGCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                   <  3:42831/149‑1 (MQ=255)
                                 cccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgc                                                                                                                  >  8:231800/1‑149 (MQ=255)
                                  ccGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgct                                                                                                                 <  5:303397/149‑1 (MQ=255)
                                   cGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttgctg                                                                                                                <  7:230817/149‑1 (MQ=255)
                                         tGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATc                                                                                                          >  8:267671/1‑149 (MQ=255)
                                           gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                        >  7:36529/1‑149 (MQ=255)
                                           gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                        >  3:170024/1‑149 (MQ=255)
                                           gTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCg                                                                                                         <  5:40697/148‑1 (MQ=255)
                                            tCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGc                                                                                                       <  3:216009/149‑1 (MQ=255)
                                                gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGcc                                                                                                                            <  8:321360/124‑1 (MQ=255)
                                                gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGcc                                                                                                                            >  7:321360/1‑124 (MQ=255)
                                                gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                   >  5:12965/1‑133 (MQ=255)
                                                gCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttgctgttg                                                                                                                   <  6:12965/133‑1 (MQ=255)
                                                        gCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTTCAGTTCATCGGCTTGCTGGTGCAc                                                                                           <  7:56662/149‑1 (MQ=255)
                                                        gCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTg                                                                                              <  7:191915/146‑1 (MQ=255)
                                                        ccTTGGGCGAGGGCGAATCGATGTACATCGACTCGCTGTTGCCAGAACCGTACACGCGGCTCGAAGCCTAAAGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGATGACAACGGCCTTGATGTTGATGTTCATCGCCTTGCTGGTg                                                                                              >  8:191915/2‑146 (MQ=255)
                                                             ggcgagggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCa                                                                                      <  4:45979/149‑1 (MQ=255)
                                                             ggcgagggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCa                                                                                      <  7:267671/149‑1 (MQ=255)
                                                               cgagggcgaATGCTTGTACATCGCCTCCCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCATGCCGCGCCAGCAGATCGGCTTCATCGAGCTGACCACGGCCTTGCTGTTGCTGTTCAACGGCTTGCTGGTGCACTGGCaaa                                                                                    >  3:108030/1‑149 (MQ=255)
                                                                   ggcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTg                                                                                <  4:170024/149‑1 (MQ=255)
                                                                    gcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttg                                                                                                                         >  7:343320/1‑107 (MQ=255)
                                                                    gcgaATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCttg                                                                                                                         <  8:343320/107‑1 (MQ=255)
                                                                               tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTg                                                                                                    <  8:245149/117‑1 (MQ=255)
                                                                               tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTg                                                                                                    >  7:245149/1‑117 (MQ=255)
                                                                               tACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAgc                                                                     <  7:231800/148‑1 (MQ=255)
                                                                                 cATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATc                                                                                                                                           <  3:21601/76‑1 (MQ=255)
                                                                                 cATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATc                                                                                                                                           >  4:21601/1‑76 (MQ=255)
                                                                                  aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCATGGCCTTGCTGTTGCTGTTCATCGGCTTGCTg                                                                                                 <  5:94387/117‑1 (MQ=255)
                                                                                  aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCATGGCCTTGCTGTTGCTGTTCATCGGCTTGCTg                                                                                                 >  6:94387/1‑117 (MQ=255)
                                                                                  aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCt                                                                 >  4:317159/1‑149 (MQ=255)
                                                                                  aTCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCt                                                                 >  2:256132/1‑149 (MQ=255)
                                                                                    cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGGTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTgg                                                               >  6:342223/1‑149 (MQ=255)
                                                                                    cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                              >  7:12712/1‑134 (MQ=255)
                                                                                    cGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                              <  8:12712/134‑1 (MQ=255)
                                                                                                gCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTg                                                                <  4:217976/136‑1 (MQ=255)
                                                                                                gCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTg                                                                >  3:217976/1‑136 (MQ=255)
                                                                                                 cccGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                        <  5:373540/95‑1 (MQ=255)
                                                                                                 cccGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCgg                                                                                                        >  6:373540/1‑95 (MQ=255)
                                                                                                   cGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCgag                                                <  3:317159/149‑1 (MQ=255)
                                                                                                    gAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCgaga                                               <  5:314708/149‑1 (MQ=255)
                                                                                                     aaCCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGc                                                                              >  4:166295/1‑117 (MQ=255)
                                                                                                             acGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                      >  6:329156/1‑149 (MQ=255)
                                                                                                             acGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                      >  4:352904/1‑149 (MQ=255)
                                                                                                               ggggCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTgg                                    >  4:300806/1‑149 (MQ=255)
                                                                                                                  cgatgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTg                                 >  3:367442/5‑149 (MQ=38)
                                                                                                                   cTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTAGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGc                                <  8:36529/149‑1 (MQ=255)
                                                                                                                    tGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGcc                               >  5:345677/1‑149 (MQ=255)
                                                                                                                    atgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                      <  6:44469/140‑1 (MQ=38)
                                                                                                                    atgAAGCCGATCTGCTGGCGCGGCAGGCGTTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCt                                      >  5:44469/3‑142 (MQ=38)
                                                                                                                      gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                             >  4:301402/1‑149 (MQ=255)
                                                                                                                      gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                             <  4:224234/149‑1 (MQ=255)
                                                                                                                      gAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGa                             >  1:193765/1‑149 (MQ=255)
                                                                                                                            gAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGgc                       >  4:363892/1‑149 (MQ=255)
                                                                                                                                  cTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGt                 <  1:256132/149‑1 (MQ=255)
                                                                                                                                         gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTg            >  6:68114/1‑147 (MQ=255)
                                                                                                                                         gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa          <  7:46702/149‑1 (MQ=255)
                                                                                                                                         gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa          <  5:329156/149‑1 (MQ=255)
                                                                                                                                         gcCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGaa          <  2:193765/149‑1 (MQ=255)
                                                                                                                                               agATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAggc    >  3:153971/1‑149 (MQ=255)
                                                                                                                                                 aTCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAggcgg  <  3:363892/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCATTATTCGCTGACCTTGGAGCCGCTGAAACCCCCGGTACTGGTCACGCAGATCCGCTTGGGCGAGGGCGAATGGTTGTACATCGCCTCGCTGTTGCCCGAACCGTACACGGGGCTGGAAGCCGAACGCCTGCCGCGCCAGCAGATCGGCTTCATCGTGCTGACCACGGCCTTGCTGTTGCTGTTCATCGGCTTGCTGGTGCACTGGCAAAGCTGGCCACTGAAGCGCCTGGCGCGGGCCGCCCGAGAGATGTCGCTGGGTGCCGATGTGGCGCCAGTGGCTGAAGCAGGCGG  >  NC_002947/1155325‑1155618

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