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.

TCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  NC_000913/595421‑595615
                                                                                                    |                                                                                               
tCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAtt                                                                                                 >  1:2043923/1‑100 (MQ=255)
 cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:626385/1‑100 (MQ=255)
 cTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  1:1240755/1‑100 (MQ=255)
 cTTCGACAATAAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                >  2:269281/1‑100 (MQ=255)
    cggcAATCAACCGGTTCAGATTTCCTTCGCAGGTTGGCCGGGGCATTTTTGGGGAGGGGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAttt                                                                                                <  1:2462466/95‑1 (MQ=255)
    cGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGt                                                                                             <  2:2506341/100‑1 (MQ=255)
      aCAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCa                                                                                           <  2:3391321/100‑1 (MQ=255)
      aCAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCa                                                                                           <  1:379828/100‑1 (MQ=255)
        aaTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAtt                                                                                         >  2:2127949/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGTCAAAATGACAATTTTGTCAtt                                                                                         >  1:942304/1‑99 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:2758308/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCAttt                                                                                        >  1:1564236/1‑100 (MQ=255)
         aTCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttt                                                                                        >  2:3741710/1‑100 (MQ=255)
          tCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAtttt                                                                                       <  1:2757229/100‑1 (MQ=255)
            aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:3786380/1‑100 (MQ=255)
            aaCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTc                                                                                     >  2:729237/1‑100 (MQ=255)
             aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:1359486/1‑100 (MQ=255)
             aCAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCt                                                                                    >  2:3212446/1‑100 (MQ=255)
              cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:3152600/100‑1 (MQ=255)
              cAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:58966/100‑1 (MQ=255)
              cAGTTTCAAATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGTATCAGATGGCAATCGCTTATTGGCAAAATGTCAAAATTGTCAttttgcca                                                                                   <  1:3544261/100‑5 (MQ=255)
               aGTTTCATATTTCCTCCGCATGTTGCCCGTGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTg                                                                                   <  1:1689267/99‑1 (MQ=255)
               aGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGt                                                                                  <  2:487953/100‑1 (MQ=255)
                gTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTc                                                                                 >  2:1989221/1‑100 (MQ=255)
                 tttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCa                                                                                >  2:1329087/1‑100 (MQ=255)
                  ttCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCAc                                                                               <  1:1454436/100‑1 (MQ=255)
                   tCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAAAATTGTCAttttgccaataag                                                                              <  1:2073831/100‑10 (MQ=38)
                     atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  2:60441/1‑100 (MQ=255)
                     atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCgg                                                                            >  1:888651/1‑100 (MQ=255)
                     atatTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTAATTTTTCTGTCACCg                                                                             >  1:3267961/1‑99 (MQ=255)
                        tttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaa                                                                         >  1:362437/1‑100 (MQ=255)
                        attCATCTGCACGTTGCCCGGGCAATTCAAAAATAGCGGGATCAGATGGCAGTCAATTATTAGAAAAATGACAATGTTGTCGTTTTACTgaaac                                                                               >  1:2309718/2‑90 (MQ=255)
                         ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  1:3085135/100‑1 (MQ=255)
                         ttCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGaaaa                                                                        <  2:2997193/100‑1 (MQ=255)
                            ctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCa                                                                     <  2:2768733/100‑1 (MQ=255)
                            ctccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAAGCa                                                                     <  1:664789/100‑1 (MQ=255)
                             tccGCATGTTGTCCGGGCAATTCTAGAGTAGTTTTATCAGATGGCAATCGCTTATTGGAAAAATGACAATTTTGTCATTTTTCTTTCACCGGAAAATCa                                                                     <  1:1229646/99‑1 (MQ=255)
                              ccGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCaga                                                                   >  2:789645/1‑100 (MQ=255)
                               cGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCagag                                                                  <  2:75110/100‑1 (MQ=255)
                                gCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGc                                                                 >  2:3258531/1‑100 (MQ=255)
                                     ttGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGc                                                            <  1:2406796/100‑1 (MQ=255)
                                      tGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCg                                                           >  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)
                                                                                                    |                                                                                               
TCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATGACAATTTTGTCATTTTTCTGTCACCGGAAAATCAGAGCCTGGCGAGTAAAGTTGGCGGCATAAAATCACCAGAAATTATGAGCCTATGTCTCCTTGTAAAC  >  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: