Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,473,882 T→G 11.7% L94R (CTG→CGG)  yfdP → CPS‑53 (KpLE1) prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,473,8820TG11.7% 170.7 / 2.8 79L94R (CTG→CGG) yfdPCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (25/41);  new base G (9/3);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACTTATTA  >  NC_000913/2473737‑2474030
                                                                                                                                                  |                                                                                                                                                    
tcactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGAc                                                                                                                                                   <  1:1507991/149‑1 (MQ=255)
tcactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCCGAc                                                                                                                                                   >  1:780819/1‑149 (MQ=255)
 cactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACt                                                                                                                                                  <  1:24641/149‑1 (MQ=255)
  actcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTa                                                                                                                                                 >  1:1773019/1‑149 (MQ=255)
  actcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTa                                                                                                                                                 >  1:1193545/1‑149 (MQ=255)
  actcaACCTGGAAATCAAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCTCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTa                                                                                                                                                 >  1:1144270/1‑149 (MQ=255)
          ttgAAATCTAACCATTTGATGAGAACAGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAg                                                                                                                                         <  1:497337/147‑1 (MQ=255)
           ggAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAgg                                                                                                                                        >  1:170955/1‑149 (MQ=255)
           ggAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAgg                                                                                                                                        <  1:1584515/149‑1 (MQ=255)
             aaagcaaaCCATTTCATGAGAACTGTCTGAAAATAAATCACAACACTACATATGTTCGACCGACTAACCGTGAGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTGTGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGaa                                                                                                                                      <  2:411979/143‑1 (MQ=255)
                gcGAACCCATTAATGAAAACCGTCTTAAAACCAAACACAAACTCTCAAATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGAc                                                                                                                                   <  2:255578/148‑1 (MQ=255)
                  gAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCCGGAAGACCa                                                                                                                                 >  1:119148/1‑149 (MQ=255)
                           gatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGCTGTATGTCAATCGCGGTGGTCGGCTGACTATTCTGCCGGAAGACCAGGGACAAtt                                                                                                                        >  1:214267/1‑149 (MQ=255)
                           gatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAAtt                                                                                                                        >  1:1564062/1‑149 (MQ=255)
                           gatgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGAAGGAAGACCAGGGCCAAtt                                                                                                                        >  2:398775/1‑149 (MQ=255)
                            atgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTa                                                                                                                       >  1:577499/1‑149 (MQ=255)
                            atgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTa                                                                                                                       >  1:1107446/1‑149 (MQ=255)
                             ggagaACCGTCTGAAAATCAAACACAAACTCGCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGTGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGGCAATTAc                                                                                                                      >  2:1325464/2‑149 (MQ=255)
                                 aaCCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCCACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGGGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGAc                                                                                                                  >  1:709848/1‑149 (MQ=255)
                                    cGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGActct                                                                                                               <  1:788560/149‑1 (MQ=255)
                                             ttCAAACACAAACACTAATATGTTAGCACGACTAACCCTGGTAAATTTTCCGAAGAAGACACCACCGAAACGCCTATGTATGTCAAGCGCGGTGGACGCCTGACTACTCTGCAGGAAGACCAGGGACAATTACTGCCTCTTGCCGGTGa                                                                                                      <  2:1729189/148‑1 (MQ=255)
                                             aTCAAACACAAACTTCCATATGTGCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGa                                                                                                      <  1:985035/149‑1 (MQ=255)
                                              tCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTACAGGAAAACCCGGGACAGTGACTGCCCCTTGCAGGTGAg                                                                                                     >  2:819970/1‑149 (MQ=255)
                                               ccaaCACAAACTCTCATATGTTCGCCCGACTAACCGTCGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGc                                                                                                    <  2:1439370/147‑1 (MQ=255)
                                               cAAACACAGACGCTCATATGCTCGCCCGACTATCCGTGGGAAACTCGCCGAAGACGACACCACCGAAACGCCGATGTATGTCAACCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAAGTCCTGACTCTGGCCGGTGAGc                                                                                                    <  2:1018874/149‑1 (MQ=255)
                                                aaaCACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGGAGACCAGGGACCAGTAGTGGCTCTTGCCGGGTAGcc                                                                                                   >  1:554618/1‑149 (MQ=255)
                                                  acacAACCTCTCATATGTTCGCCCGACTAACCGAGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTg                                                                                                 <  1:43022/149‑1 (MQ=255)
                                                    ataaaCTCGAATATCTTCGCCCGGCTACCCGTGGAAATTTCGCCCAAGAAGCCACCACCGAGCCGCCGATGTATGCCACTCGCGGTGTTCGCCTTACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGAc                                                                                               <  2:1772450/147‑1 (MQ=255)
                                                         cctctCATATGTTCGCCCGACTAACCGGGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaa                                                                                           <  1:350053/148‑1 (MQ=255)
                                                           tctcATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaaaa                                                                                         <  1:1528067/149‑1 (MQ=255)
                                                           gtctacctGTTCCCCCGATTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaaaa                                                                                         <  2:348419/142‑1 (MQ=255)
                                                                 atGCTCCCCCGACTAACCGTCGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcgc                                                                                   <  1:764673/149‑1 (MQ=255)
                                                                  tgtgCGACCGACTAACCGCGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcgcg                                                                                  <  1:1510391/149‑1 (MQ=255)
                                                                  gcttCGGCCGACTCCCCGTGTGAATATTACCGAAGCATACACCACCGAACCGCCGATGTATGTCAATCGCCGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcgcg                                                                                  <  2:1483737/147‑1 (MQ=255)
                                                                     tCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCgcag                                                                               <  1:637374/149‑1 (MQ=255)
                                                                     ccGCCCGACTAACCGTGGCAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCgcag                                                                               <  1:1056664/148‑1 (MQ=255)
                                                                       gccggccTAAGCGTCGGAAAATTGCCGAAGCAGACCCCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCgcagca                                                                             <  2:945370/143‑1 (MQ=255)
                                                                       gCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGACCATTCTACAGGAAAACCAGAGACAATTACTGACTCTTGCCGGGGAGCCTGACGGCAAGCTCCCGgcagcc                                                                             >  2:1418081/1‑148 (MQ=255)
                                                                        cccGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGcagcag                                                                            <  1:1379421/149‑1 (MQ=255)
                                                                            aCTAACCGTGGGATAATTTCCGTAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCATGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCa                                                                        <  2:36395/149‑1 (MQ=255)
                                                                              tAACCGTGGGAAAATTTCCGAGGAAGACACGACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCAtt                                                                      >  1:941547/1‑149 (MQ=255)
                                                                                     gggAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCa                                                               <  1:1121522/149‑1 (MQ=255)
                                                                                     gggAAAACTTCCGAAGAAGACACCACCGAAACGCCGATGTATGCCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCa                                                               <  1:170879/149‑1 (MQ=255)
                                                                                      ggAAAATTTCCGAAGAAGACCCCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAt                                                              <  1:108801/149‑1 (MQ=255)
                                                                                       gAAGATTTCCGAAGAAGACACCACCGAAACGACGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAtt                                                             <  1:1028792/149‑1 (MQ=255)
                                                                                         aaaTTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCGTTCt                                                           >  1:506510/1‑149 (MQ=255)
                                                                                           aTTTCCGAAGAAGACACCACCGAAACGCCGTTGGATGTCAATCGCAATACGGGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCtta                                                         <  1:1313145/149‑1 (MQ=255)
                                                                                             ttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGTACCATTAAGCAAAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAACACTCCGCGCGGCAGGTCATTAATATCATTCttaat                                                       >  2:545431/1‑149 (MQ=255)
                                                                                                 gaatcagaCAACACCGAACCGACCATGTATGTCAATCGAGGTGGTCGTCGGCAAATTCTGCAGGCAGACCAGGGACAGTTACTGACTCTGGCCGGGGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAAc                                                   <  1:791947/144‑1 (MQ=255)
                                                                                                  aagaacACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACt                                                  <  2:1444015/149‑1 (MQ=255)
                                                                                                    gaagaCACCACCGAAACGCCGATGTTTGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTaa                                                >  1:901880/1‑149 (MQ=255)
                                                                                                    gaagaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGTAAGACCAGGGACAATTACTGACTCGTGCCGGTGAGCCTGACGGAAAACTCGGCGCAGCAGGTCATTAATATCATTCTTAATTAACTaa                                                <  1:1174343/149‑1 (MQ=255)
                                                                                                    gaagaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGGCTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCCGGTCATTAATATCCTTCTTAATTACCTaa                                                >  1:387195/1‑149 (MQ=255)
                                                                                                       gaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCAGAAAGAAGACCAGGGACAATTACTGACTCTTGCCGATGAGCCTGACGGAAAACTCCACGCAGCAGGTCATTAATATCATTCTTAATTAACTActta                                             >  2:927345/1‑149 (MQ=255)
                                                                                                         caccacCGAAACGCCGATGTATGGCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAattatt                                           <  1:490618/149‑1 (MQ=255)
                                                                                                          accaccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAAttattt                                          <  1:923162/149‑1 (MQ=255)
                                                                                                           ccaacGAAACGCCGATATATGTCAATCGCGGTGGTCGACTGACTATTCTGCAGTAAGACCAGGGATAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGTAGTTCATTAATATCATTCTTAATTAACTAAttattta                                         <  2:912356/149‑1 (MQ=255)
                                                                                                             accGGAACCCCGATGTATGGCAATCGCGGTGGTCGCCGGAAGATTCTGCAGGCAGACCAGGGACAATTACTCACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTAtc                                       <  1:849994/149‑1 (MQ=255)
                                                                                                                gAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCAGAAGAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtca                                    >  2:1381947/1‑149 (MQ=255)
                                                                                                                  aaCGACGATGTATCCAAAGCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAACTACTGACTCTAGCCGGTGAGCCTTACGGTATGCTCCGCGCAACAGGTCAGTAATATCATTCTTAATTAACTAATTATTTATCtcatc                                  <  2:1223617/149‑1 (MQ=255)
                                                                                                                   aCGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAACCTCCGCTCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtcatca                                 <  1:1706622/149‑1 (MQ=255)
                                                                                                                     gCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTGTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACt                               >  1:1152338/1‑149 (MQ=255)
                                                                                                                     gCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTAAAAAGAAAAACCAGGGACACTTACTGACTCTTGCCGCTGAGCCTCACGGAAAACCACGAGCAGCAGGTCATTAATATAATTCTTAATCAACTAATTCTTTCTATCATCACt                               >  2:90579/1‑149 (MQ=255)
                                                                                                                     gCCGATGTATGTCAATCGCGGTGGTCGCCGGATCATTCAAAAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCGTTAATATCATTCTTAATTAACTCATTATTTTTCTCAacacg                               >  2:112982/1‑148 (MQ=255)
                                                                                                                     gCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGCAAGAAGACCAGGGACAATGACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACt                               >  2:1334256/1‑149 (MQ=255)
                                                                                                                        gATGTATGACAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGaa                            <  2:1213640/149‑1 (MQ=255)
                                                                                                                         atgtatgtCAATCGCGGTGGTCGCCTGACCATTCAGCAAGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAat                           >  2:433075/1‑149 (MQ=255)
                                                                                                                              tgtCAATCTCGGTGGTCGCCGTATTATTAAAAAAGAAGACAAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTACATATCtt                      >  2:889301/1‑149 (MQ=255)
                                                                                                                              tgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCtt                      >  1:1119994/1‑149 (MQ=255)
                                                                                                                              tgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCtt                      >  1:1571936/1‑149 (MQ=255)
                                                                                                                               gtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTa                     <  1:194512/149‑1 (MQ=255)
                                                                                                                                     cgccGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAg               <  1:585946/149‑1 (MQ=255)
                                                                                                                                         gtgggCGCCTGACTATTCGGCAGGAAGACCCTGGACAATTACTGCCTCTTGCCGGTGAGCCTGACGTAAAACTCCGAGCAGAAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAg           <  2:811192/149‑1 (MQ=255)
                                                                                                                                           ggtCCCCGGACTATTCTGGAGGAAGACCAGGGACAATTACTGACTCTTGCAGGTGAGCCTGACGGAAAATTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGa         <  2:595607/149‑1 (MQ=255)
                                                                                                                                             tCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACt       >  1:636611/1‑149 (MQ=255)
                                                                                                                                               gCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGACTATCTTAATATAGTGAGGACtta     >  1:455966/1‑149 (MQ=255)
                                                                                                                                               gCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACtta     >  1:1071842/1‑149 (MQ=255)
                                                                                                                                                ccTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAGCTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCGCATCACTGAGTATCTTAATATGGTGAGGACtgag    >  1:1509715/1‑146 (MQ=255)
                                                                                                                                                 cTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACttatt   <  2:1564062/149‑1 (MQ=255)
                                                                                                                                                  tGACTATTCTGCAGGAAGACCAGGGACAATTACTTACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACttatta  <  1:1615961/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCACTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACTTATTA  >  NC_000913/2473737‑2474030

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: