Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 595,520 T→A 19.0% 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,5200TA19.0% 344.9 / 67.2 150intergenic (‑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 (54/66);  new base A (19/10);  total (74/76)
Fisher's exact test for biased strand distribution p-value = 6.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 158 total aligned reads displayed.

TCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595421‑595615
                                                                                                    |                                                                                               
tCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAtt                                                                                                 >  1:2043923/1‑100 (MQ=255)
 cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:626385/1‑100 (MQ=255)
 cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:1240755/1‑100 (MQ=255)
 cTTCGACAATAAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  2:269281/1‑100 (MQ=255)
    cggcAATCAACCGGTTCAGATTTCCTTCGCAGGTTGGCCGGGGCATTTTTGGGGAGGGGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                <  1:2462466/95‑1 (MQ=255)
    cGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGt                                                                                             <  2:2506341/100‑1 (MQ=255)
      aCAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCa                                                                                           <  2:3391321/100‑1 (MQ=255)
      aCAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCa                                                                                           <  1:379828/100‑1 (MQ=255)
        aaTCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtt                                                                                         >  2:2127949/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGTCAAAATGACAATTTTGTCAtt                                                                                         >  1:942304/1‑99 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:2758308/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:1564236/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttt                                                                                        >  2:3741710/1‑100 (MQ=255)
          tCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAtttt                                                                                       <  1:2757229/100‑1 (MQ=255)
            aaCAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:3786380/1‑100 (MQ=255)
            aaCAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:729237/1‑100 (MQ=255)
             aCAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:1359486/1‑100 (MQ=255)
             aCAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:3212446/1‑100 (MQ=255)
              cAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:3152600/100‑1 (MQ=255)
              cAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:58966/100‑1 (MQ=255)
              cAGTTTCAAATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGTATCAGATGGCAATCGCTTATTGGCAAAATGTCAAAATTGTCAttttgcca                                                                                   <  1:3544261/100‑5 (MQ=255)
               aGTTTCATATTTCCTCCGCATGTT‑GCCCGTGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:1689267/99‑1 (MQ=255)
               aGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGt                                                                                  <  2:487953/100‑1 (MQ=255)
                gTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTc                                                                                 >  2:1989221/1‑100 (MQ=255)
                 tttCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCa                                                                                >  2:1329087/1‑100 (MQ=255)
                  ttCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                               <  1:1454436/100‑1 (MQ=255)
                   tCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgccaataag                                                                              <  1:2073831/100‑10 (MQ=38)
                     atatTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  2:60441/1‑100 (MQ=255)
                     atatTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  1:888651/1‑100 (MQ=255)
                     atatTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTAATTTTTCTGTCACCg                                                                             >  1:3267961/1‑99 (MQ=255)
                        tttCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaa                                                                         >  1:362437/1‑100 (MQ=255)
                        attCATCTGCACGTT‑GCCCGGGCAATTCAAAAATAGCGGGATCAGATGGCAGTCAATTATTAGAAAAATGACAATGTTGTCGTTTTACTgaaac                                                                               >  1:2309718/2‑90 (MQ=255)
                         ttCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  1:3085135/100‑1 (MQ=255)
                         ttCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  2:2997193/100‑1 (MQ=255)
                            ctccGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                     <  2:2768733/100‑1 (MQ=255)
                            ctccGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAGCa                                                                     <  1:664789/100‑1 (MQ=255)
                             tccGCATGTT‑GTCCGGGCAATTCTAGAGTAGTTTTATCAGATGGCAATCGCTTATTGGAAAAATGACAATTTTGTCATTTTTCTTTCACCGGAAAATCa                                                                     <  1:1229646/99‑1 (MQ=255)
                              ccGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCaga                                                                   >  2:789645/1‑100 (MQ=255)
                               cGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCagag                                                                  <  2:75110/100‑1 (MQ=255)
                                gCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGc                                                                 >  2:3258531/1‑100 (MQ=255)
                                     tt‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                            <  1:2406796/100‑1 (MQ=255)
                                      t‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCg                                                           >  2:101764/1‑100 (MQ=255)
                                          ccGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGt                                                        <  2:3480724/100‑1 (MQ=255)
                                          ccGGGCAATTCTAAAGTAGCGGGATCAGATGGCCATCACTTATTGGAAAAATGAAAATATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAg                                                         >  1:3860316/1‑99 (MQ=255)
                                           cGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTa                                                       >  2:3439059/1‑100 (MQ=255)
                                            gggCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGCAAAATCAGAGCCTGGCGAGTaa                                                      <  1:1908708/100‑1 (MQ=255)
                                              gCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAg                                                    <  1:3786348/100‑1 (MQ=255)
                                                  ttCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGt                                                   <  1:1418898/97‑1 (MQ=255)
                                                   tCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggc                                               <  1:1007915/100‑1 (MQ=255)
                                                    cTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcg                                              >  2:3879471/1‑100 (MQ=255)
                                                      agagTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTggcggc                                            <  2:1918874/100‑1 (MQ=255)
                                                          tAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaa                                        <  2:3300784/100‑1 (MQ=255)
                                                           aGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaa                                       >  1:867267/1‑100 (MQ=255)
                                                            gCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATaaaa                                      <  1:3342737/100‑1 (MQ=255)
                                                               ggATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCa                                   <  2:3669421/100‑1 (MQ=255)
                                                                 aTCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                 <  2:1411939/100‑1 (MQ=255)
                                                                  tCAGGTGGCAATCGCTTATTTGCAAAATGACAAGTTTGTCATTTTTCTGTCACCGGAAAATCAGAGTCTGGCGAGTAAAGTTGGCGGCATAAAATCAcc                                 <  1:920052/99‑1 (MQ=255)
                                                                  tCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGGCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCa                                <  1:1583141/100‑1 (MQ=255)
                                                                   cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                               >  1:1163469/1‑100 (MQ=255)
                                                                   cAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAg                               >  1:2738790/1‑100 (MQ=255)
                                                                    agctGGCACTCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGCCGAGTAAAGTTTACGGCCTAAAATCACCa                                >  1:752612/4‑98 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGCGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2607587/1‑100 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:1042382/1‑100 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2011732/1‑100 (MQ=255)
                                                                     gATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGaa                             >  1:2876750/1‑100 (MQ=255)
                                                                       tGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAAt                           <  2:3056676/100‑1 (MQ=255)
                                                                         gCAATCGCTTATTGGCAAAATGACAATTTTGTGATTTTTATGGCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTa                         <  1:2724724/100‑1 (MQ=255)
                                                                         gCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAACGTTGGCGGCATAAAATCACCAGAAATTa                         >  1:456221/1‑100 (MQ=255)
                                                                         gCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTa                         <  2:2841220/100‑1 (MQ=255)
                                                                          cAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTAt                        <  2:1908607/100‑1 (MQ=255)
                                                                             tCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATg                       <  1:1158156/98‑1 (MQ=255)
                                                                                   nacagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg               >  1:3702549/7‑100 (MQ=38)
                                                                                   aTTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATg               <  2:948134/100‑1 (MQ=255)
                                                                                    ttGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGt              >  1:678727/1‑100 (MQ=255)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:2666789/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:3366515/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:439755/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:3659892/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  2:263094/5‑100 (MQ=39)
                                                                                     cagaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTCGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:859085/5‑100 (MQ=39)
                                                                                      ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtct            >  1:2897760/1‑100 (MQ=255)
                                                                                      ggCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtc             >  1:2292796/1‑99 (MQ=255)
                                                                                       gtaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:2624885/3‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAACTCAGAGCCTGGCGAGTAAAGTTAGCGGCATAAAAACACCAGAAATAATAAGCCTATGtct            >  1:446620/3‑99 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:1706348/3‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:2673006/3‑100 (MQ=255)
                                                                                       gaaaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGAATAAAATCCCCAGAAATTATGAGCCTATGtct            >  1:1308789/3‑99 (MQ=255)
                                                                                       gCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGtctc           >  2:3442975/1‑100 (MQ=255)
                                                                                        cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          >  1:742739/1‑100 (MQ=255)
                                                                                        cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          >  2:5341/1‑100 (MQ=255)
                                                                                        cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          <  1:1784161/100‑1 (MQ=255)
                                                                                        cAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTcc          <  2:2340729/100‑1 (MQ=255)
                                                                                         aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         <  2:1871145/100‑1 (MQ=255)
                                                                                         aaaaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         <  2:888659/100‑1 (MQ=255)
                                                                                         aaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         >  2:2364344/1‑100 (MQ=255)
                                                                                         aaaaTGACAAAATTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCt         >  2:2247499/1‑100 (MQ=255)
                                                                                           aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  2:3781303/100‑1 (MQ=255)
                                                                                           aaTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  1:1658827/100‑1 (MQ=255)
                                                                                            aTGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTg       <  2:466128/99‑1 (MQ=255)
                                                                                              gACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTaa    >  1:1821868/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTT‑GCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595421‑595615

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