Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 138,952 T→C 6.8% R518G (AGA→GGA) ‡ kinB ← alginate biosynthesis sensor protein KinB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947138,9520TC 6.8% 214.3 / 4.3 73R518G (AGA→GGA) ‡kinBalginate biosynthesis sensor protein KinB
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGCTT  >  NC_002947/138805‑139096
                                                                                                                                                   |                                                                                                                                                 
tttCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTg                                                                                                                                                  >  1:213921/1‑149 (MQ=255)
ctaccTTCCCCAGCGCTAGCCCCATCCCACCCCCGCCTTTCATGCTACATTCCTGAACTAAAGGCGCTAAGATGCTCGCTTGATGCCCATAGGAAATGCCCTCTCCTTTTTCCTCGACGCTGATGATAACCCGTTCTGCATGCCGTCTg                                                                                                                                                  <  4:251432/146‑1 (MQ=255)
     ttGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTgggagg                                                                                                                                             >  7:10426/1‑145 (MQ=255)
             cgcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGa                                                                                                                                     <  2:213921/149‑1 (MQ=255)
              gcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGAt                                                                                                                                    <  1:79180/149‑1 (MQ=255)
              gcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCAg                                                                                                                                          >  8:86081/1‑143 (MQ=255)
              gcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCAg                                                                                                                                          <  7:86081/143‑1 (MQ=255)
                      cagcccagcGCCGCCTTTCATGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGcgcccgc                                                                                                                            >  5:228240/1‑149 (MQ=255)
                         cccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTAGTTGCCCAAAGGAAATGCCCTCTCCGTTGTCCTCGACGCAGATCATCACCCGTTCTGCATGCCGTCTCGCATGCAGCCGGATGCGCCCGCCGt                                                                                                                         >  5:286014/1‑149 (MQ=255)
                          ccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTTCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGtt                                                                                                                        <  2:245563/149‑1 (MQ=255)
                            agcGCCGCCGTTCTTGCGCCCTACCTGCACGACAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGGCCTCGACGCTGATGATCACCCTTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTgg                                                                                                                      <  8:275812/149‑1 (MQ=255)
                              cgccgcCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCg                                                                                                                    <  3:248921/149‑1 (MQ=255)
                               gccgccTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCgg                                                                                                                   >  8:267123/1‑149 (MQ=255)
                                         ttGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCACATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTATGCATGCCGTCTGTCATGAAGGCGGATGCGCCCGCCGTTGGCGTTATGGCGGAt                                                                                                        >  6:261160/1‑149 (MQ=255)
                                           gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATTGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                         <  6:286014/146‑1 (MQ=255)
                                           gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                         >  6:74772/1‑146 (MQ=255)
                                           gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                         <  7:267123/146‑1 (MQ=255)
                                           gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                         <  5:74772/146‑1 (MQ=255)
                                                    ctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGACCGCCGATGGCGGTATGGAGGATGGCGTtatg                                                                                               >  8:221213/1‑147 (MQ=255)
                                                       cACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATTGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                          <  2:339785/149‑1 (MQ=255)
                                                       cACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                          <  5:261160/149‑1 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGtt                                                                                                  <  1:282366/138‑1 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGtt                                                                                                  >  2:282366/1‑138 (MQ=255)
                                                               gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTc                                                                                   >  2:72782/1‑148 (MQ=255)
                                                                cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGCTATGGCGGATGGCGTTGTGCAGCAGGTTGTcc                                                                                  <  1:216524/149‑1 (MQ=255)
                                                                cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGATATGGCGGATGGCGTTGTGCAGCAGGTTGTcc                                                                                  <  1:211576/149‑1 (MQ=255)
                                                                          cgggcTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCTGGTTGTCCAGCACTCGGTc                                                                       <  6:228240/146‑1 (MQ=255)
                                                                          cgcgcTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGt                                                                        <  6:122460/148‑1 (MQ=255)
                                                                           gcgcTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTc                                                                       >  1:158711/1‑148 (MQ=255)
                                                                                 gttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGcagcag                                                                                         >  6:141647/1‑124 (MQ=255)
                                                                                 gttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGcagcag                                                                                         <  5:141647/124‑1 (MQ=255)
                                                                                  ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgt                                                                                                <  2:148981/116‑1 (MQ=255)
                                                                                  ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgt                                                                                                >  1:148981/1‑116 (MQ=255)
                                                                                  ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtg                                                                                                 >  1:124566/1‑115 (MQ=255)
                                                                                  ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtg                                                                                                 <  2:124566/115‑1 (MQ=255)
                                                                                    gCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGc                                                              >  5:336257/1‑148 (MQ=255)
                                                                                      ccATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCaa                                                                    <  2:70240/140‑1 (MQ=255)
                                                                                      ccATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCaa                                                                    >  1:70240/1‑140 (MQ=255)
                                                                                        aTAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGt                                                         >  2:235320/1‑149 (MQ=255)
                                                                                              aaTGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCt                                                   >  7:273415/1‑149 (MQ=255)
                                                                                                  cccTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGAt                                               <  2:158711/149‑1 (MQ=255)
                                                                                                   ccTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATa                                              >  3:15336/1‑149 (MQ=255)
                                                                                                   ccTCTCCGTTGTCCTCGACGCCGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATa                                              >  4:209135/1‑149 (MQ=255)
                                                                                                    ctctCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACCCGGTCCAACTGGGCGACGTCGGCCTGGATAc                                             >  2:34175/1‑149 (MQ=255)
                                                                                                       tcCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                          >  1:162423/1‑149 (MQ=255)
                                                                                                        ccGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                          <  4:15336/148‑1 (MQ=255)
                                                                                                        ccGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCgg                                         >  1:322150/1‑149 (MQ=255)
                                                                                                         cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                          <  6:336257/147‑1 (MQ=255)
                                                                                                              tCCTCGACGCTGATGATCACCCGTGCTGCATGCCGTCTGGCATGCAGGCGTATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcgg                                   <  8:197477/149‑1 (MQ=255)
                                                                                                               ccTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGATGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGACGGCCTGCATACCCGGTAgcggc                                  >  1:337085/1‑149 (MQ=255)
                                                                                                                 tCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGcggcgg                                >  2:311339/1‑149 (MQ=255)
                                                                                                                   gACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATTGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGcggcgg                                <  7:221213/147‑1 (MQ=255)
                                                                                                                   gACGCCGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGcggcgg                                <  3:209135/147‑1 (MQ=255)
                                                                                                                             atCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTTTGGCGGATGGCGTGGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGGCTGGATACGCTGTAGCGGCGGTTCCAGTTCCtt                    <  6:248416/149‑1 (MQ=255)
                                                                                                                              tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGTAGGATCGCGTTCTGCAGCAGGTTGTCCAGCACTCGGTCGAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTGCCGTTACTTg                   >  5:219475/1‑149 (MQ=255)
                                                                                                                              tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                   <  2:162423/149‑1 (MQ=255)
                                                                                                                              tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                   <  7:274060/149‑1 (MQ=255)
                                                                                                                               cACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTTGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGa                  >  3:182785/1‑149 (MQ=255)
                                                                                                                               cACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                   <  8:273415/148‑1 (MQ=255)
                                                                                                                                aCCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGa                  <  1:311339/148‑1 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCg                                                                          >  3:141769/1‑91 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCg                                                                          <  4:141769/91‑1 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTCGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATc                <  4:182785/149‑1 (MQ=255)
                                                                                                                                      ccTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTc           >  2:226775/2‑149 (MQ=255)
                                                                                                                                       cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTggg                                                               <  4:332025/96‑1 (MQ=255)
                                                                                                                                       cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTggg                                                               >  3:332025/1‑96 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCGTGGATACGCGGTAGCGGCGGGGCCAGTTCCTTGCTCAGCTCGAt        >  7:294552/1‑149 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGt                              >  2:101156/1‑127 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGt                              <  1:101156/127‑1 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGGGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAAc             >  5:235461/1‑144 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGCGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAAc             <  6:235461/144‑1 (MQ=255)
                                                                                                                                             gCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGc    >  1:249048/1‑149 (MQ=255)
                                                                                                                                               cGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGCtt  >  2:251125/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGCTT  >  NC_002947/138805‑139096

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: