Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,160 G→T 18.0% 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,1600GT18.0% 329.2 / 56.7 129D87Y (GAC→TAC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (52/53);  new base T (10/13);  total (63/66)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 137 total aligned reads displayed.

AAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569061‑569259
                                                                                                   |                                                                                                   
aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                     <  1:541074/100‑1 (MQ=255)
aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                     <  1:736427/100‑1 (MQ=255)
aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTgg                                                                                                     <  2:236294/100‑1 (MQ=255)
  aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGAc                                                                                                   >  1:2174916/1‑100 (MQ=255)
   ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                  <  1:2559123/100‑1 (MQ=255)
    gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTGAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                 <  2:2041174/100‑1 (MQ=255)
       tgCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACt                                                                                              <  1:2375339/99‑1 (MQ=255)
        gCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACta                                                                                             >  2:2375350/1‑99 (MQ=255)
         cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCATGCATTTCTTTGTTAGTTGAAAGTCTGGACAAACttt                                                                                            >  1:2437519/1‑100 (MQ=255)
          aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                           <  2:1176919/100‑1 (MQ=255)
          aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                           <  1:416120/100‑1 (MQ=255)
            gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                         <  2:2719004/100‑1 (MQ=255)
            gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          <  2:2800430/99‑1 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCatc                                                                                        <  1:506881/100‑3 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:1461490/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:786211/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:2360736/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:2356687/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:2690754/100‑2 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGaagatcgga                                                                                        >  2:48637/1‑92 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:181073/1‑100 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        <  1:2497139/100‑1 (MQ=255)
             cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGAAAGGACCTACATTGTAGGATGAAAGTCTGAAAAGACTTTcc                                                                                          >  2:535121/1‑97 (MQ=255)
             cacaATCAGGAGATTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:1913360/1‑100 (MQ=255)
              acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaact                                                                                       <  1:1161363/100‑3 (MQ=255)
              acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAaga                                                                                       <  2:1830211/100‑1 (MQ=255)
              acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAaga                                                                                       >  1:1924924/1‑100 (MQ=255)
               caATCAGGAGCTTTTTCGGAATTTTTCGATGCTATAGAGGATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGCCAGACTTTCaac                                                                                        <  1:761546/98‑2 (MQ=255)
                  tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaa                                                                                   <  1:2867233/100‑1 (MQ=255)
                    aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                 <  2:789418/100‑1 (MQ=255)
                    aGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                 >  2:1659925/1‑100 (MQ=255)
                     ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                <  2:2476281/100‑1 (MQ=255)
                     ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  2:641860/1‑100 (MQ=255)
                     ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  2:2498317/1‑100 (MQ=255)
                          tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                           <  1:2829103/100‑1 (MQ=255)
                          tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                           <  2:276420/100‑1 (MQ=255)
                           ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGta                                                                          >  2:2829113/1‑99 (MQ=255)
                              tCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAg                                                                       >  2:777086/1‑100 (MQ=255)
                                ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCg                                                                     >  1:243143/1‑100 (MQ=255)
                                ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTAAATTGTTAGTTGAAAGTCTGGAAAGACTTTCAAGAGAAAAAGTCGGTGAATCg                                                                     >  1:277848/1‑100 (MQ=255)
                                  aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                   <  2:1267853/100‑1 (MQ=255)
                                   aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                  >  2:1132961/1‑100 (MQ=255)
                                    tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                 <  1:2096358/100‑1 (MQ=255)
                                     ttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                <  1:2539894/100‑1 (MQ=255)
                                        tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                             <  2:860843/100‑1 (MQ=255)
                                         aGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGt                                                            >  2:369294/1‑100 (MQ=255)
                                              cTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                       >  1:1063360/1‑100 (MQ=255)
                                              cTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                       <  1:1890547/100‑1 (MQ=255)
                                                ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTGGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:1323362/1‑100 (MQ=255)
                                                ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:1893774/1‑100 (MQ=255)
                                                ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  1:424904/1‑100 (MQ=255)
                                                    agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  1:1132959/100‑1 (MQ=255)
                                                    agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  1:1168165/100‑1 (MQ=255)
                                                    agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 >  2:2181747/1‑100 (MQ=255)
                                                        tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTg                                             <  2:287725/100‑1 (MQ=255)
                                                        tGGTTATATATTGCCAGGAACTACATTGTGAGTTGGAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                              <  1:1634542/99‑1 (MQ=255)
                                                         ggTTATATATTGCCAGGAATTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                            >  2:2354827/1‑100 (MQ=255)
                                                         ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                            <  1:1024281/100‑1 (MQ=255)
                                                          gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGaa                                           >  2:2213483/1‑100 (MQ=255)
                                                           ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                          <  2:2260841/100‑1 (MQ=255)
                                                            tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATc                                         <  2:1323999/100‑1 (MQ=255)
                                                              tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                       >  2:1358232/1‑100 (MQ=255)
                                                               atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                      >  2:250693/1‑100 (MQ=255)
                                                               atatTGCCAGGAACTACATTGTTAGTTGAAAGTCCGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAGTTTGAATCACg                                      <  1:1525989/100‑1 (MQ=255)
                                                                 atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  1:1201884/1‑100 (MQ=255)
                                                                 atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  1:755385/1‑100 (MQ=255)
                                                                 atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  2:1843724/1‑100 (MQ=255)
                                                                 atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGt                                    >  2:1216676/1‑100 (MQ=255)
                                                                  ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGAATATAAAAAAAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                   <  1:899769/100‑1 (MQ=255)
                                                                    gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                 >  1:386917/1‑100 (MQ=255)
                                                                     ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                >  1:2784129/1‑100 (MQ=255)
                                                                     ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                <  1:181071/100‑1 (MQ=255)
                                                                               cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                      <  1:369292/100‑1 (MQ=255)
                                                                                 ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActc                    <  2:2006780/100‑1 (MQ=255)
                                                                                   ttctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                  >  1:335770/6‑100 (MQ=38)
                                                                                   gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                  >  2:740700/1‑100 (MQ=255)
                                                                                     tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:392601/1‑100 (MQ=255)
                                                                                     ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:2227962/4‑100 (MQ=39)
                                                                                     ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                >  2:631303/4‑100 (MQ=39)
                                                                                      tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               >  1:1993138/3‑100 (MQ=255)
                                                                                      aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               <  1:460268/100‑1 (MQ=255)
                                                                                       gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATTTTATAACTCTTTGCg              >  2:2861763/1‑100 (MQ=255)
                                                                                       gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  1:2866829/1‑100 (MQ=255)
                                                                                       gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATATTATAACTCTTTGCg              >  2:2864271/1‑100 (MQ=255)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGGGATTGAACGTCTGAAGTTGATTTTGAATCACGGTATTGATGTTCTAACTCttt                 >  1:2370016/2‑97 (MQ=38)
                                                                                       cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:1401692/2‑100 (MQ=255)
                                                                                         tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa             >  1:2621736/1‑99 (MQ=255)
                                                                                         tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  2:883407/1‑100 (MQ=255)
                                                                                         tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            <  1:250691/100‑1 (MQ=255)
                                                                                         tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  1:1929819/1‑100 (MQ=255)
                                                                                         tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAGGTTATAACTCTTTGCGAc            >  1:2784759/1‑100 (MQ=255)
                                                                                          gAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa           >  2:154687/1‑100 (MQ=255)
                                                                                           aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAACGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa          >  1:1788394/1‑100 (MQ=255)
                                                                                            aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt         <  1:1358230/100‑1 (MQ=255)
                                                                                              gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc       <  2:1063362/100‑1 (MQ=255)
                                                                                              gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc       <  1:2852747/100‑1 (MQ=255)
                                                                                              gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc       <  1:2852451/100‑1 (MQ=255)
                                                                                                cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg     <  2:386919/100‑1 (MQ=255)
                                                                                                 tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    >  2:2453993/1‑100 (MQ=255)
                                                                                                 nggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    >  2:1585580/2‑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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: