Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,582,272 G→A 52.7% C137Y (TGT→TAT)  yhhZ → putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,582,2720GA52.7% 12.5 / 214.6 161C137Y (TGT→TAT) yhhZputative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base G (43/32);  new base A (48/36);  total (92/69)
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

Only 100 of 163 total aligned reads displayed.

TACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCG  >  NC_000913/3582173‑3582371
                                                                                                   |                                                                                                   
taCAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTg                                                                                                     >  2:854974/1‑100 (MQ=255)
  cttcTAATAGGCGCGTTTTGATCGGCTAGTCATCACCAGATGAGTTTAAACCAACTGTATACATTAACAATAACTATTAGGTGTGAATGTATCCTCTGt                                                                                                    >  1:2488168/4‑99 (MQ=255)
  cAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGTAAACCACCTGGATCCATAAACACTAACTCTTAGTTATGACTTTATCCTCTGTc                                                                                                   >  1:1807808/1‑100 (MQ=255)
  cAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGTATACAGAAACTATAACTATTAGTTATGATTTTATCCTCTGTc                                                                                                   <  2:2880579/100‑1 (MQ=255)
  cAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTc                                                                                                   >  1:2062703/1‑100 (MQ=255)
   aaCTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCa                                                                                                  <  2:3920291/100‑1 (MQ=255)
   aaCTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCa                                                                                                  <  1:4092622/100‑1 (MQ=255)
   aaCTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATCAATTTATCCTCTATCa                                                                                                  >  1:3036587/1‑100 (MQ=255)
    aCTAAGAGGCGCTTTTTTATCGGTTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCaa                                                                                                 <  1:161624/100‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACa                                                                                               <  2:92677/100‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACGATTAGTTATGAATTTATCCTCTATCAACa                                                                                               <  2:2755091/100‑1 (MQ=255)
          agGCGCTTTTTTATTTTCTTTTCATCAACAGATCATTGAAAAAAAACTGGATATTGAAACAATAACTATTAGTTATAAATTTATACTATATACACATCtt                                                                                           >  1:406764/1‑100 (MQ=255)
          agGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCtt                                                                                           >  1:1769657/1‑100 (MQ=255)
              gcTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCg                                                                                       >  1:655849/1‑100 (MQ=255)
              gcTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCg                                                                                       >  2:637272/1‑100 (MQ=255)
               cTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGc                                                                                      >  1:3584200/1‑100 (MQ=255)
                ttttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCa                                                                                     >  1:1423694/1‑100 (MQ=255)
                 tttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCa                                                                                     <  1:1884833/99‑1 (MQ=255)
                  ttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCaaa                                                                                   >  2:1797167/1‑100 (MQ=255)
                  ttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCaaa                                                                                   >  2:1557903/1‑100 (MQ=255)
                  ttttATCGGCTATTCATCACCACATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGCTATGAATTTATCCTCTATCAACATCTTATCGCaaa                                                                                   >  1:4019558/1‑100 (MQ=255)
                     tATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATAc                                                                                >  2:1348624/1‑100 (MQ=255)
                       tCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCg                                                                              >  2:3644494/1‑100 (MQ=255)
                        cGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGa                                                                             <  1:4068798/100‑1 (MQ=255)
                        cGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGa                                                                             >  2:4243306/1‑100 (MQ=255)
                        cGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGa                                                                             <  1:3392338/100‑1 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAGTTTATCCTCTGTCAACATCTTATCGCAAATACCGAg                                                                            >  2:2998489/1‑100 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAg                                                                            >  2:532175/1‑100 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAAATATTAGTTGTGAATTTATCCTCTGTCATCATCTTATCGCAAATACCGAg                                                                            >  1:449610/1‑100 (MQ=255)
                          gCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGt                                                                           <  1:376864/100‑1 (MQ=255)
                            tATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTc                                                                         <  2:2686393/100‑1 (MQ=255)
                              ttCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAg                                                                       <  2:93412/100‑1 (MQ=255)
                               tcatcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGc                                                                      >  1:2428956/1‑100 (MQ=255)
                               tcatcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGc                                                                      <  2:1744990/100‑1 (MQ=255)
                               tcatcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTCTCGCAAATACCGAGTTCAGc                                                                      <  1:2606173/100‑1 (MQ=255)
                               tcatcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGc                                                                      >  1:1174920/1‑100 (MQ=255)
                               tcatcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGc                                                                      <  2:2921663/100‑1 (MQ=255)
                                catcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTTTTTTTTTTTTTTATCTTATCGCAAATACCGAGTTCAGCt                                                                     >  1:259324/1‑100 (MQ=255)
                                catcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCt                                                                     >  2:3036191/1‑100 (MQ=255)
                                 atcaCCACATCATTGAAAACCAACTGTATACAGAAACAATAACTATTAGTTATGAAATTATCCGCTATCAACATCTTATCGCAAATACCGACTTCAGc                                                                      <  1:270517/98‑1 (MQ=255)
                                   caCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTAtt                                                                  >  1:479864/1‑100 (MQ=255)
                                    aCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTAttt                                                                 <  2:4113481/100‑1 (MQ=255)
                                       aGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGc                                                              <  2:1941849/100‑1 (MQ=255)
                                        gATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCa                                                             >  1:3671208/1‑100 (MQ=255)
                                          tCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACt                                                           >  2:1105246/1‑100 (MQ=255)
                                           cATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTc                                                          >  2:1957956/1‑100 (MQ=255)
                                             ttGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTccc                                                        >  1:4108916/1‑100 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTCATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTg                                                      <  1:2286006/100‑1 (MQ=255)
                                                    ccAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAAt                                                 <  1:1421271/100‑1 (MQ=255)
                                                    ccAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAAt                                                 >  1:2874396/1‑100 (MQ=255)
                                                     cAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCCACATCTTATGGCAAATAACAAGTTCAGCTATTTGGCACTCCCTGaaa                                                   >  1:3678008/1‑97 (MQ=255)
                                                     cAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAAtt                                                <  2:1423698/100‑1 (MQ=255)
                                                        cTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATtata                                             <  1:3365711/100‑1 (MQ=255)
                                                          ggATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAAc                                           >  1:1971697/1‑100 (MQ=255)
                                                             tACAGAAACAATACCTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCg                                         >  1:2641029/1‑99 (MQ=255)
                                                               cAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgttt                                      <  2:1873945/100‑1 (MQ=255)
                                                               cAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCACCTATTTGGCACTCCCTGAAAATTATAACCgttt                                      >  1:3851529/1‑100 (MQ=255)
                                                                aGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgttt                                      >  1:1077687/1‑99 (MQ=255)
                                                                 gAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGATTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgctcat                                    >  1:1099372/1‑95 (MQ=255)
                                                                 gAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgt                                    >  2:2701100/1‑100 (MQ=255)
                                                                  aaaCAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgtt                                   >  1:878272/1‑100 (MQ=255)
                                                                  aaaCAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgt                                    >  2:3859511/1‑99 (MQ=255)
                                                                   aaCAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgttt                                  <  1:3861902/100‑1 (MQ=255)
                                                                     cAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGttttt                                <  1:2027369/100‑1 (MQ=255)
                                                                      aaTAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTa                               >  1:2261749/1‑100 (MQ=255)
                                                                      aaTAACTATTAGTTATGAATTTATCCGCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTa                               <  1:1511812/100‑1 (MQ=255)
                                                                       aTAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTAc                              <  2:3722175/100‑1 (MQ=255)
                                                                          aCTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCaa                           >  1:215026/1‑100 (MQ=255)
                                                                          aCTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCaa                           >  2:3860261/1‑100 (MQ=255)
                                                                          aCTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGCCACTCCCTGAAAATTATAACCGTTTGTTTTTActa                            >  1:2356234/1‑97 (MQ=255)
                                                                           cTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCaaa                          >  2:3669244/1‑100 (MQ=255)
                                                                           cTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCaaa                          <  2:4159335/100‑1 (MQ=255)
                                                                            tATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAAt                         <  2:4277812/100‑1 (MQ=255)
                                                                                aGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaa                     >  1:3379084/1‑100 (MQ=255)
                                                                                aGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaa                     >  2:121693/1‑100 (MQ=255)
                                                                                aGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaa                     >  2:3278373/1‑100 (MQ=255)
                                                                                 gTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaaa                    <  2:3129096/100‑1 (MQ=255)
                                                                                 gTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaaa                    >  1:1594867/1‑100 (MQ=255)
                                                                                 gTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaaa                    <  2:2802116/100‑1 (MQ=255)
                                                                                   tATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAAGTCaaaaa                  <  2:2238054/100‑1 (MQ=255)
                                                                                   tATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGACACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaaaaa                  >  2:313096/1‑100 (MQ=255)
                                                                                      gAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAcc               <  2:1594871/100‑1 (MQ=255)
                                                                                      gAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAcc               >  2:3714007/1‑100 (MQ=255)
                                                                                       aaTTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACca              >  1:2507968/1‑100 (MQ=255)
                                                                                       aaTTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACca              >  1:3128096/1‑100 (MQ=255)
                                                                                       aaTTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACca              >  2:2229901/1‑100 (MQ=255)
                                                                                        aattagcttctATCAACTTCTTATCGCAAATACGGAGGTCAGCTATTTGCAATTCCCTGAAATTTGTAACCGTTTGTTTTTAAAAAATTAAAAAAAcc               <  1:782644/90‑1 (MQ=255)
                                                                                        aTTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaa             >  1:3039183/1‑100 (MQ=255)
                                                                                        aTTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaa             >  1:579043/1‑100 (MQ=255)
                                                                                         tttATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaa            >  1:1021809/1‑100 (MQ=255)
                                                                                         tttATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaa            >  1:1021096/1‑100 (MQ=255)
                                                                                          ttATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaac           <  2:1294995/100‑1 (MQ=255)
                                                                                           tATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaaca          >  2:2499312/1‑100 (MQ=255)
                                                                                           tATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaaca          <  2:698068/100‑1 (MQ=255)
                                                                                            aTCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaacaa         >  1:3581064/1‑100 (MQ=255)
                                                                                              ccTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAaat       >  1:2980446/1‑100 (MQ=255)
                                                                                               ctctGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAaata      <  1:4216521/100‑1 (MQ=255)
                                                                                               ctctGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAaata      >  2:96218/1‑100 (MQ=255)
                                                                                                  tGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATc   <  2:3765225/100‑1 (MQ=255)
                                                                                                   gTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCg  <  2:1626954/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCG  >  NC_000913/3582173‑3582371

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: