Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,520 T→A 5.1% intergenic (‑77/‑80) 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,5200TA 5.1% 219.0 / 2.3 78intergenic (‑77/‑80)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/34);  new base A (2/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAT  >  NC_000913/595382‑595668
                                                                                                                                          |                                                                                                                                                    
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)
                                                                                                                                          tttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCAATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAt  >  2:421529/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAT  >  NC_000913/595382‑595668

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: