Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,160 G→T 17.6% D87Y (GAC→TAC) ‡ ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1600GT17.6% 425.1 / 66.8 168D87Y (GAC→TAC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (69/68);  new base T (16/14);  total (86/82)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

AAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569061‑569259
                                                                                                   |                                                                                                   
aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                     <  2:4166910/100‑1 (MQ=255)
aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                     <  1:3062757/100‑1 (MQ=255)
  aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCACGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGAc                                                                                                   <  2:250295/100‑1 (MQ=255)
     aaaGCATGCACAATCAGGAGCTTTTTCGTAATTTTTAAATGCGATAGCGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                 <  1:1767749/99‑1 (MQ=255)
     aaaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                >  2:3547748/1‑100 (MQ=255)
       aGCATGCACAATCAGGAGCTTTTTCGGAATGTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTACAGa                                                                                                >  1:4261792/1‑98 (MQ=255)
         cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                            >  1:3546412/1‑100 (MQ=255)
          aTGCACAATCAGGAGCTTTTTCGGATTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                           <  1:1821643/100‑1 (MQ=255)
           tGCACAATCAGGAGCTTTTTTGTTTTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          >  1:1501265/1‑100 (MQ=255)
           tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          <  2:1935512/100‑1 (MQ=255)
           tGCACAATCAGGAGCTTTTTCGGAATGTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          <  2:2557283/100‑1 (MQ=255)
           tGCACAATCAAGAGCTTTTTCTGTATTTTTAGATGATATAGATCATGGTTATATATTGACATTAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          <  1:1545505/100‑1 (MQ=255)
            gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                         <  2:2200797/100‑1 (MQ=255)
            gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                         <  1:2060615/100‑1 (MQ=255)
            gCACAATCAGGAGCTTTTTCGGAATTTGTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                         <  1:3834836/100‑1 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:2282896/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:904446/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:599188/1‑100 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        <  1:3714495/100‑1 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:719140/1‑100 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  1:1115300/1‑100 (MQ=255)
              acaATGAGGAGCTTTTTCGGATTTTTTAGAGGCTATTGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        <  1:2874485/99‑1 (MQ=255)
              acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCAaga                                                                                       >  2:1354126/1‑100 (MQ=255)
               caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:1275816/98‑2 (MQ=255)
               caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                      >  2:92537/1‑100 (MQ=255)
                 aTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaa                                                                                    <  1:762838/100‑1 (MQ=255)
                  tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTATATTGTTAGTTGAAAGTCTGTCCAGTCTTTCaactaaca                                                                                   <  1:2633605/100‑7 (MQ=38)
                  tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaa                                                                                   >  1:1435701/1‑100 (MQ=255)
                    aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                 <  2:289052/100‑1 (MQ=255)
                    aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                 <  2:140304/100‑1 (MQ=255)
                         cTTTTTCGGAATTTTTAGATGCTATAGAGCATTGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                             <  1:2822631/99‑1 (MQ=255)
                         cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                            <  2:3222155/100‑1 (MQ=255)
                         cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                            >  2:2171935/1‑100 (MQ=255)
                          tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                           >  2:143313/1‑100 (MQ=255)
                              tCGGAATTTTTAGATGCTATAGAGCATGGTTATATCTTGCCAGGAACTAAATTGTTAGTTGAAAGTCTGTAAAAACTTTCAAAAGAAAAAGtggaagaag                                                                       >  1:4258189/1‑92 (MQ=255)
                                ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCg                                                                     >  1:3336959/1‑100 (MQ=255)
                                 gAATTTTTCGATGCTAAAGAGAATGGTTATATATTGCCAGGGCCTCGCTTGTTAGTTGAACGTCTGGACGGAATTTCGTGAGAAGTCGTCGGTgccgcaa                                                                    <  2:1911396/100‑7 (MQ=255)
                                  aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                   <  2:2478237/100‑1 (MQ=255)
                                   aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                  <  2:696892/100‑1 (MQ=255)
                                   aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                  <  1:1987753/100‑1 (MQ=255)
                                     ttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                >  2:3672007/1‑100 (MQ=255)
                                      tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                               <  2:496774/100‑1 (MQ=255)
                                      tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                               >  1:2723697/1‑100 (MQ=255)
                                       aaccaTGCTCTATAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAAc                                                              >  2:1625899/5‑100 (MQ=255)
                                         aGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGt                                                            >  1:872344/1‑100 (MQ=255)
                                          gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                           >  1:1257319/1‑100 (MQ=255)
                                          gATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                           <  2:3392328/100‑1 (MQ=255)
                                               tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaaa                                                      >  2:4188097/1‑100 (MQ=255)
                                               tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaaa                                                      >  1:2835897/1‑100 (MQ=255)
                                                    agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  2:2976198/100‑1 (MQ=255)
                                                     gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                <  1:1354122/100‑1 (MQ=255)
                                                     gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                <  2:151250/100‑1 (MQ=255)
                                                       aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                              >  1:3408199/1‑100 (MQ=255)
                                                        tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTg                                             >  2:3072347/1‑100 (MQ=255)
                                                         ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                            >  1:693990/1‑100 (MQ=255)
                                                          gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGaa                                           <  1:3229879/100‑1 (MQ=255)
                                                             atatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                        >  1:3028917/1‑100 (MQ=255)
                                                              tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAGGAGAAAAAGTCGGTGCAGCGATTGACCGTCTGAAATTGATTTTGAATCAc                                       >  2:2072161/1‑100 (MQ=255)
                                                              tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                       <  2:3087680/100‑1 (MQ=255)
                                                              tatatTGCCAGGAACTACATGGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTAAACGTCTGAAATTGATTTTGAATCAc                                       >  1:4005388/1‑100 (MQ=255)
                                                                tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                     <  1:2377498/100‑1 (MQ=255)
                                                                tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                     >  2:2923523/1‑100 (MQ=255)
                                                                tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                     >  2:560436/1‑100 (MQ=255)
                                                                tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                     <  2:3441584/100‑1 (MQ=255)
                                                                       aGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAt                              <  2:3519107/100‑1 (MQ=255)
                                                                        ggAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                             <  2:4170781/100‑1 (MQ=255)
                                                                         gAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGt                            <  1:1094316/100‑1 (MQ=255)
                                                                           aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTta                          <  2:3798222/100‑1 (MQ=255)
                                                                            cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                         >  2:1778299/1‑100 (MQ=255)
                                                                            cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                         >  1:2404949/1‑100 (MQ=255)
                                                                            cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                         >  2:2686789/1‑100 (MQ=255)
                                                                             tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                        <  2:652110/100‑1 (MQ=255)
                                                                              tcATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAATAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTCGAATCACGGTATTGATGTTATaa                       >  1:122772/2‑100 (MQ=255)
                                                                               cATTGTTAGTTGAAAGTCTGGTCAGCCTTTCAAGAGAAAAAGTCGGTGAAGCAATTGAACGTCTGAAATTGATATTGAATAAAGGTATTGATGTta                          >  1:3384101/1‑96 (MQ=255)
                                                                                    ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                 <  1:2171930/100‑1 (MQ=255)
                                                                                     tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                <  1:1033459/100‑1 (MQ=255)
                                                                                     tAGTTGAAAGTCTGGACAGAAATAAAAGAAAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                <  1:2339739/100‑1 (MQ=255)
                                                                                      tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTTTAACTCTTTGc               >  1:2918182/3‑100 (MQ=39)
                                                                                      tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               >  2:1102105/3‑100 (MQ=255)
                                                                                      tcttGAAAGTCTGTCAAGACTTTAAAGAGAAAAAGTCGGTGAAGCTATTGAACGTCTGAAAATGATTTTGAATCACGTTATTGATGTTATAACACTTTGc               >  1:3912182/3‑100 (MQ=255)
                                                                                      aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               <  1:1579465/100‑1 (MQ=255)
                                                                                      aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               <  2:1224322/100‑1 (MQ=255)
                                                                                      aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               <  2:763886/100‑1 (MQ=255)
                                                                                       gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:1272054/1‑100 (MQ=255)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  1:2202631/2‑100 (MQ=255)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  1:2454496/2‑100 (MQ=255)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:2022258/2‑100 (MQ=255)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:3862382/2‑100 (MQ=255)
                                                                                        ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa             >  1:3590830/1‑100 (MQ=255)
                                                                                         tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGGCTGAAATTGCTTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  2:410887/1‑100 (MQ=255)
                                                                                          gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAGTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa           <  1:2987948/100‑1 (MQ=255)
                                                                                               tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa      >  1:556037/1‑100 (MQ=255)
                                                                                                cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg     <  2:3916780/100‑1 (MQ=255)
                                                                                                cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg     <  2:3471823/100‑1 (MQ=255)
                                                                                                 tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    >  2:1415069/1‑100 (MQ=255)
                                                                                                 tGGACAGACTTGCAAGAGAAAGAGTCGGTGAAGCGATTGGACGTCTTAAATTGATTTTGAATCACGGTATTGATGTTATAACTATTTGCGACAATACAGt    <  1:3466612/100‑1 (MQ=255)
                                                                                                  ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   >  1:3216634/1‑100 (MQ=255)
                                                                                                  ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   >  1:857005/1‑100 (MQ=255)
                                                                                                   gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt  <  2:1250956/100‑1 (MQ=255)
                                                                                                   gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt  >  2:2049539/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569061‑569259

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: