Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,995,018 T→A 6.3% intergenic (+146/+31) 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,0180TA 6.3% 239.0 / 5.0 79intergenic (+146/+31)PP_5251/PP_5252transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/33);  new base A (3/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGCT  >  NC_002947/5994873‑5995165
                                                                                                                                                   |                                                                                                                                                     
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                                                                            <  6:81558/127‑1 (MQ=255)
                                                                                              gCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                            >  5:81558/1‑127 (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                                                     <  6:336654/112‑1 (MQ=255)
                                                                                                                                    aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCt                                                     >  5:336654/1‑112 (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)
                                                                                                                                                   |                                                                                                                                                     
TGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACAAGGCCCGCT  >  NC_002947/5994873‑5995165

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: