Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,965,662 T→A 7.9% intergenic (+68/‑146) PP_3496 → / → ydcP endoribonuclease L‑PSP/peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,965,6620TA 7.9% 177.7 / 4.6 63intergenic (+68/‑146)PP_3496/ydcPendoribonuclease L‑PSP/peptidase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base A (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATATAACGAC  >  NC_002947/3965518‑3965807
                                                                                                                                                |                                                                                                                                                 
gtgcgGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCgg                                                                                                                                               >  6:374342/1‑149 (MQ=255)
   cggcggTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGccc                                                                                                                                            <  4:135985/149‑1 (MQ=255)
        gtgggtgTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTcg                                                                                                                                       <  1:7466/149‑1 (MQ=255)
                 gCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGa                                                                                                                                    >  4:74463/1‑143 (MQ=255)
                 gCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGa                                                                                                                                    <  3:74463/143‑1 (MQ=255)
                 gCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAg                                                                                                                              >  5:186300/1‑149 (MQ=255)
                 gCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAg                                                                                                                              <  3:285686/149‑1 (MQ=255)
                   ttCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGc                                                                                                                            >  3:71747/1‑149 (MQ=255)
                        tGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTc                                                                                                                       >  7:119605/1‑149 (MQ=255)
                          gCCTATGAAGGGATCCTTGTGGATGGGACGGCTTTGGCTGCCTTGCGCTTAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCTACATGGGTGTTTGGTGGGAAAATCGACGGGGCCGCTCTGCGGCCCTTTCGCGACGCTTGGCCGCTc                                                                                                                       <  4:9714/147‑1 (MQ=255)
                            cTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCca                                                                                                                    >  8:50690/1‑148 (MQ=255)
                            cTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCAGAGCGGccccg                                                                                                                                          <  1:245331/126‑3 (MQ=255)
                            cTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCAGAGCGGccccg                                                                                                                                          >  2:245331/1‑124 (MQ=255)
                             tttGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCAGAGCGGccccgg                                                                                                                                         <  2:111378/126‑4 (MQ=255)
                             tttGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGATGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCAGAGCGGccccgg                                                                                                                                         >  1:111378/1‑123 (MQ=255)
                               tGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCcaca                                                                                                                  <  3:344374/147‑1 (MQ=255)
                                  aGGCATGGGTGTGGAAGTGACGTGGATGGCTTCCAAGCGCTTAGCTTCGCTTGTGGGGGTGCGCAGCGCCCCGACTTCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTg                                                                                                              <  4:71747/148‑1 (MQ=38)
                                    gCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTgcgc                                                                                                           <  7:50690/149‑1 (MQ=255)
                                    gCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTgcgc                                                                                                           <  7:496/149‑1 (MQ=255)
                                    gCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTgcgc                                                                                                           <  5:374342/149‑1 (MQ=255)
                                        gCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTAc                                                                                                       >  4:309600/1‑149 (MQ=255)
                                           tgtgGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCg                                                                                                    >  2:72773/1‑149 (MQ=255)
                                             tgGAAGTGCCGGCTTTGGCTGCCAAGCGCTGAGCGGCGCTTGTGGGGGTGCTCAGGGCCCCGACATCTCTTATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCa                                                                                                  <  8:119605/149‑1 (MQ=38)
                                              ggAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAg                                                                                                 <  2:240301/149‑1 (MQ=255)
                                               gAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAg                                                                                                 <  8:350692/148‑1 (MQ=255)
                                                      cGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATccc                                                                                         >  3:237579/1‑149 (MQ=255)
                                                       ggCGAGGTCTGCCAAGCGCGGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGAGTGGTCGCAAAATCGCCGGGGCCGCTATGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACATTTGCGCGTACCCGCAGCCAATCCCg                                                                                        <  8:4660/149‑1 (MQ=255)
                                                          gATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCAc                                                                                     >  1:274094/1‑149 (MQ=255)
                                                               cTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCtctc                                                                                <  6:186300/149‑1 (MQ=255)
                                                               cTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCtctc                                                                                >  6:325043/1‑149 (MQ=255)
                                                                tGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGtt                                                                                                               >  4:116007/1‑117 (MQ=255)
                                                                tGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGtt                                                                                                               <  3:116007/117‑1 (MQ=255)
                                                                  ccAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTAccc                                                                                                     >  5:260825/1‑125 (MQ=255)
                                                                  ccAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTAccc                                                                                                     <  6:260825/125‑1 (MQ=255)
                                                                  ccAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCgg                                                                              <  2:17997/148‑1 (MQ=255)
                                                                  ccAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCg                                                                               <  2:274094/147‑1 (MQ=255)
                                                                    aaGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCAc                                                                           >  6:89198/1‑149 (MQ=255)
                                                                        gcTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCAc                                                                                     >  4:90065/1‑135 (MQ=255)
                                                                        gcTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCAc                                                                                     <  3:90065/135‑1 (MQ=255)
                                                                         cTGAGCGGCGCTGGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCaa                                                                      <  6:280613/149‑1 (MQ=255)
                                                                            aGCTGCGCTTGTGGGGTTGCTCAGCGCCGCGAGATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGAAGCAAGGCCGCTCGCTCATTTTCGCGTACCCTCAGCCAATCCCGGACGTCTCGGGAGCCCTat                                                                     <  2:106853/147‑1 (MQ=25)
                                                                                    ttGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCa                                                           >  3:106259/1‑149 (MQ=255)
                                                                                            gTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCg                                                                               >  7:228608/1‑121 (MQ=255)
                                                                                            gTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCg                                                                               <  8:228608/121‑1 (MQ=255)
                                                                                                         cGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATg                                      >  4:28986/1‑149 (MQ=255)
                                                                                                          gACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGc                                     >  5:350344/1‑149 (MQ=255)
                                                                                                           aCATCGGTGTTTGGTCGGAAATTCGCCGGGTCGGCTCTGCGGCACTTTCGCGACGCTATGCCGCTCCCACATTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGGCCGTTTCACCGCTACAATGcc                                    <  1:45627/149‑1 (MQ=255)
                                                                                                            cATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGccc                                   <  4:106259/149‑1 (MQ=255)
                                                                                                              tCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGc                                 <  3:333288/149‑1 (MQ=255)
                                                                                                               cGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCAcccc                                                                        >  5:343837/1‑109 (MQ=255)
                                                                                                               cGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCAcccc                                                                        <  6:343837/109‑1 (MQ=255)
                                                                                                                  tGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTa                              >  2:119166/1‑148 (MQ=255)
                                                                                                                        gTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGt                        >  3:242058/1‑148 (MQ=255)
                                                                                                                           gCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGAc                     <  3:309600/148‑1 (MQ=255)
                                                                                                                           gCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACa                    >  8:171558/1‑149 (MQ=255)
                                                                                                                               aaTCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAACGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGccc                <  4:237579/149‑1 (MQ=255)
                                                                                                                                   gaaagggCCGCAGAGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGaa             >  3:31280/5‑148 (MQ=255)
                                                                                                                                       gggCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGGTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAata         >  5:203218/1‑148 (MQ=255)
                                                                                                                                        ggCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGaaa            >  1:343674/1‑144 (MQ=255)
                                                                                                                                        ggCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGaaa            <  2:343674/144‑1 (MQ=255)
                                                                                                                                          ccGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATATaa      <  7:171558/148‑1 (MQ=255)
                                                                                                                                          ccGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATATAAc     <  5:325043/149‑1 (MQ=255)
                                                                                                                                           cGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGCCCCCCCAATGATGCCGCCAGCCCGTTTCACAGCTACAATGCCCGCCTAAAACGTGACAGCCCGAAAATATAACg    >  8:279078/1‑149 (MQ=255)
                                                                                                                                             ctctGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATATAACGAc  <  5:89198/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATATAACGAC  >  NC_002947/3965518‑3965807

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: