Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,633,923 G→A G101E (GGG→GAG)  pinQ → Qin prophage; putative site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,633,9230GA99.0% 216.4 / ‑4.1 96G101E (GGG→GAG) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (43/52);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTT  >  NC_000913/1633789‑1634058
                                                                                                                                      |                                                                                                                                       
gcttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCa                                                                                                                                    <  2:140731/140‑1 (MQ=14)
gcttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCa                                                                                                                                    <  1:1407569/140‑1 (MQ=14)
gcttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCa                                                                                                                                    <  1:1397669/140‑1 (MQ=14)
gcttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCa                                                                                                                                    <  1:1264351/140‑1 (MQ=14)
gcttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCa                                                                                                                                    <  1:1254966/140‑1 (MQ=14)
  ttgctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCAtt                                                                                                                                  >  2:454993/1‑140 (MQ=17)
     ctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgac                                                                                                                               <  2:1366488/140‑1 (MQ=17)
     ctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgac                                                                                                                               <  2:457574/140‑1 (MQ=18)
     ctCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgac                                                                                                                               >  2:967003/1‑140 (MQ=17)
          cTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgacctgac                                                                                                                          <  2:981742/140‑1 (MQ=18)
          cTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgacctgac                                                                                                                          <  2:14119/140‑1 (MQ=18)
          cTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATtgacctgac                                                                                                                          <  2:188740/140‑1 (MQ=18)
                  tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCaa                                                                                                                  <  1:151237/140‑1 (MQ=32)
                  tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCaa                                                                                                                  >  2:989753/1‑140 (MQ=17)
                  tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCaa                                                                                                                  >  2:229079/1‑140 (MQ=17)
                  tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCaa                                                                                                                  >  2:463989/1‑140 (MQ=17)
                   ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAAc                                                                                                                 <  2:809158/140‑1 (MQ=18)
                    gtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACa                                                                                                                >  2:1075801/1‑140 (MQ=17)
                    gtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACa                                                                                                                <  2:670803/140‑1 (MQ=18)
                                 tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatga                                                                                                   <  1:785524/140‑1 (MQ=18)
                                 tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatga                                                                                                   <  1:669854/140‑1 (MQ=18)
                                 tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatga                                                                                                   >  2:683995/1‑140 (MQ=17)
                                 tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatga                                                                                                   >  1:1323932/1‑140 (MQ=17)
                                   tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatgatg                                                                                                 >  1:1315390/1‑140 (MQ=17)
                                   tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAATGATGATGc                                                                                                >  2:457352/1‑140 (MQ=12)
                                      cAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCaa                                                                                              <  2:415266/140‑1 (MQ=18)
                                        aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGt                                                                                            <  1:454993/140‑1 (MQ=18)
                                             ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAAttt                                                                                       >  1:747000/1‑140 (MQ=17)
                                              gATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTc                                                                                      >  2:1377544/1‑140 (MQ=17)
                                                   ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTcagcag                                                                                 >  1:488552/1‑140 (MQ=17)
                                                     ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTc                                                                               >  2:926897/1‑140 (MQ=17)
                                                      tGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                              <  1:787816/140‑1 (MQ=18)
                                                      tGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                              <  2:277034/140‑1 (MQ=17)
                                                       ggTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                             <  2:1315390/140‑1 (MQ=18)
                                                       ggTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                             <  2:524510/140‑1 (MQ=18)
                                                        gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                            <  2:335145/140‑1 (MQ=17)
                                                        gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                            >  2:384159/1‑140 (MQ=17)
                                                            tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAt                                                                        <  1:1028664/140‑1 (MQ=17)
                                                            tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAt                                                                        <  1:1075801/140‑1 (MQ=17)
                                                                 cAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGaa                                                                   <  1:683995/140‑1 (MQ=17)
                                                                  aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAAc                                                                  >  1:726511/1‑140 (MQ=17)
                                                                   aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                            >  2:1076099/1‑129 (MQ=17)
                                                                   aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                            <  1:1076099/129‑1 (MQ=17)
                                                                   aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                 <  2:1409774/140‑1 (MQ=17)
                                                                    tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGGCCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                <  1:764503/140‑1 (MQ=17)
                                                                    tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                <  1:384159/140‑1 (MQ=17)
                                                                         atCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCt                                                           <  1:989753/140‑1 (MQ=17)
                                                                         atCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCt                                                           <  1:918414/140‑1 (MQ=17)
                                                                           cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                         >  2:551874/1‑140 (MQ=17)
                                                                              gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                >  1:1459324/1‑130 (MQ=17)
                                                                              gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                <  2:1459324/130‑1 (MQ=17)
                                                                               aaaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACt                                                     <  2:1323932/140‑1 (MQ=17)
                                                                                aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACtt                                                    >  2:999714/1‑140 (MQ=17)
                                                                                 aactgTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTg                                                   >  1:818047/5‑140 (MQ=16)
                                                                                 aaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTg                                                   >  1:462389/1‑140 (MQ=17)
                                                                                    aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                <  1:602366/140‑1 (MQ=17)
                                                                                     gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                               <  2:726511/140‑1 (MQ=17)
                                                                                      tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                              >  1:1220531/1‑140 (MQ=17)
                                                                                      tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                              >  2:940444/1‑140 (MQ=18)
                                                                                      tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                              >  2:1497154/1‑140 (MQ=17)
                                                                                       ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCa                                             >  2:1357237/1‑140 (MQ=18)
                                                                                        gAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                            >  2:1124381/1‑140 (MQ=18)
                                                                                         aacaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                           <  1:463989/140‑1 (MQ=17)
                                                                                             acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                       <  1:926897/140‑1 (MQ=17)
                                                                                             acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                       <  1:229079/140‑1 (MQ=17)
                                                                                             acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                       <  1:365848/140‑1 (MQ=17)
                                                                                              cTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                      <  1:967003/140‑1 (MQ=17)
                                                                                               tGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTcc                                     <  1:1377544/140‑1 (MQ=17)
                                                                                                gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCg                                    >  2:342025/1‑140 (MQ=17)
                                                                                                  ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                       >  2:196059/1‑135 (MQ=18)
                                                                                                  ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                       <  1:196059/135‑1 (MQ=17)
                                                                                                  ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCggg                                  >  2:1193101/1‑140 (MQ=18)
                                                                                                  ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCggg                                  >  2:61231/1‑140 (MQ=18)
                                                                                                   cGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                 <  1:1193101/140‑1 (MQ=17)
                                                                                                   cGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                 <  1:61231/140‑1 (MQ=17)
                                                                                                      aaCAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAg                              <  1:1124381/140‑1 (MQ=17)
                                                                                                      aaCAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAg                              >  2:1402595/1‑140 (MQ=18)
                                                                                                        cAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTa                            >  2:984301/1‑140 (MQ=18)
                                                                                                        cAGGTATCAGAGTGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTa                            >  2:81868/1‑140 (MQ=18)
                                                                                                                    tGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggc                <  2:462389/140‑1 (MQ=17)
                                                                                                                    tGCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggc                <  1:999714/140‑1 (MQ=17)
                                                                                                                     gCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggcg               <  1:551874/140‑1 (MQ=17)
                                                                                                                     gCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggcg               <  2:1220531/140‑1 (MQ=17)
                                                                                                                     gCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggcg               >  1:351021/1‑140 (MQ=18)
                                                                                                                     gCATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggcg               <  1:1497154/140‑1 (MQ=17)
                                                                                                                      cATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc              >  2:541655/1‑140 (MQ=18)
                                                                                                                      cATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc              <  2:991700/140‑1 (MQ=17)
                                                                                                                      cATTGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc              >  1:991700/1‑140 (MQ=17)
                                                                                                                         tGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg           >  2:289952/1‑140 (MQ=18)
                                                                                                                         tGCTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg           >  1:741469/1‑140 (MQ=18)
                                                                                                                           cTTAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa         >  1:351583/1‑140 (MQ=18)
                                                                                                                            ttAGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTGCAGCAGTCGGTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCCGGATAGTAAGAGCCCGCGGCGCAGGGaa        >  2:1463213/1‑140 (MQ=255)
                                                                                                                              aGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                    >  2:1068912/1‑126 (MQ=18)
                                                                                                                              aGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                    <  1:1068912/126‑1 (MQ=17)
                                                                                                                              aGCATTGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc      >  1:429068/1‑140 (MQ=18)
                                                                                                                                  ttGGAGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGttt  <  2:741469/140‑1 (MQ=17)
                                                                                                                                      |                                                                                                                                       
GCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTT  >  NC_000913/1633789‑1634058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: