Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,519 T→A 6.0% intergenic (‑76/‑81) cusR ← / → cusC response regulator in two‑component regulatory system with CusS/copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913595,5190TA 6.0% 215.5 / 2.8 78intergenic (‑76/‑81)cusR/cusCresponse regulator in two‑component regulatory system with CusS/copper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base T (40/33);  new base A (2/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGA  >  NC_000913/595382‑595667
                                                                                                                                         |                                                                                                                                                    
gTTAACCCTTTGGTCAAGTATTGTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGAGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtttt                                                                                                                                           >  1:62156/1‑149 (MQ=255)
 ttAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATAAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttg                                                                                                                                          <  1:706394/149‑2 (MQ=255)
  tAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                                                                         <  1:493598/149‑3 (MQ=255)
   aaCCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                                                                        <  1:1052992/149‑1 (MQ=255)
   aaCCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                                                                        <  1:698831/149‑1 (MQ=255)
             tCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaa                                                                                                                              >  1:1129267/1‑149 (MQ=255)
              cAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaa                                                                                                                             >  1:858393/1‑149 (MQ=255)
                 gTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCATTCGCTTATTGGCAAAATGAGAATTTTGTGATTTTGCTGTCACCGGAAAATc                                                                                                                          >  2:216107/1‑149 (MQ=255)
                 gTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATc                                                                                                                          <  1:749523/149‑1 (MQ=255)
                 gTATTCACCGGTTTTCTTTTCATCTTCCACAATCAACAGTTTCATATTACCTCCGCAGGTTGCCAGGGAAAGTCTAGAGTAGCGGGATCAGATGGCAATCGCATAATGGCAAAATGAAAATTTTGTCATTTTTCTATCACCCGAAAAta                                                                                                                          >  2:962876/1‑148 (MQ=255)
                     tctcCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACAGGAAAATCagag                                                                                                                      >  1:62137/1‑149 (MQ=255)
                       tcCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCTGGAGCAGATGGGAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGcc                                                                                                                    >  1:1116248/1‑149 (MQ=255)
                       tcCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCagagac                                                                                                                    >  2:168657/1‑147 (MQ=255)
                        ccGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATCGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCt                                                                                                                   <  1:258302/149‑1 (MQ=255)
                         cGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTg                                                                                                                  <  2:764338/149‑1 (MQ=255)
                          ggTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTgg                                                                                                                 >  1:727980/1‑149 (MQ=255)
                          ggTGTTCTTTTCATCTTCGACATTCAACAGTTTCATATTTCCTCCGCCTGTTGCCCGGGCAATTCTAGAGTAGCGGGATGAGATGGCAATCGCTTATTGGCAATATGACATTTTTGTCATTTTTCTGTCACCGGTAAATCAGAGCCTgg                                                                                                                 <  1:772123/149‑1 (MQ=255)
                           gTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                                                                                >  2:57458/1‑149 (MQ=255)
                           gTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                                                                                <  2:41596/149‑1 (MQ=255)
                                cttttcATCTTCGACAATCACCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTa                                                                                                           >  1:729685/1‑149 (MQ=255)
                                cttttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTa                                                                                                           >  2:527040/1‑149 (MQ=255)
                                  tttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGCTGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTaaa                                                                                                         >  1:894182/1‑149 (MQ=255)
                                   ttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTaaat                                                                                                        >  1:237511/1‑148 (MQ=255)
                                   ttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                                                                        >  1:851421/1‑149 (MQ=255)
                                    tcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTGGCCCGGGCAATTCTAGAGTAGCCGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                                                                       <  2:382413/149‑1 (MQ=255)
                                     cATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGtt                                                                                                      >  2:1032480/1‑149 (MQ=255)
                                         ttCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                  <  1:582070/149‑1 (MQ=255)
                                         ttCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                  <  2:41342/149‑1 (MQ=255)
                                         ttCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                  >  2:509930/1‑149 (MQ=255)
                                             aCAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCAt                                                                                              <  2:983662/149‑1 (MQ=255)
                                              cAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATa                                                                                             >  2:456936/1‑149 (MQ=255)
                                                  cAACAGTTTCATATTTCCTCAGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                                                                         <  1:665300/149‑1 (MQ=255)
                                                   aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCTGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATc                                                                                        <  2:905466/149‑1 (MQ=255)
                                                    aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCa                                                                                       >  2:502010/1‑149 (MQ=255)
                                                        tttCATATTTCCTCCGCATGTTGCCCCGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCat                                                                                   >  2:347355/1‑148 (MQ=255)
                                                              atTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTAATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGAGGCATAAAATCACAAGAAATTa                                                                             >  1:775913/1‑149 (MQ=255)
                                                                tgccgacGCATGGTGCCCGGGCAATTCTAGAGTAGCGGGATGAGAGGGCAATCGCTGATTGGCAAACTGACAAATATGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATg                                                                           <  2:786652/143‑1 (MQ=255)
                                                                  cctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAg                                                                         <  2:603714/149‑1 (MQ=255)
                                                                    tccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATAACAATTTTGTCATTTTTCTGTCACAGGAAAATCATAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACAAGAAATTATGAGcc                                                                       >  2:1083125/1‑149 (MQ=255)
                                                                       gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTAt                                                                    >  1:1127400/1‑149 (MQ=255)
                                                                       gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCCAGTAAAGTTGGCGGGATAAAATCACCTGAAATTATGAGCCTAt                                                                    >  1:790793/1‑149 (MQ=255)
                                                                         aTGTTGCCCGGGCAATTCTAGAGTCGCGGGATCAGATCGCAATCGCTTTATGGCAAAATGACAATTTTGGCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt                                                                  <  1:574916/149‑1 (MQ=255)
                                                                              gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt                                                             >  1:955164/1‑149 (MQ=255)
                                                                              gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt                                                             >  1:454869/1‑149 (MQ=255)
                                                                              gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTATCCt                                                             >  1:916362/1‑149 (MQ=255)
                                                                               cccGTGCAATTCTAGAGTAGCGGGATCAGATGGCAATCTCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCtt                                                            <  1:214823/149‑1 (MQ=255)
                                                                                ccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg                                                           >  1:368847/1‑149 (MQ=255)
                                                                                      aaTTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACt                                                     >  2:1156057/1‑149 (MQ=255)
                                                                                          ttatagTAGCGTGATCTGATTTCAATCGCTTTTTGTCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGCGCCTGGCGAGAAAAGCTGTCGGCATAAAATCACCAGAAATTATGACCCTATGTCTCCTTGTTAACTTCTg                                                 <  2:1167128/145‑1 (MQ=255)
                                                                                           aagagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCTTTTTTCTGTCACGGGAATATCAGAGCCTGGCGAGTAAAGTTGGCGGCATTATATCACCAGAAATTATGAGCCTTTGTCTCCTTGTAAACTTCTGc                                                <  2:190919/148‑1 (MQ=255)
                                                                                              agTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCAt                                             <  2:622162/149‑1 (MQ=255)
                                                                                              agTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCAt                                             <  1:57458/149‑1 (MQ=255)
                                                                                                  gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACATTTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGAAAACTTCTGCCATTTtg                                         >  1:1168587/1‑149 (MQ=255)
                                                                                                   cGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgt                                        <  1:1002268/149‑1 (MQ=255)
                                                                                                   cGGGATCAGATGGCAATCGATTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgt                                        <  1:59867/149‑1 (MQ=255)
                                                                                                     gtaTCATATGGTAATCGCTTATTTGCAAAATTTCAATTTTTTAATTTTTCTGTCACCGGATAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgt                                      <  1:539569/147‑1 (MQ=255)
                                                                                                     ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgt                                      <  2:851421/149‑1 (MQ=255)
                                                                                                         cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGcc                                  <  1:527040/149‑1 (MQ=255)
                                                                                                           gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCt                                <  2:916362/149‑1 (MQ=255)
                                                                                                             tGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTg                              >  2:84448/1‑149 (MQ=255)
                                                                                                              ggCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGACTCCTTGTAAACTTCTGCAATTTTGTGTGGCCCTTgc                             >  1:718130/1‑149 (MQ=255)
                                                                                                                cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCATTgcgc                           >  2:939155/1‑149 (MQ=255)
                                                                                                                 aaTCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCt                          <  1:347355/149‑1 (MQ=255)
                                                                                                                       ttATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCgg                    <  1:509930/149‑1 (MQ=255)
                                                                                                                         ataggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCCCCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgtt                  <  2:62156/146‑1 (MQ=255)
                                                                                                                         aTTGGCTAAATGACAATTTTGTCATTTTTCTGTCTCCGGAAAATCAGTGCCTGGCGAGTAATGTTGGCGGCATAAAATGACCAGATATTATGAGCCTTTGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgtt                  <  1:962876/149‑1 (MQ=255)
                                                                                                                          ttGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttg                 >  2:7174/1‑149 (MQ=255)
                                                                                                                           tGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCTTAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttgt                <  2:894182/149‑1 (MQ=255)
                                                                                                                            agaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttgtt               >  2:139351/4‑149 (MQ=255)
                                                                                                                             gCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTTAAATTCTGCCATTTTGTGTGGCCCTTGCGCTAACCTGTTTTTc              >  2:263910/1‑149 (MQ=255)
                                                                                                                             gCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTc              >  1:1157382/1‑149 (MQ=255)
                                                                                                                              aaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGACATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCAc            >  2:677392/2‑149 (MQ=255)
                                                                                                                                 aaTGACAATTTTGTCATTTTTCTGTCACAGGAAAATCAGAGCGTGGCGAGTAAAGTTGGCGGCATAAAATCCCCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTg          <  2:397823/149‑1 (MQ=255)
                                                                                                                                    gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCa       <  1:636304/149‑1 (MQ=255)
                                                                                                                                    gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATATTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCa       >  2:80318/1‑149 (MQ=255)
                                                                                                                                    gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGACTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCa       >  2:83969/1‑149 (MQ=255)
                                                                                                                                      cAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGTCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTATACTTCTGCCAGTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCCCTGGCAcc     >  1:245345/1‑149 (MQ=255)
                                                                                                                                       aaTTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCCGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGACATTTTGGTTGGCCCTTGCGCTAACCGGTTGTGCACTGGCACCg    >  1:1134195/1‑149 (MQ=255)
                                                                                                                                         ttttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGa  <  2:858393/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
GTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGA  >  NC_000913/595382‑595667

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: