Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,520 T→A 6.2% 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 6.2% 455.4 / 12.0 152intergenic (‑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 (74/68);  new base A (6/4);  total (80/72)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 157 total aligned reads displayed.

ACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAT  >  NC_000913/595372‑595668
                                                                                                                                                    |                                                                                                                                                    
aCCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                                                                     >  2:2008648/1‑150 (MQ=255)
 ccGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAAGTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAtttt                                                                                                                                                    >  1:2055114/1‑150 (MQ=255)
   ggCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGt                                                                                                                                                  >  2:704504/1‑150 (MQ=255)
        cGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                                                                             >  2:2143990/1‑150 (MQ=255)
         ggTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtttt                                                                                                                                            >  2:219393/1‑150 (MQ=255)
         ggTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtttt                                                                                                                                            >  2:1555947/1‑150 (MQ=255)
          gTTAACCCTATGGTCAAGTATTATAAGGTTTTCTTATCATCTTCGACAATCAACAGTTTCATTTTTCCTCCTCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCACTCGCTTATTGGCAAAATGACAAAATTGTCAttttg                                                                                                                                           <  1:497686/150‑2 (MQ=255)
           ttAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                                                                          <  1:255674/150‑1 (MQ=255)
           ttAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                                                                          <  1:1264008/150‑1 (MQ=255)
           ttAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgc                                                                                                                                          <  2:1708898/150‑3 (MQ=255)
           ttAACCCTTTGGTCAAGTATTATCCGGATTTCATTTCATCTTCGCCAATCCACAGTTTCATATTACCTCCGCATGTTGCCCGGGCAATTCTAGAGTCGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                                                                          <  2:346949/150‑1 (MQ=255)
            tAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCTGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgcc                                                                                                                                         <  1:108261/150‑4 (MQ=255)
             taCCCTTTTGTAAAGTATTCTGCGGTTTTCTTTTCATCTTCGACAATCAACATTATCATATTTCCTCCGCCTGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGTAAATTTTGTCATTTTTCTg                                                                                                                                        <  2:125861/149‑1 (MQ=255)
             aaCCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                                                                        <  2:946632/150‑1 (MQ=255)
                 cTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                                                                                    <  1:1715375/150‑1 (MQ=255)
                   ttGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                                                                  <  1:727647/150‑1 (MQ=255)
                   ttGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCg                                                                                                                                  <  2:1888887/150‑1 (MQ=255)
                    tGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCAGGGCAATTCTAGAGTAGAGGGATCAGATGGCAATAGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                                                                                 >  1:364969/1‑150 (MQ=255)
                      gTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaa                                                                                                                               >  2:510274/1‑150 (MQ=255)
                      gTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaa                                                                                                                               >  2:844281/1‑150 (MQ=255)
                           gTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                                                                          >  1:541246/1‑150 (MQ=255)
                           gTATTCTCCGGTTTTCTTTTCAGCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCCCCGGAAAATCa                                                                                                                          >  2:1412787/1‑150 (MQ=255)
                            tATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATAAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCag                                                                                                                         >  1:1423529/1‑150 (MQ=255)
                                    ggTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                                                                                 <  1:1025877/150‑1 (MQ=255)
                                       tttcttttctTCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTTTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAg                                                                                                              <  1:832350/150‑1 (MQ=255)
                                        ttcttttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGt                                                                                                             <  1:1483248/150‑1 (MQ=255)
                                            tttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                                                                         >  2:1314794/1‑150 (MQ=255)
                                            tttcATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                                                                         <  2:1692949/150‑1 (MQ=255)
                                               cATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCAAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTg                                                                                                      >  1:288070/1‑150 (MQ=255)
                                                  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTTGAGGGATCAGATTGCAATCGCTTATTGGCAAAATGACATTTTTGTCATTTTTCTGTCAACGGAAAATCAGAGCCTGGAGAGTAAAGTTggcg                                                                                                   >  1:1616864/1‑150 (MQ=255)
                                                  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                   >  2:869667/1‑150 (MQ=255)
                                                  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                   <  1:713609/150‑1 (MQ=255)
                                                  cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                                                                                   >  2:898338/1‑150 (MQ=255)
                                                    tCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcggc                                                                                                 <  2:723965/150‑1 (MQ=255)
                                                     cGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCa                                                                                                <  1:1989004/150‑1 (MQ=255)
                                                        cAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAATTTGGCGGCATaa                                                                                             >  2:1348769/1‑150 (MQ=255)
                                                        cAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaa                                                                                             >  2:791699/1‑150 (MQ=255)
                                                         aaTGAACTGTGTCATATATAGTCCGCATGTTGCACGGGCAATTCTCGATTAGCGGGATAAGATGGCAACCTATTATTGGCAAAATGAAAATTTTGTCATTTTTCTGTCACCTGAAAATCAGAGACTGGCGAGTAAAGTTGGCGGCATata                                                                                            <  1:292128/150‑3 (MQ=255)
                                                           tCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTATCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTATCTGTCACCGGAATATCAGAACCTTGCGAGTAAAGTTGGCGGCATAAAAt                                                                                          >  2:330202/1‑150 (MQ=255)
                                                           tCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCTTTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAAt                                                                                          <  1:1854149/150‑1 (MQ=255)
                                                            cAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATc                                                                                         >  1:1976927/1‑150 (MQ=255)
                                                               cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                                                                      <  2:2093314/150‑1 (MQ=255)
                                                                aGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCa                                                                                     <  2:550390/150‑1 (MQ=255)
                                                                 gattCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                                                                                    <  1:1739764/148‑1 (MQ=255)
                                                                 gTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                                                                                    <  1:1204953/150‑1 (MQ=255)
                                                                 gTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACAAg                                                                                    >  2:1478464/1‑150 (MQ=255)
                                                                    tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaaa                                                                                 >  2:53994/1‑150 (MQ=255)
                                                                     cATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAAt                                                                                <  2:1138554/150‑1 (MQ=255)
                                                                         tttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATg                                                                            >  1:1879155/1‑150 (MQ=255)
                                                                              tccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAACGTTGGCGGCATAAAATCACCAGAAATTATGAGCCt                                                                       >  2:428277/1‑150 (MQ=255)
                                                                              tccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATACAATCACCAGAAATTATGAGCCt                                                                       >  2:1947509/1‑150 (MQ=255)
                                                                              tccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCt                                                                       <  2:1960501/150‑1 (MQ=255)
                                                                               ccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTa                                                                      >  2:1126734/1‑150 (MQ=255)
                                                                               ccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTa                                                                      <  2:1414740/150‑1 (MQ=255)
                                                                                cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGTTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTAt                                                                     <  2:779790/150‑1 (MQ=255)
                                                                                 gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                                                                    <  1:698844/150‑1 (MQ=255)
                                                                                 gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg                                                                    <  1:2080170/150‑1 (MQ=255)
                                                                                   aTGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc                                                                  <  1:2008648/150‑1 (MQ=255)
                                                                                       tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGCCCATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTGGGCGGCATAAAATCACCAGAAATTCTGAGCCTATGTCTCCt                                                              <  1:1808374/150‑1 (MQ=255)
                                                                                       tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt                                                              <  2:813972/150‑1 (MQ=255)
                                                                                        gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCtt                                                             >  1:1443851/1‑150 (MQ=255)
                                                                                        gCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGGCATTTTTCTGTCGCCGGAAGATCAGAGCGTGGCGAGGTATGTGGGCGGCTTTAAATGCGCAGAAATTATGAGGCTATGTCGCCtg                                                             >  2:2056660/1‑149 (MQ=255)
                                                                                          ccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGt                                                           >  1:819496/1‑150 (MQ=255)
                                                                                           cGGGCTATTCTAGAGTTGCTGGATGAGATGGCAATAGCTTATTGGCAAAATGACAATTTAGTCATTTTTCTGTCACCGGATAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTa                                                          <  1:152155/150‑1 (MQ=255)
                                                                                             ggCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaaa                                                        <  2:676978/150‑1 (MQ=255)
                                                                                               cAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACt                                                      >  1:1079911/1‑150 (MQ=255)
                                                                                                   tCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAACCTTCTg                                                  >  2:1605035/1‑150 (MQ=255)
                                                                                                   tCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTg                                                  >  1:625490/1‑150 (MQ=255)
                                                                                                    cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGc                                                 >  1:211916/1‑150 (MQ=255)
                                                                                                     tAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGcc                                                >  1:652845/1‑150 (MQ=255)
                                                                                                          tAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTATGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTTTAAACTTCTGCCAtttt                                           >  2:2065789/1‑150 (MQ=255)
                                                                                                            gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTATCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACAAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGACATTTtga                                         >  1:437121/1‑149 (MQ=255)
                                                                                                               ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgtg                                      >  2:1827562/1‑150 (MQ=255)
                                                                                                               ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgtg                                      <  2:1729410/150‑1 (MQ=255)
                                                                                                               ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgtg                                      >  1:1902101/1‑150 (MQ=255)
                                                                                                               ggATCAGATGGCAATAGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTtgtgtg                                      <  1:1885638/150‑1 (MQ=255)
                                                                                                                 aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGc                                    >  2:625765/1‑150 (MQ=255)
                                                                                                                 aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGGCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGc                                    <  2:1564968/150‑1 (MQ=255)
                                                                                                                        ggCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTgcg                             <  1:1915139/150‑1 (MQ=255)
                                                                                                                         gCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCTGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTgcgc                            <  2:544698/150‑1 (MQ=255)
                                                                                                                         gCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTgcgc                            <  1:1950766/150‑1 (MQ=255)
                                                                                                                          cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCt                           >  1:1853433/1‑150 (MQ=255)
                                                                                                                              cGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAAcc                       <  1:919467/150‑1 (MQ=255)
                                                                                                                                 ttATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgt                    <  1:53994/150‑1 (MQ=255)
                                                                                                                                   aTTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttg                  >  1:1197964/1‑150 (MQ=255)
                                                                                                                                     tGGCAAAATGACAAATTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATAACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCAGTTTGTGTGGCCCTTGCGCTAACCGgttgat                >  2:1672213/1‑148 (MQ=255)
                                                                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAACCTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttgtt                >  1:1106071/5‑150 (MQ=255)
                                                                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGgttgtt                >  2:133864/5‑150 (MQ=255)
                                                                                                                                      ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTc               >  2:2191078/1‑150 (MQ=255)
                                                                                                                                      agaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACAAGAAATTATGACCCTATGTATCCTTGTAAACTTCTGCCATTTTTTGTGGCACTTGCGCTAACCGGTTGTTc               >  2:418673/4‑150 (MQ=255)
                                                                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAATGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCa              >  1:1098777/3‑150 (MQ=255)
                                                                                                                                         aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGATCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACt            <  1:219393/150‑1 (MQ=255)
                                                                                                                                         aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACt            >  1:1071405/1‑150 (MQ=255)
                                                                                                                                          aaaTGATAATTTTGTCATTTTTGTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGAAAACTTCTGCCATTTAGTGTGGCCCTTGCGCTAACCGGTTGTTCACTg           <  2:102430/150‑1 (MQ=255)
                                                                                                                                          aaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCTGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTg           <  1:869667/150‑1 (MQ=255)
                                                                                                                                          aaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTg           <  1:176521/150‑1 (MQ=255)
                                                                                                                                           aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTgg          <  2:1986301/150‑1 (MQ=255)
                                                                                                                                            aTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGc         <  1:2115464/150‑1 (MQ=255)
                                                                                                                                                   ttttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAt  >  2:1466204/1‑150 (MQ=255)
                                                                                                                                                   ttttGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAt  <  2:1265500/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAACTTCTGCCATTTTGTGTGGCCCTTGCGCTAACCGGTTGTTCACTGGCACCGGAT  >  NC_000913/595372‑595668

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