Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 595,953 C→T T118T (ACC→ACT)  cusC → copper/silver export system outer membrane channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913595,9530CT96.1% 147.5 / ‑1.8 51T118T (ACC→ACT) cusCcopper/silver export system outer membrane channel
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/31);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

AACGGCCTGGTTAACGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTG  >  NC_000913/595717‑596197
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
aaCGGCCTGGTTAACGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCGTGTGAAGTCGCTGATTAGCGAGGCGCTGGTGAATAATCGTGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGTGAATACTGCCACTACACGGGa                                                                                                                                                                                                                                        <  2:300327/251‑1 (MQ=255)
 aCGGCCTGGTTAACGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAg                                                                                                                                                                                                                                       <  1:259463/251‑1 (MQ=255)
    gCCTGGTTATCGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCCGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGttt                                                                                                                                                                                                                                    >  1:181271/1‑251 (MQ=255)
      cTGGTTAACGCAGCGGATAACTATCAGAACGCGGGCTGGCGCACATTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGCTGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGAGGATACTGCCACGACACGGGGGTTTTc                                                                                                                                                                                                                                  >  2:288749/1‑251 (MQ=255)
        ggTTAACGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGGGCGGAAAACTTCAAGTGAATACTGCCACGACACGGGGGTTTTCg                                                                                                                                                                                                                                 >  2:275187/1‑250 (MQ=255)
               gcagcaGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCt                                                                                                                                                                                                                         <  1:224675/251‑1 (MQ=255)
                 agcagATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGGTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAAGCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCGCAATGGCGAGGGCAGCGGCAGCTGGAGCGTTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACTGGCCTTa                                                                                                                                                                                                                       <  2:46394/251‑1 (MQ=255)
                  gcagATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGTAGCGGTAGTCTTAAAGTGAATAATGCCGCGACACGGGAGTTTTCGGCCGGCCTTat                                                                                                                                                                                                                      >  1:16278/1‑250 (MQ=255)
                      aTAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCTCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCTATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGcc                                                                                                                                                                                                                  <  2:156651/251‑1 (MQ=255)
                             tCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTg                                                                                                                                                                                                           >  1:385322/1‑251 (MQ=255)
                              cAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTg                                                                                                                                                                                                           <  1:359491/250‑1 (MQ=255)
                                 aactcgGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCt                                                                                                                                                                                                       >  1:63076/5‑251 (MQ=255)
                                   cgcgGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCg                                                                                                                                                                                                     >  1:210192/1‑251 (MQ=255)
                                     cgGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGGGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGAt                                                                                                                                                                                                   >  2:301271/1‑251 (MQ=255)
                                          tggagcACCTTTTTTGTTGATAATCAGGTGACGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATATTAAAGGGAAAACTGCCACGACACGGGGGTTTTCGACCGGCCTTAACGACAGCTTTGACCTCGGtttttt                                                                                                                                                                                              >  2:300831/5‑251 (MQ=255)
                                                cccTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCg                                                                                                                                                                                        <  2:127965/250‑1 (MQ=255)
                                                                agcAGGTGAAGACGCTGATTAGCGAGGGGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGc                                                                                                                                                                        <  2:324493/249‑1 (MQ=255)
                                                                     gTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGc                                                                                                                                                                   <  2:210192/251‑1 (MQ=255)
                                                                           aCGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGGGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCGTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCg                                                                                                                                                             <  2:115569/251‑1 (MQ=255)
                                                                            cGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTc                                                                                                                                                                                             <  2:75532/218‑1 (MQ=255)
                                                                            cGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTc                                                                                                                                                                                             >  1:75532/1‑218 (MQ=255)
                                                                                   tAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAatt                                                                                                                                                     >  1:164205/1‑251 (MQ=255)
                                                                                                      aatCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGc                                                                                                                                  <  1:297635/251‑1 (MQ=255)
                                                                                                         cGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTc                                                                                                                                <  1:219301/250‑1 (MQ=255)
                                                                                                          gggAATTGCGCATGGCGACGCTGAAAGTGCAGGAGGAACGGGCGCAATATTGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGTGCGGGAATCTTAAAGGGAATCCTGCCACGACACGGGAGTTTTCGCCCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAACAAACATGAGCGCAGCCGACCGACAAATTTTTTTAGACACTGCGGAAGCTc                                                                                                                                >  2:238886/1‑249 (MQ=255)
                                                                                                                    cATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGc                                                                                                                    <  1:313214/251‑1 (MQ=255)
                                                                                                                       ggctaCGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTGTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCTGTGCata                                                                                                                 <  2:209521/247‑1 (MQ=255)
                                                                                                                               tGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTctgctg                                                                                                         >  1:401832/1‑251 (MQ=255)
                                                                                                                                gAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTgg                                                                                                        <  1:301271/251‑1 (MQ=255)
                                                                                                                                gAAAGGGCAGGAGGCCCGGGCGCAATTGCGTCGGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCGGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCGCTGAGGAAGCTCAGCGCGCGGTGCATATTctgctg                                                                                                         <  2:181271/250‑1 (MQ=255)
                                                                                                                                                cGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCa                                                                                        <  1:259268/251‑1 (MQ=255)
                                                                                                                                                          atCGTCTGACCGATGCCGACCGCTCCCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTAttt                                                                               <  2:368716/250‑1 (MQ=255)
                                                                                                                                                             gTCTGACCGATGACGACCGCTCCACACAGCTCACTGGCGAGGGCAGCGGCAGCTGGAGCGGTCATCTTAGAGGGAATACTGCCAAGACACGGGAGTTTTCGGCCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGCAGCAGAGCGACAAAATTATTTAGCCACTGAGGAAGGTCAGCGCGCGGGGCATGTTTTGCTGGTTTCTACTGTCGCGCAAAGCTATTTCatt                                                                           >  2:244134/1‑249 (MQ=255)
                                                                                                                                                                 gACCGATGCCGACATCTACACACAGCTCAATGGCGAGGGCAGCGTCAGCGTGAGCGGTAATCTTAAAGGGAATACTGCCACGTCCCGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTTACCTCGATTTTTACGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTTCATAGTCTGCTGGTTTCTAATGTCGCGTAAAGCTATTTCAATCAGc                                                                       <  2:225718/251‑1 (MQ=255)
                                                                                                                                                                     gATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACt                                                                   >  1:61228/1‑251 (MQ=255)
                                                                                                                                                                           gACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTa                                                             <  2:61228/251‑1 (MQ=255)
                                                                                                                                                                               gCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATgcg                                                         <  1:100445/251‑1 (MQ=255)
                                                                                                                                                                                  aCCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGTAGTTTTCGACCGGCCTTAACGCAAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCa                                                       >  1:336287/1‑250 (MQ=255)
                                                                                                                                                                                    ccccagatcaagggggAGGGCAGCGGCGGCTGGGGCGGTAATCTTAAAGGGAATACCGCCACGACACGGTAGTTTTCGACCGGCCTTAACGCAAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAAtt                                                    <  2:336287/236‑1 (MQ=255)
                                                                                                                                                                                       caGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTGAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGCACATGAGCGAAGCCGAGAGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCACTCAGAAACTGGCGAATACGCAATTGc                                                  >  2:286694/1‑250 (MQ=255)
                                                                                                                                                                                                   gAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTGGACCTCGATTTTTTCGGTCGCTTAAAGAACAGGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCgaaga                                     <  1:244134/251‑1 (MQ=255)
                                                                                                                                                                                                     gggCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACTGTCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCCCTTAAAGATCATGAGCGAAGCCGAGCGACAAAATTAGTTTGCCACTGAGGTAGCTCAGCGCGCGGTGCAGATTCTGCTGGTTTCTAATTTCGCGCAAATCTATTTCAATCAGCAACTGTCGTATGCGCAATTGCAAATAGCCGAAGaaa                                   <  1:405475/251‑1 (MQ=255)
                                                                                                                                                                                                         cgcggcagtAGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACTGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACAGGAGCGAAGCCGAGCGTCAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCATCTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCt                               <  2:13500/250‑1 (MQ=255)
                                                                                                                                                                                                         agcggcagcTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAg                                            >  1:367909/1‑238 (MQ=255)
                                                                                                                                                                                                         agcggcagcTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAg                                            <  2:367909/238‑1 (MQ=255)
                                                                                                                                                                                                         agcggcagcTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCt                               >  1:392480/1‑251 (MQ=255)
                                                                                                                                                                                                               agcTGGAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTaa                         <  1:253780/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggAGCGGTAATCTTAAAGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAAtt                       <  1:158529/249‑1 (MQ=255)
                                                                                                                                                                                                                                tAAAGGGAATACTGCCACGACTCGGGAGTTTTCGACCGGCCTTAACTCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCat           <  1:224001/248‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaGGGAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGc       <  1:346495/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gAATACTGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTg  <  1:307158/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
AACGGCCTGGTTAACGCAGCAGATAACTATCAGAACGCGGGCTGGCGCACCTTTTTTGTTGATAATCAGGTGAAGACGCTGATTAGCGAGGCGCTGGTGAATAATCGGGATTTGCGCATGGCGACGCTGAAAGTGCAGGAGGCACGGGCGCAATATCGTCTGACCGATGCCGACCGCTACCCACAGCTCAATGGCGAGGGCAGCGGCAGCTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTG  >  NC_000913/595717‑596197

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