Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,792,063 T→C 6.9% Y159C (TAC→TGC)  PP_2444 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,792,0630TC 6.9% 166.8 / 2.6 58Y159C (TAC→TGC) PP_2444LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

AGGCGGCTGCGTACCACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGTGCCAG  >  NC_002947/2791922‑2792211
                                                                                                                                             |                                                                                                                                                    
agGCGGCTGCGTACAACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCAAGAGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACtt                                                                                                                                               >  2:290066/1‑149 (MQ=255)
         cGTACCACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCa                                                                                                                                       <  1:171716/148‑1 (MQ=255)
          gTACCACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAgg                                                                                                                                     >  4:132660/1‑149 (MQ=255)
          gTACCACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAgg                                                                                                                                     <  1:330960/149‑1 (MQ=255)
             ccACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGa                                                                                                                                             >  5:131661/1‑138 (MQ=255)
             ccACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGa                                                                                                                                             <  6:131661/138‑1 (MQ=255)
             ccACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAggtg                                                                                                                                   <  5:2654/148‑1 (MQ=255)
              cACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAggtggt                                                                                                                                 >  6:141874/1‑149 (MQ=255)
                  gTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGTAGCTCAGCCGGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAgg                                                                                                                             <  6:64892/149‑1 (MQ=255)
                   tctcCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAggt                                                                                                                                    <  6:239173/141‑1 (MQ=255)
                   tctcCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAggt                                                                                                                                    >  5:239173/1‑141 (MQ=255)
                   tctcCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGc                                                                                                                            <  8:118834/149‑1 (MQ=255)
                     tcCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCAt                                                                                                                          >  2:50953/1‑149 (MQ=255)
                          gCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATggg                                                                                                                       <  7:145803/147‑1 (MQ=255)
                          gCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCAAGGTGGTTAGGCATGGGGc                                                                                                                     >  3:278423/1‑149 (MQ=255)
                           ccGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCt                                                                                                                    >  2:182119/1‑149 (MQ=255)
                              ggAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTc                                                                                                                 >  2:257264/1‑149 (MQ=255)
                                aaCCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCg                                                                                                                >  2:104593/1‑148 (MQ=255)
                                 aCCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGa                                                                                                              >  8:190666/1‑149 (MQ=255)
                                  ccAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTa                                                                                                                               <  8:181049/131‑1 (MQ=255)
                                  ccAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTa                                                                                                                               >  7:181049/1‑131 (MQ=255)
                                           gAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAgg                                                                                                    <  1:50953/149‑1 (MQ=255)
                                                  cTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGat                                                                                                             >  1:217679/1‑133 (MQ=255)
                                                       ggCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTg                                                                                         >  4:239710/1‑148 (MQ=255)
                                                           aGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAg                                                                                    <  1:257264/149‑1 (MQ=255)
                                                           aGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAg                                                                                    <  8:131076/149‑1 (MQ=255)
                                                            gCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGc                                                                                   >  2:42882/1‑149 (MQ=255)
                                                                gAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAg                                                                                <  1:42882/148‑1 (MQ=255)
                                                                     ggACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGc                                                                          >  5:106181/1‑149 (MQ=255)
                                                                     ggACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGc                                                                          <  7:295247/149‑1 (MQ=255)
                                                                     ggACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGc                                                                          <  3:132660/149‑1 (MQ=255)
                                                                      gACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGcc                                                                         <  7:190666/149‑1 (MQ=255)
                                                                      gACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGCAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCgg                                                                                                               <  3:69219/111‑1 (MQ=255)
                                                                      gACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGCAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCgg                                                                                                               >  4:69219/1‑111 (MQ=255)
                                                                      gACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGCAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCgg                                                                                                               >  4:73726/1‑111 (MQ=255)
                                                                      gACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGCAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCgg                                                                                                               <  3:73726/111‑1 (MQ=255)
                                                                          ttCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGt                                                                      >  5:304197/1‑148 (MQ=255)
                                                                           tCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTaa                                                                    >  1:339766/1‑149 (MQ=255)
                                                                             aaGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAAcc                                                                  <  3:291726/149‑1 (MQ=255)
                                                                                             gcCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCATGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTc                                                  <  1:290066/149‑1 (MQ=255)
                                                                                             gcCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCa                                                                        >  7:165079/1‑127 (MQ=255)
                                                                                             gcCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCa                                                                        <  8:165079/127‑1 (MQ=255)
                                                                                                ggggTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTcgcg                                               <  3:253058/149‑1 (MQ=255)
                                                                                                   gTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCa                                            <  5:41832/149‑1 (MQ=255)
                                                                                                    tGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAg                                           <  1:104593/149‑1 (MQ=255)
                                                                                                              gCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATgcg                                   >  2:155705/1‑147 (MQ=255)
                                                                                                                   gCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCAcc                            >  1:12797/1‑149 (MQ=255)
                                                                                                                         aCTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCtcctcc                      <  2:339766/149‑1 (MQ=255)
                                                                                                                          cTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCa                     <  3:239710/149‑1 (MQ=255)
                                                                                                                                cAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCa               <  5:141874/149‑1 (MQ=255)
                                                                                                                                   cTTTTGGCTGTAGCACTTTATGTGCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATg             <  4:278423/148‑1 (MQ=255)
                                                                                                                                    ttttGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTc                                                     >  6:234693/1‑107 (MQ=255)
                                                                                                                                    ttttGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGAGATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTc                                                     <  5:234693/107‑1 (MQ=255)
                                                                                                                                        ggCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGt       >  1:239831/1‑149 (MQ=255)
                                                                                                                                         gCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCgcag                                >  4:24311/1‑123 (MQ=255)
                                                                                                                                         gCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCgcag                                <  3:24311/123‑1 (MQ=255)
                                                                                                                                           tGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGTGc     >  2:313224/1‑148 (MQ=255)
                                                                                                                                            gTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGTGc     >  1:323131/1‑147 (MQ=255)
                                                                                                                                             tAGGACTTGTTGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAg                                                                                <  6:184493/71‑1 (MQ=255)
                                                                                                                                             tAGCCCGTGATGTCCAGGTGGGTAGGCCGGGGGCGGGCGGATATCACCAGGTCCAGGCGCGGGATAGCCAGGTCTGCCAGTCCCCGCGCAAGTTGGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGGCCATGGTCGGTGCcag  <  5:347795/149‑1 (MQ=255)
                                                                                                                                             tAGCACTTGATGTCCAGGTGGTTAGGCATTGGTCTGTCGGATATCACAAGGTCCAGGCGCTGGATAGCCAg                                                                                >  5:184493/1‑71 (MQ=255)
                                                                                                                                             tAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGTGCcag  <  6:304197/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGGCGGCTGCGTACCACGGTCTCCTGGCCGGGAACCAGCAGCGGAGCATCCTGCAGGCAAGCTGGAAAAGGACCTTCAAGCCGTTGTGCCAGCGCCGGGGTGGCGAAGAAGCTCAGCCCGCACTCGCCCAGCTTTTGGCTGTAGCACTTGATGTCCAGGTGGTTAGGCATGGGGCTGTCGGATATCACCAGGTCCAGGCGCTGGATAGCCAGGTCTGCCAGTAACCGCTCAAGTTTGTCCTCGCGGCAGTTGATGCGCAGCACCTCCTCCAGGTCCATGGTCGGTGCCAG  >  NC_002947/2791922‑2792211

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