Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,622,2500AG36.0% 77.1 / 13.2 51E161G (GAG→GGG) yjjJtoxin YjjJ
Reads supporting (aligned to +/- strand):  ref base A (4/28);  new base G (18/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.05e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATT  >  NC_000913/4622011‑4622497
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
gtgCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGAtttat                                                                                                                                                                                                                                                      >  1:444228/1‑248 (MQ=255)
  gCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGTAGGGGAATTATCAg                                                                                                                                                                                                                                                    <  2:670440/251‑1 (MQ=255)
  gCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTGGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCTGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATTATATACGTCTCGGTCAGGAAGAAGATGTGCTTTACTCGCTGACCGTATTTATCGGTGAATATACTGGCTGTTGGTGGGTTGGGGAGGGGAAGTATCa                                                                                                                                                                                                                                                     <  2:348261/250‑1 (MQ=255)
       ccGCAGGGAAGTTGTCTGGTTACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGACGATGTGCTCTACGCCCTGACCGTATTTAACGGAGAATATACTGGCGGTGGGTTGGTTGGGGTGGGGAGTTATAAGGGATg                                                                                                                                                                                                                                               >  2:519566/1‑251 (MQ=255)
         gCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTTGTCGGGGGGGGGGAtttttcaggatgtt                                                                                                                                                                                                                                             >  1:120129/1‑239 (MQ=255)
               aaGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTGTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTTCTCTACGCCCTGACCGTATTTAACGGTGAATATTCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGTATTACtc                                                                                                                                                                                                                                      <  2:42471/251‑2 (MQ=255)
                    gTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGcaca                                                                                                                                                                                                                                  <  1:513581/250‑1 (MQ=255)
                      ccggTAACAGGCGGAGATGGCGTCGAACAGTGGTTTGATGGTTTGCCCTGGTAGTTGACCGATATCCGAACGTAGGGCTTGTTAGGGGGGGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTGTGGCAGGAAGAAGATGTGCGCTGCGGCCGGACCGTAGTTAACGGTGAAGATACTGGCGGTTGGTTGGGCGGGGAGGGGAATTATCAGCGATGGATTTCTGCacaac                                                                                                                                                                                                                                <  2:329528/248‑1 (MQ=255)
                        gggaaCAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGGGAATATACTGGCGGTTGGTGGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaac                                                                                                                                                                                                                                <  2:644905/245‑1 (MQ=255)
                        ggTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTGTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACAcc                                                                                                                                                                                                                             <  2:429188/251‑1 (MQ=255)
                            aCAGGCGCAGATGGTGGCGAACAGTGGTTTGGTGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGc                                                                                                                                                                                                                         <  2:653167/251‑1 (MQ=255)
                              aGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGTCGGGGAGGGGGGTTATCAGCGATGGGTTACTTCCCAACACCCTGCgg                                                                                                                                                                                                                       >  1:608646/1‑251 (MQ=255)
                                gcgcAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTTTCAGCGATGGATTACTGCacaaacacctgcggg                                                                                                                                                                                                                      >  1:322559/1‑239 (MQ=255)
                                   cAGATGGCGACGAACAGTGGGTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGGCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGAGGGATTACTGCGCAACGCCCTGCGGAAAtt                                                                                                                                                                                                                 <  2:131903/251‑1 (MQ=255)
                                    aGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGAGCGTATTTAACGGTGAATATACTGGCGGGTGGGTGGTCGGGGAGGGGAATGATCAGCGATGGATTACTGCacaaca                                                                                                                                                                                                                               <  2:414459/237‑1 (MQ=255)
                                    aGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGGTTATCAGCGATGGATTACTGCacaacaagatcggaaaagg                                                                                                                                                                                                                  >  1:414459/1‑237 (MQ=255)
                                     gATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGGCGGGGGGGGGGATTTTTTGCGAGGGATTACTTCacaaacccccgcggagattac                                                                                                                                                                                                                >  2:101780/1‑234 (MQ=255)
                                        ggCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGGCCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCctct                                                                                                                                                                                                            <  2:395081/251‑1 (MQ=255)
                                        ggCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTTGGGGGGGGGAAATATTAACGATGGATTACTGCACAAAACCCTGCGGAAATTTctct                                                                                                                                                                                                            >  1:513789/1‑251 (MQ=255)
                                                              gTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGAATATCAGCGAAGGGGTATTGCACAACCACCTTCGGGAATACTCCTGGGTCAAAAAATCAccccatac                                                                                                                                                                                      >  2:49756/1‑246 (MQ=255)
                                                                         tATTTGACCGATCTCCGACCGCAGGGCTTTTTAGAGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGTGCGGGGGGGGGTATTTTCTACGATGGAGTTATGCCCCACACCCTGCCGAAAATACTCTGGATCAAAAAATCAcccccttccacacaccggg                                                                                                                                                                           >  2:542166/1‑235 (MQ=255)
                                                                                        cGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGc                                                                                                                                                            <  1:4608/251‑1 (MQ=255)
                                                                                              cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGCCCGTATTTAACGGTGAATATACTGGCGGCTGGTTGGGCGGGGGGGGGAATTATCCGCGATGGGTTTCTGCCCAAAACCCTGCGGGAAATCCTCTTGATCAAAAACTCACCCACTACCAAAAGCCTGCAAGTGGTGCAATGGCAgggg                                                                                                                                                       >  1:289928/1‑248 (MQ=255)
                                                                                              cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAAACCCCTGCGGGAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAggggg                                                                                                                                                      >  1:322113/1‑248 (MQ=255)
                                                                                                      tttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGCTGGGTTTCTGCCCAAACCCCCGCGGAAATTTCTCTTGATCCAAAAATCACCCATTACGAACCGCTGGCAAGGGATGCAATGGCAGGAGGAATTTTgg                                                                                                                                           >  1:36607/1‑251 (MQ=255)
                                                                                                          gggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTAGTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCGGCGGAAATTCCTCTGGATCAAAAACTCACCCAGTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGttc                                                                                                                                       <  1:75377/251‑1 (MQ=255)
                                                                                                           ggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGCTGGTCGGGGGGGGGGATTATCAGAGATGGATTAATGCCCCACACCCTGCGGAAATTCCTCTGGGTCAAAAACTCACCCCTTACGAACAGCTGGCAAGTGATTCACTGGCAGGTGAAATTGTGGGttct                                                                                                                                      >  1:82024/1‑251 (MQ=255)
                                                                                                              gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTtcttct                                                                                                                                   <  1:613376/251‑1 (MQ=255)
                                                                                                              gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTtcttct                                                                                                                                   <  1:689867/251‑1 (MQ=255)
                                                                                                                                     gcgcTACTGAAGTTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTTCGCCATGTCCGTATTTAACGGTGAATATACTGGCGGTTGGTTGATCGGGGAGGGGAATTATCAGCGATGTATTACTGCACAACACCCTGCTGAGATTCCTCTGGATCAAAAACTCACCCATTGCGAACAGCTGGCAAGTGATGCACTGGCAGGATACATTTTGGGTTCTTCTGCGGGCGGCGAGCAGCCaaaacc                                                                                                           <  2:20977/251‑3 (MQ=255)
                                                                                                                                          acggAATCTGACTGTGGATATCCGTTTTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGTTTCTTCTGCGGGCGGCGAGCAGCCAAAAGTTACcta                                                                                                    <  2:459075/248‑1 (MQ=255)
                                                                                                                                            tGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATCAGCGATGGATTAATGCACAACACCCTGCGGGAATTCCTCTGGAACAAAAAATCCCCCATTACGAACAGCTGGCAAGTGATGCACCTGCCGGAGAAATTGTGGGGTCTTCTGCGGGGGGGGGGCAGCCAAAAATTACctac                                                                                                   >  1:477323/1‑251 (MQ=255)
                                                                                                                                              aaTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGGTGGGTTACTGCACAACACCCTGCGGGAATTTCTCTGGATCAAAAACTCACCCACTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGGAAATTTGGGTTCTTCTTCGGGCGGCGGGCAGCCAAAATTTACctaat                                                                                                  >  1:30738/1‑248 (MQ=255)
                                                                                                                                               aTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTTGTTGGTCGGGGGGGGGGATTATCAGGGGTGGGTTTCTGCACAACACCCTGCGGAAAATTCCCTGGATCAAAAAAACCCCCACTACACAAAGATGGCAAGGGAGACACTGGCAGGGAGAAATGTGGGGTCTTTTGCGGGGGGGGAGGAGACCAAATTTAcaatctat                                                                                                >  1:523326/1‑244 (MQ=255)
                                                                                                                                                   gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATg                                                                                               <  2:30738/248‑1 (MQ=255)
                                                                                                                                                    actgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTCTTTAACGGTGAATATACTGGCGGATGGTTGGTCGGGGGGGGGAATTAACAGCGATCGATTACTGCACAACACCCTGCGGGAATTCCTCTGGATCAAAAACCCACCCATTACGAACAGCTGGCAAGTGATGCACCGGCAGGAGGAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATTcac                                                                                            >  1:108616/1‑250 (MQ=255)
                                                                                                                                                                cGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTAGGCACAAACGCCGTCAgg                                                                               <  1:519566/251‑1 (MQ=255)
                                                                                                                                                                     ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAGTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGTCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCa                                                                             <  1:590129/248‑1 (MQ=255)
                                                                                                                                                                       ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGGTCTGCTGCTGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAATCGCCGTCAGGCGATaaa                                                                        <  2:305345/251‑1 (MQ=255)
                                                                                                                                                                               gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGCTGGGCGGGGGGGGGAATTATCAGCGATGGTTTACTGCCAAACACCCCGCGGGAATTCCTCTGGATCAAAAAATCAACCATTACGAACAGCTGGCAAGGGATGCACTCGCAGGAAGAAATTTGGGTTCTTCTGCGGGCGGCGAGCAGACCAAATTTACCTACTATGGACAAACGACGCCACGCAAAAAACAtaggtt                                                                >  2:578556/1‑246 (MQ=255)
                                                                                                                                                                                            aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCa                                                    <  2:7234/249‑1 (MQ=255)
                                                                                                                                                                                             cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGTCGGGGAGGGGTAATATCCGCGATGGATTACTGGACAACACCCTGCGGGAAATCCTCTGGGTTCAAAACTCCCCCCTTACGAACAGCCGGCAAGAGATGCACTGGCAGGAGAAAATGTGGGTTTTTCTTCGGGGGGGGGGGCGCCAAAAGTTTACTACTATTCCCCAACGACGTCAGGCAATAAACCTGTTTTGGGGAAAttacccc                                                 >  2:525681/1‑246 (MQ=255)
                                                                                                                                                                                             cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTTTTAGCGGTGGATTAATGCACAACACCCCGCGGAAATTTCTCTGGGTCAAAAAATCACCCATTACGAACAGCCGGCCAGAGATTCACTGGGAGGGGAAAATTTGGGGTTTTCTTCGGGCGGCGAGCCAACAAAAATTTCCCACTATGAAAAAACGCCGGAAGGCAATAAACATGTGTTTGTGGAAATTCAcc                                                  >  1:112726/1‑251 (MQ=255)
                                                                                                                                                                                                   gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACcac                                           <  1:467288/251‑1 (MQ=255)
                                                                                                                                                                                                           ttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGAATATCAACGATTGCTTACTGCCCAACAACCCGCGGAAAATCCTCTGGAACCAAAACCCACCCAATACCGACAGCTGGCAAGGGATGCACCGGAAGGAGGAAATTTGGGGTTTTCTGCGGGCGGCGAGGCGGCAAAAATTACCTACCTTTCACAAAAGCCATCAAGCAATAAACAATTGTTGGTGGAATTCACCggacacccccaaacc                                   >  1:185134/1‑237 (MQ=255)
                                                                                                                                                                                                                         gcTGGCGGTGGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCGGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAAGTCACCGTACCACAGCAAACCGCGGTCAGCCAACg                     <  2:120129/250‑1 (MQ=255)
                                                                                                                                                                                                                          cTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGCGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCaa                                                                            >  1:367286/1‑196 (MQ=255)
                                                                                                                                                                                                                          cTGGCGGTTGGTTGGGCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAGTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGAGGCACTGGCAGTAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTAGGCACAAACGCCGTCAGGCaa                                                                            <  2:367286/196‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGt                    <  1:392595/249‑1 (MQ=255)
                                                                                                                                                                                                                                gttggttgGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGt              <  1:49756/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATGCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCGTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGc       <  1:364119/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAAtt  <  1:251009/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
GTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATT  >  NC_000913/4622011‑4622497

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