Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,344,319 A→G 12.3% intergenic (+50/‑106) PP_3817 → / → PP_3818 polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,344,3190AG12.3% 156.8 / 13.3 66intergenic (+50/‑106)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base G (5/3);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

GGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGC  >  NC_002947/4344173‑4344458
                                                                                                                                                      |                                                                                                                                           
gggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGcaca                                                                                                                                           <  3:290799/149‑1 (MQ=255)
 ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGcac                                                                                                                                            >  3:382345/1‑147 (MQ=255)
 ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGcac                                                                                                                                            <  4:382345/147‑1 (MQ=255)
   cGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                        <  1:220610/149‑1 (MQ=255)
                   cGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                             <  1:275440/144‑1 (MQ=255)
                   cGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                             >  2:275440/1‑144 (MQ=255)
                    gTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCa                                                                                                                        >  3:63760/1‑148 (MQ=255)
                         gCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAAAAAAGTCTTGTACGTTGAGAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTa                                                                                                                     >  1:126506/1‑149 (MQ=255)
                               gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTgg                                                                                                                                   <  4:300332/126‑1 (MQ=25)
                               gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTgg                                                                                                                                   >  3:300332/1‑126 (MQ=25)
                                 aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                 >  4:147445/1‑126 (MQ=255)
                                 aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                 <  3:147445/126‑1 (MQ=255)
                                 aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTa                                                                                                          >  8:202385/1‑149 (MQ=255)
                                     gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcc                                                                                                      >  3:374105/1‑149 (MQ=255)
                                       ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCATCCAGATAGAGCACCCCCCGCGCGGCCCACGATATGGCAGGCGCTGCACTATTTGa                                                                                                                >  1:412031/1‑137 (MQ=255)
                                            tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                               <  7:202385/149‑1 (MQ=255)
                                              tataCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcaca                                                                                              <  4:374105/148‑1 (MQ=255)
                                                  ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                         <  3:147506/149‑1 (MQ=255)
                                                      attggCCGAGAGCGACCTGAGCAGGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGTAGGCGCTGCACTATTggt                                                                                                                <  3:257334/119‑2 (MQ=255)
                                                      aTGGGCCGAGAGCGACCTGAGCAGGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGCCACACGATCTGGTAGACGCTGCACTATTGGa                                                                                                                >  4:257334/1‑122 (MQ=255)
                                                         ggCCGATAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcgc                                                                                                                            <  7:57650/107‑1 (MQ=255)
                                                                agCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                           >  6:124847/1‑149 (MQ=255)
                                                                     ccTGAGCAAGTCGGCCGGGTCGCATA‑‑ACCGCTCTTTTACGTTCCAGCGCCTGCCAGATCGCCCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACACCCCGTCACATTATCATCACGCCACTGTCGTa                                                                        <  4:63760/147‑1 (MQ=255)
                                                                     ccTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                                       >  4:293629/1‑148 (MQ=255)
                                                                     ccTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCAAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                                       >  7:356677/1‑148 (MQ=255)
                                                                        gAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCATGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACAtt                                                                                            >  1:90157/1‑125 (MQ=255)
                                                                        gAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACAtt                                                                                            <  2:90157/124‑1 (MQ=255)
                                                                           cAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGACCACCGCCCGCGCGGCACACGAACTGGCAGGCGCAGCACTATTGGAGGCGTACGCCCCGTGACATTATCATCACGCCACAGTCGTAAACTTCCg                                                                >  3:75313/1‑149 (MQ=255)
                                                                              gTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGAAGGAGCACCCCCCGCGCGGCACACGCTCTGGCAGGCGCTGCACAATTGGCGGCGTACGCCCCGGCACATTTTCATCACCCCACAGTCCTAAACTTCCGGc                                                              >  1:359128/1‑148 (MQ=255)
                                                                                  gCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                         >  7:113289/1‑149 (MQ=255)
                                                                                    cGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                       <  2:126506/149‑1 (MQ=255)
                                                                                          gcATT‑‑ACCGCTCTTGTACGTTGCAGCGCCAGCCAGATCGACCACCCCCCGCGCGGCACACGCACTTGCAGGCGCTGCATTAATGGCGGCGTACCCCCCGTCACATTATACTCCAGCCACAGTCGTAATCTGCCGGCCAACTCAATGccc                                                 >  6:235833/1‑148 (MQ=255)
                                                                                          gcATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCAAGATCGAGCACGGCCCGCGGGGCAAACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                                 <  2:359128/149‑1 (MQ=255)
                                                                                            aTT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                                 >  2:94740/1‑147 (MQ=255)
                                                                                            aTT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                                 >  1:235091/1‑147 (MQ=255)
                                                                                                 aaCGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGt                                            >  3:64808/1‑149 (MQ=255)
                                                                                                   cGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGccc                                                                                                     <  1:68792/90‑1 (MQ=255)
                                                                                                   cGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGccc                                                                                                     >  2:68792/1‑90 (MQ=255)
                                                                                                    gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                                 <  5:420425/125‑1 (MQ=255)
                                                                                                    gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                                 >  6:420425/1‑125 (MQ=255)
                                                                                                    gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACAAGATCTGGCAGGCGCTGCACTATTGGAGGCCTACGCCCCGTCACATTATCATCACCCCACTGTAGTAATCTGACGGCCAACTCTATTCCACAGGTCAc                                         >  6:234977/1‑149 (MQ=255)
                                                                                                      tctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTc                                           <  8:37417/145‑1 (MQ=255)
                                                                                                      tctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTc                                           >  7:37417/1‑145 (MQ=255)
                                                                                                          gTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGc                                    <  5:124847/148‑1 (MQ=255)
                                                                                                              gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCt                               >  7:199024/1‑149 (MQ=255)
                                                                                                                 gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTcc                             <  8:113289/148‑1 (MQ=255)
                                                                                                                 gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCg                            >  5:414101/1‑149 (MQ=255)
                                                                                                                   cAGCGCCGGCACGATCGCCCCCGGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCg                            <  6:234552/148‑1 (MQ=255)
                                                                                                                    aGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                          >  7:68969/1‑149 (MQ=255)
                                                                                                                    aGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGGCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                          >  7:168731/1‑149 (MQ=255)
                                                                                                                         cTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                     <  6:414101/149‑1 (MQ=255)
                                                                                                                           tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTg                    >  5:138762/1‑149 (MQ=37)
                                                                                                                           tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcaca                                                                                              >  8:134687/1‑75 (MQ=255)
                                                                                                                           tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcaca                                                                                              <  7:134687/75‑1 (MQ=255)
                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc                  >  8:3212/1‑149 (MQ=255)
                                                                                                                              cAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                      >  4:396210/1‑112 (MQ=25)
                                                                                                                              cAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                      <  3:396210/112‑1 (MQ=25)
                                                                                                                               aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                <  8:68969/149‑1 (MQ=255)
                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                                    <  7:23628/111‑1 (MQ=255)
                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                                    >  8:23628/1‑111 (MQ=255)
                                                                                                                                   cgtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg            >  1:167545/7‑149 (MQ=38)
                                                                                                                                    gtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTcc                             >  4:245982/6‑131 (MQ=38)
                                                                                                                                    gtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTcc                             <  3:245982/126‑1 (MQ=38)
                                                                                                                                     aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct          >  4:144478/1‑149 (MQ=255)
                                                                                                                                          cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                       >  4:185901/1‑99 (MQ=255)
                                                                                                                                          cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                       <  3:185901/99‑1 (MQ=255)
                                                                                                                                          cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAccc     >  8:127382/1‑149 (MQ=255)
                                                                                                                                           gcccccGCGGCACACGATATGTCTGGCGCTTCAATATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCt    <  5:235833/149‑1 (MQ=255)
                                                                                                                                             ccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGc  <  2:136644/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
GGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATT‑‑AACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGC  >  NC_002947/4344173‑4344458

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