Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,532,326 T→G 6.8% E286A (GAA→GCA)  PP_3121 ← transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,532,3260TG 6.8% 184.0 / 4.1 59E286A (GAA→GCA) PP_3121transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base G (2/2);  total (28/31)
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.40e-01

CCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGG  >  NC_002947/3532178‑3532467
                                                                                                                                                    |                                                                                                                                             
ccTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACtt                                                                                                                                               >  5:283041/1‑149 (MQ=255)
   gTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTc                                                                                                                                              >  7:231111/1‑147 (MQ=255)
      gCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAg                                                                                                                                           <  3:369188/147‑1 (MQ=255)
      gCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAg                                                                                                                                           >  4:369188/1‑147 (MQ=255)
       cAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAg                                                                                                                                        >  7:358201/1‑149 (MQ=255)
              gcgACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTGCTTGAAGTGCttct                                                                                                                                 <  8:323053/149‑3 (MQ=255)
               cgACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTtgtg                                                                                                                                >  4:142471/1‑149 (MQ=255)
                 aCTTAGGCAGCCTTGCGTCCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGtgct                                                                                                                              <  4:253694/149‑1 (MQ=255)
                      ggCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCgg                                                                                                                          <  2:376503/148‑1 (MQ=255)
                       gCAGCCTTGCGACCTGTGTCTTCGCACATGTACCAGGTGCCTTGGCCGATGGCTGCAACCCTCGTGCATGCCAGATTTACGTAAAGCATGTCACCTCCGGGATACGGGTATTGTAAAGAAGCACTTCAAGCAGTCCTTGTGCTGCGGTg                                                                                                                        >  6:100369/1‑149 (MQ=255)
                        cAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTgg                                                                                                                       <  7:112295/149‑1 (MQ=255)
                           ccTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGaaa                                                                                                                    <  1:198515/149‑1 (MQ=255)
                               gCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                >  3:163940/1‑149 (MQ=255)
                               gCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCt                                                                                                                 <  4:202382/148‑1 (MQ=255)
                                cGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                <  6:283041/148‑1 (MQ=255)
                                cGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                <  2:162535/148‑1 (MQ=255)
                                 gACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtgt                                                                                                              <  8:358201/149‑1 (MQ=255)
                                   ccTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGc                                                                                                            <  3:158493/149‑1 (MQ=255)
                                      gggTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCt                                                                                                           <  2:193484/147‑1 (MQ=255)
                                      gggTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCt                                                                                                           >  1:193484/1‑147 (MQ=255)
                                         tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCt                                                                                                       >  5:76973/1‑148 (MQ=255)
                                         tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTg                                                                                                      <  3:142471/149‑1 (MQ=255)
                                         tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCGTGTGCTGCGGTGGAAAGCTTGTGCTCGCTg                                                                                                      >  4:71728/1‑149 (MQ=255)
                                              gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtg                                                                                                               >  7:47787/1‑135 (MQ=255)
                                              gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtg                                                                                                               <  8:47787/135‑1 (MQ=255)
                                              gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAggg                                                                                                  >  5:322987/1‑148 (MQ=255)
                                              gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGACTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAggg                                                                                                  >  2:17945/1‑148 (MQ=255)
                                               cccATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGc                                                                                                 >  5:7695/1‑148 (MQ=255)
                                                  aTGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATc                                                                                             <  4:132033/149‑1 (MQ=255)
                                                       ccAGCTGCCTTGGCCGATGGCTGCAACGCTCCTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                         <  4:27061/148‑1 (MQ=255)
                                                           gTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                         <  1:131480/144‑1 (MQ=255)
                                                           gTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                         >  2:131480/1‑144 (MQ=255)
                                                              ccTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATgcgc                                                                                 <  6:7695/149‑1 (MQ=255)
                                                                       aTGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACt                                                                        <  8:231111/149‑1 (MQ=255)
                                                                           cTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTgg                                                                    >  2:87050/1‑149 (MQ=255)
                                                                             gCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTggg                                                                   >  5:113831/1‑148 (MQ=255)
                                                                             gCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTggg                                                                   <  6:113831/148‑1 (MQ=255)
                                                                               aaCGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTc                                                                <  6:76973/149‑1 (MQ=255)
                                                                                  gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTc                                                                            >  5:15808/1‑134 (MQ=255)
                                                                                  gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTc                                                                            <  6:15808/134‑1 (MQ=255)
                                                                                  gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACtgt                                                                      <  4:198344/140‑1 (MQ=255)
                                                                                  gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACtgt                                                                      >  3:198344/1‑140 (MQ=255)
                                                                                       tgcctgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAgg                                                        <  1:87050/149‑1 (MQ=255)
                                                                                         cctgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGTAAAGCTTGTGCTCGCTGAGGGCAATCACGCCTATGCGCCGCTCGACTGTGGGCTCGACCAAGGTg                                                      <  5:100369/149‑1 (MQ=255)
                                                                                           tgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGAc                                                    >  7:24773/1‑149 (MQ=255)
                                                                                                  tttACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCgg                                             <  4:163940/149‑1 (MQ=255)
                                                                                                          aCGCATGTCACCTCCCGGATCCGGGTATTGGAAAGAACCACTTCAACCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAg                                     <  6:322987/149‑1 (MQ=255)
                                                                                                                     cTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCagat                          >  1:100604/1‑146 (MQ=255)
                                                                                                                     cTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCt                          >  3:158442/1‑149 (MQ=255)
                                                                                                                           gATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGa                    <  1:145634/149‑1 (MQ=255)
                                                                                                                                       ggAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcg        >  1:69853/1‑149 (MQ=255)
                                                                                                                                        gAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcgg       >  2:36868/1‑149 (MQ=255)
                                                                                                                                        gAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcgg       >  4:169863/1‑149 (MQ=255)
                                                                                                                                        gAAAGAAGCACTTCAAGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCAGATGGTTGATGCACAGGGCcgg       >  2:160324/1‑149 (MQ=255)
                                                                                                                                          aagaagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCCTCTGGTTGAAGCACAGGGCcggca     >  6:264528/1‑149 (MQ=255)
                                                                                                                                           agaagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTTCATCTGGTTGATGCACAGGGCcggcac    <  1:17945/149‑1 (MQ=255)
                                                                                                                                            caagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGAt                   >  8:90775/2‑133 (MQ=255)
                                                                                                                                            caagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGAt                   <  7:90775/132‑1 (MQ=255)
                                                                                                                                             aagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGACGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCTGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcggcacgg  <  1:160324/149‑1 (MQ=255)
                                                                                                                                             aagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcggcacgg  >  1:233043/1‑149 (MQ=255)
                                                                                                                                                   ttCAAGCAGTGCTTGTGCTGCGGTGTAAAGCTTGTGCTCTCTGATGGCAATCACGCCGATTCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGtt                      <  5:230934/123‑1 (MQ=255)
                                                                                                                                                   ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCTGGCCCCCAGCTCCTGCATCttgat                      >  6:230934/1‑119 (MQ=255)
                                                                                                                                                   ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCt                                   <  4:87852/110‑1 (MQ=255)
                                                                                                                                                   ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCt                                   >  3:87852/1‑110 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGG  >  NC_002947/3532178‑3532467

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