Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,995,020 T→C 6.3% intergenic (+148/+29) PP_5251 → / ← PP_5252 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,995,0200TC 6.3% 233.1 / 4.3 80intergenic (+148/+29)PP_5251/PP_5252transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base C (3/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGCTGC  >  NC_002947/5994873‑5995167
                                                                                                                                                     |                                                                                                                                                     
tGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTAtt                                                                                                                                                      >  8:273232/1‑149 (MQ=255)
     gCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAGGGGGCCAGTTATTGGccc                                                                                                                                                 <  5:373598/149‑1 (MQ=255)
     gCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGccc                                                                                                                                                 >  6:98957/1‑149 (MQ=255)
       ccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGGAGACACACCCCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGccccc                                                                                                                                               >  1:294961/1‑149 (MQ=255)
           aCGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCTGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGt                                                                                                                                             <  6:273577/147‑1 (MQ=255)
            cGGCCCTTTAGCGTGTGAACCCGCTCCTACACAGATCTTCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAATGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTc                                                                                                                                          <  8:316911/148‑1 (MQ=37)
            cGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAGAGGGCCGGCACTGAAGAGACAAGAAAGGGGCCAGGTAGGGGCCCCCGGTCg                                                                                                                                         >  6:408868/1‑149 (MQ=37)
             ggCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGt                                                                                                                                        <  7:125766/149‑1 (MQ=255)
               cccATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTc                                                                                                                                       >  4:103709/1‑148 (MQ=255)
                ccATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGACAGGGGCCAGTTATTGGCCCCCGTTCGTcac                                                                                                                                     >  5:292053/1‑149 (MQ=255)
                 cATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTcaca                                                                                                                                    >  3:321898/1‑149 (MQ=255)
                  aTTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTcaca                                                                                                                                    >  6:42296/1‑148 (MQ=255)
                   ttCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACAtc                                                                                                                                  <  6:38950/149‑1 (MQ=255)
                      gcgGGTGAACCCGCTCCTACACAGATCTGCGGTACACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGACTAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATCGGCCCCCGTTCGTCACAtcag                                                                                                                                <  1:320969/149‑1 (MQ=38)
                         ggTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccctt                                                                                                                                              <  6:62430/132‑3 (MQ=255)
                         ggTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccctt                                                                                                                                              >  5:62430/1‑130 (MQ=255)
                            gAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATcagctcagct                                                                                                                         >  3:223745/1‑149 (MQ=255)
                               cccGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAg                                                                                                                      <  3:103709/149‑1 (MQ=255)
                                   cTCCTACACAGATCTGCGGTAGACACTCACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAg                                                                                                                  <  4:321898/149‑1 (MQ=255)
                                   cTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAg                                                                                                                  <  7:173208/149‑1 (MQ=255)
                                        acacAGATCTGCGGTAGACACACACCAGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCtt                                                                                                                        >  1:4476/1‑139 (MQ=255)
                                        acacAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCtt                                                                                                                        <  2:4476/138‑1 (MQ=255)
                                         cacaGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGGCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAg                                                                                                             >  1:388789/1‑148 (MQ=255)
                                             gATCTGCGGTAGACACAAACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACAAGCGAAAGGGCCTGAACTGAAGACACAAGAAAGGGGACAGTTATTGGCCCCCGTTAGTAACATAAGATCAGATTAGAAAGCGAAGCACg                                                                                                         >  5:333015/1‑148 (MQ=255)
                                                   cGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGaa                                                                                                       <  4:223745/144‑1 (MQ=255)
                                                   cGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGcc                                                                                                  >  8:140426/1‑149 (MQ=255)
                                                     gTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCa                                                                                                <  8:51826/149‑1 (MQ=255)
                                                      tAGACACACACCGGACGTTGTAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGccc                                                                                                 <  5:98957/147‑1 (MQ=255)
                                                        gACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccct                                                                                                                                               <  6:6160/100‑2 (MQ=25)
                                                        gACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccct                                                                                                                                               >  5:6160/1‑99 (MQ=25)
                                                           acacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCAGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAc                                                                                          <  7:140426/149‑1 (MQ=255)
                                                             acacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                                        <  7:273232/149‑1 (MQ=255)
                                                                    gACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAg                                                                                  >  1:231160/1‑149 (MQ=255)
                                                                     aCGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGc                                                                                 >  7:301901/1‑149 (MQ=255)
                                                                        ttGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGGCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                              <  2:388789/149‑1 (MQ=255)
                                                                        ttGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                              <  2:294961/149‑1 (MQ=255)
                                                                         tGGAACCAGGCGGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                              <  4:48096/148‑1 (MQ=255)
                                                                          ggAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTg                                                                             >  2:353580/1‑148 (MQ=255)
                                                                             aaCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                                          <  2:186167/116‑1 (MQ=255)
                                                                             aaCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                                          >  1:186167/1‑116 (MQ=255)
                                                                             aaCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAg                                                                         <  5:42296/149‑1 (MQ=255)
                                                                              aCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAg                                                                         >  1:160014/1‑148 (MQ=255)
                                                                                  gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCAtt                                                                    >  2:18899/1‑149 (MQ=255)
                                                                                     gTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTAc                                                                 <  8:301901/149‑1 (MQ=255)
                                                                                        gTAGGAGCGGGTTTACCCGCAAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCACCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCgg                                                              >  6:117660/1‑149 (MQ=255)
                                                                                          tttcgCGGGTAAACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCa                                                           >  1:400246/5‑149 (MQ=255)
                                                                                          aGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCa                                                           >  6:11300/1‑149 (MQ=255)
                                                                                              gCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                              >  5:81558/1‑127 (MQ=255)
                                                                                              gCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                              <  6:81558/127‑1 (MQ=255)
                                                                                                  gTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAACGGGGCCAGTTATTGGCCACCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGaa                                                   >  6:409651/1‑149 (MQ=255)
                                                                                                  gTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGGTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCACTGCAAGTCATTTACCGGCACCCTcca                                                    >  6:309135/1‑146 (MQ=255)
                                                                                                      aCCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGcc                                               >  5:368264/1‑149 (MQ=255)
                                                                                                        ccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGt                                                                                                                                             >  8:237194/1‑55 (MQ=255)
                                                                                                        ccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGt                                                                                                                                             <  7:237194/55‑1 (MQ=255)
                                                                                                        ccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCcgc                                             <  8:44794/149‑1 (MQ=255)
                                                                                                        ccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCcgc                                             >  7:44794/1‑149 (MQ=255)
                                                                                                          gcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTGTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCcgcg                                            >  1:106232/1‑149 (MQ=255)
                                                                                                          gcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAACCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCgcgc                                           >  3:309176/1‑149 (MQ=255)
                                                                                                           agAAAGGTCCGGCACTGAAGACACAAGAAAGTTGCCAGTTATTGGCACCCGTTCGTAACATCAGCTCAGCTTAGAAAGCGAAGAAAGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGTCACGCTCGAATGCCgcgca                                          <  6:333015/148‑2 (MQ=255)
                                                                                                               aaGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCttt                                       >  5:412729/1‑149 (MQ=255)
                                                                                                                 gggCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTc                                     >  4:165702/1‑149 (MQ=255)
                                                                                                                     cGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTgg                                 <  6:163650/149‑1 (MQ=255)
                                                                                                                      ggCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCTACGCTCGAATGCCGCGCCTTTTCGTgg                                 <  1:18899/149‑1 (MQ=255)
                                                                                                                          cTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCtgtg                            >  4:398727/1‑149 (MQ=255)
                                                                                                                             aaGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTgc                           <  2:106232/147‑1 (MQ=255)
                                                                                                                             aaGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCAAATGCCGCGCCTTTTCGTGGCTGTgc                           <  5:117660/147‑1 (MQ=255)
                                                                                                                                  acaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCCCCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCCTGca                    >  7:373253/1‑149 (MQ=255)
                                                                                                                                  acaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGca                    >  5:41881/1‑149 (MQ=255)
                                                                                                                                    aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGc                                                        <  4:225090/111‑1 (MQ=255)
                                                                                                                                    aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCt                                                       >  5:336654/1‑112 (MQ=255)
                                                                                                                                    aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCt                                                       <  6:336654/112‑1 (MQ=255)
                                                                                                                                     aGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCacac                 >  1:280267/1‑149 (MQ=255)
                                                                                                                                       aaaGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTgc                           <  1:194786/137‑1 (MQ=255)
                                                                                                                                       aaaGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTgc                           >  2:194786/1‑137 (MQ=255)
                                                                                                                                        cgggggCCAATAACTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCAcacca               >  7:345885/3‑148 (MQ=255)
                                                                                                                                           gggCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAg           >  3:275612/1‑149 (MQ=255)
                                                                                                                                                aGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAgg          >  7:312613/1‑145 (MQ=255)
                                                                                                                                                aGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAgg          <  8:312613/145‑1 (MQ=255)
                                                                                                                                                 gTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCAGTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGc     >  3:269949/1‑149 (MQ=255)
                                                                                                                                                  ttATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGCt    <  5:408868/149‑1 (MQ=255)
                                                                                                                                                    aTTGGCCCCAGTTCGTCAAATCAGCACAGCTTAGAAAGCTAAGCACGCACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCAGCGCCTTTTCGTGGCTGTGCGCATCCACACCACAAGGCCCGCTGc  >  6:210586/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGCTGC  >  NC_002947/5994873‑5995167

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: