Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,582,272 G→A 20.8% 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,2720GA20.8% 385.7 / 87.7 187C137Y (TGT→TAT) yhhZputative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base G (80/68);  new base A (20/19);  total (100/87)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 189 total aligned reads displayed.

TACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATC  >  NC_000913/3582173‑3582370
                                                                                                   |                                                                                                  
taCAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACTAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTg                                                                                                    >  2:899221/1‑100 (MQ=255)
  cAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTc                                                                                                  <  1:327213/100‑1 (MQ=255)
  cAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAGTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTc                                                                                                  <  2:2668991/100‑1 (MQ=255)
    aCTAAGAGGCGCTTTTTTTTTTTCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCaa                                                                                                >  2:3603814/1‑100 (MQ=255)
    aCTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCaa                                                                                                >  2:1844515/1‑100 (MQ=255)
    aCTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCa                                                                                                 <  2:1341156/99‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACa                                                                                              <  2:2965867/100‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACa                                                                                              <  1:2896098/100‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACa                                                                                              <  1:3374904/100‑1 (MQ=255)
      tAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACa                                                                                              <  2:3224256/100‑1 (MQ=255)
       aaGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACag                                                                                             >  1:4096742/1‑99 (MQ=255)
       aaGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACAt                                                                                             <  1:866660/100‑1 (MQ=255)
       aaGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACAt                                                                                             >  1:806457/1‑100 (MQ=255)
       aaGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACAt                                                                                             <  1:3772067/100‑1 (MQ=255)
            gcgcTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTAt                                                                                        >  1:3915419/1‑100 (MQ=255)
                ttttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCa                                                                                    <  1:3716537/100‑1 (MQ=255)
                ttttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCa                                                                                    <  1:2259581/100‑1 (MQ=255)
                ttttttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCa                                                                                    <  1:31339/100‑1 (MQ=255)
                   tttATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAAt                                                                                 >  1:717494/1‑100 (MQ=255)
                     tATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATAc                                                                               >  2:100327/1‑100 (MQ=255)
                      aTCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATAcc                                                                              >  1:205600/1‑100 (MQ=255)
                      aTCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATAcc                                                                              >  1:2262514/1‑100 (MQ=255)
                       tCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCg                                                                             <  2:249262/100‑1 (MQ=255)
                       tCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCg                                                                             >  2:312847/1‑100 (MQ=255)
                        cgtcTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTAATGAATTTATCCTCTGTCAACATCTTATCGCAAATAccaa                                                                            >  1:3105878/4‑98 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAg                                                                           <  2:175783/100‑1 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAg                                                                           >  2:3115494/1‑100 (MQ=255)
                         ggCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAg                                                                           <  1:3218735/100‑1 (MQ=255)
                            tATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTc                                                                        >  2:2777743/1‑100 (MQ=255)
                             aTTCATCACCAGATCATTGAAAACCAACTGGATACAGACACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCa                                                                       <  2:2262500/100‑1 (MQ=255)
                                catcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCt                                                                    >  1:1033916/1‑100 (MQ=255)
                                catcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCt                                                                    >  1:1133101/1‑100 (MQ=255)
                                catcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCt                                                                    >  1:126565/1‑100 (MQ=255)
                                 atcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCTATTTCATCTa                                                                   >  1:776489/1‑100 (MQ=255)
                                  tcaCCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTAt                                                                  <  2:3637720/100‑1 (MQ=255)
                                      cAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTgg                                                              <  2:3426200/100‑1 (MQ=255)
                                       aGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGc                                                             <  1:1879695/100‑1 (MQ=255)
                                       aGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGc                                                             >  1:987304/1‑100 (MQ=255)
                                           cATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTc                                                         >  2:3394859/1‑100 (MQ=255)
                                             ttGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTccc                                                       >  2:1312247/1‑100 (MQ=255)
                                             ttGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTccc                                                       <  2:1622672/100‑1 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTg                                                     >  2:2439953/1‑100 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTg                                                     >  2:4065550/1‑100 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTg                                                     >  1:1496320/1‑100 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAATAACTATTAGTAATGAATTTATCCTCTGTCAAGATATTATCGCAAATACCAAGTTCAGCTATTTAGCACTCACTg                                                     >  2:3107759/1‑100 (MQ=255)
                                               gAAAACCAACTGGATACAGAAACAAGAAATATTAGTTGTGAATTTATCCTCTGTCAACAACTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCt                                                      >  2:3476380/1‑99 (MQ=255)
                                                 aaaCCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGaa                                                   <  2:1531804/100‑1 (MQ=255)
                                                  aaCCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGaaa                                                  >  1:354253/1‑100 (MQ=255)
                                                  aaCCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGaaa                                                  <  1:3599409/100‑1 (MQ=255)
                                                   aCCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGaaaa                                                 <  1:1000303/100‑1 (MQ=255)
                                                    ccAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAAt                                                >  1:464379/1‑100 (MQ=255)
                                                      aaCTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATta                                              <  2:1171117/100‑1 (MQ=255)
                                                      aaCTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATta                                              <  1:1918282/100‑1 (MQ=255)
                                                      aaCTGGATACAGAAACAATAACTATCAGTTATGAATTTATCCTCTGTCAACATCTTATCGCTAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATta                                              <  1:1771632/100‑1 (MQ=255)
                                                       aCTGGATACAGAAACAATAACTATTAGTTATCAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATtat                                             >  2:220255/1‑100 (MQ=255)
                                                        cTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATtata                                            <  2:1910266/100‑1 (MQ=255)
                                                          ggATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAAc                                          >  2:1953526/1‑100 (MQ=255)
                                                            aTACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCg                                        <  1:3041078/100‑1 (MQ=255)
                                                             tACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgt                                       <  2:54702/100‑1 (MQ=255)
                                                             tACAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAAATATAACAgt                                       >  2:1189550/1‑100 (MQ=255)
                                                               cAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgttt                                     >  1:215839/1‑100 (MQ=255)
                                                               cAGAAACAATAACTATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgttt                                     >  1:2498250/1‑100 (MQ=255)
                                                                aGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttg                                    <  1:3719268/100‑1 (MQ=255)
                                                                aGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttg                                    >  2:180245/1‑100 (MQ=255)
                                                                 aaattcaTAACTAATAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgt                                   >  1:2468222/7‑100 (MQ=255)
                                                                  aattcaTAACTAATAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgtt                                  >  2:809426/6‑100 (MQ=255)
                                                                   aaCAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgttt                                 >  2:857088/1‑100 (MQ=255)
                                                                   aaCAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgttt                                 <  2:2902304/100‑1 (MQ=255)
                                                                   aaCAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCgtttgttt                                 <  1:623432/100‑1 (MQ=255)
                                                                     cAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGttttt                               >  2:434895/1‑100 (MQ=255)
                                                                      aaTAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTa                              <  2:2094448/100‑1 (MQ=255)
                                                                      aaTAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTa                              <  1:2196638/100‑1 (MQ=255)
                                                                      aaTAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTa                              <  1:1846640/100‑1 (MQ=255)
                                                                       aTAACTATTAGTTATGAATATATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGTTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTAc                             <  2:3966497/100‑1 (MQ=255)
                                                                        tAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATATCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTAcc                            <  2:2060723/100‑1 (MQ=255)
                                                                        tAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCATCTATTTGGCACTCCCTGAAAATTATAACCGTTTTTTTTTAc                             >  2:3065937/1‑99 (MQ=255)
                                                                         aaCTATTAGTTATGAATTTATCCCCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCa                           <  1:3318803/100‑1 (MQ=255)
                                                                            tATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCaatt                        <  1:2813715/100‑3 (MQ=255)
                                                                            tATTAGTTATGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTaagatcg                        >  2:318936/1‑94 (MQ=255)
                                                                            tATTAGTTATGAATTTATCCACTGTCAACATCTTATCGCAAATACCTAGTTCAGCTATTTAGTACTCCCTGAAAATTATAACCCTTTGTTTTTACCaaa                         >  2:1442085/1‑99 (MQ=255)
                                                                            gatgaGTGATGAAGTTATCCTTTGTTAACATCCTATCGAGAATACCGAGTTCTGCTAGTTGGCAATCCCCGAAAACTATAACCGTTTGTTTTTa                              <  2:227249/90‑1 (MQ=255)
                                                                             aTTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAAtt                       >  1:1007766/1‑100 (MQ=255)
                                                                             aTTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAAtt                       >  1:1431521/1‑100 (MQ=255)
                                                                              ttAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTc                      <  1:1682120/100‑1 (MQ=255)
                                                                                aGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAAtt                       >  2:4019384/1‑97 (MQ=255)
                                                                                aGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCaa                    >  2:1969333/1‑100 (MQ=255)
                                                                                 gTTATGAATTTATCCTCTATCACCACCCTATCGCAAATACAGAACTCAGCTTATTTGCACTACCTGAAAATTATAACCGTTTGTTTTAACCAAATTc                      <  2:109099/97‑1 (MQ=255)
                                                                                    aTGAATTTATCCTCTGTAAACATCGTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCTAATTCaaaaaa                >  2:744463/1‑100 (MQ=255)
                                                                                     tGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAc               >  1:2726304/1‑100 (MQ=255)
                                                                                     tGAATTTATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAc               <  1:2610121/100‑1 (MQ=255)
                                                                                      gAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAcc              >  1:4009858/1‑100 (MQ=255)
                                                                                       aaTTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCATCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAAcc              <  2:2511964/99‑1 (MQ=255)
                                                                                          ttATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaac          <  2:3931835/100‑1 (MQ=255)
                                                                                           tATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaaca         >  2:1340413/1‑100 (MQ=255)
                                                                                           tATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaaca         <  1:2076878/100‑1 (MQ=255)
                                                                                           tATCCTCTATCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaaca         <  1:1628728/100‑1 (MQ=255)
                                                                                             tCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACcaaacaaa       <  1:1088727/100‑1 (MQ=255)
                                                                                              ccTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAaat      >  2:3291434/1‑100 (MQ=255)
                                                                                                  tGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATc  <  1:3451043/100‑1 (MQ=255)
                                                                                                   gTCAATATGTTATCGCAAATGCCGAGTTAAGCTATTTGGCGCTCCCTGAAAATGCTAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATc  <  2:2172142/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
TACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATC  >  NC_000913/3582173‑3582370

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: