Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,136 G→C 10.0% D22H (GAC→CAC)  ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,1360GC10.0% 271.8 / 19.4 80D22H (GAC→CAC) ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base C (4/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGG  >  NC_000913/1835997‑1836228
                                                                                                                                           |                                                                                            
ccGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCaa                                                                                      <  1:45418/148‑1 (MQ=255)
          tCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaaa                                                                                  <  4:1031421/142‑1 (MQ=255)
          tCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaaa                                                                                  >  3:1031421/1‑142 (MQ=255)
             gCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATg                                                                         <  3:775948/148‑1 (MQ=255)
                       ttCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAgg                                                               >  3:392887/1‑148 (MQ=255)
                       ttCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     >  3:1054764/1‑124 (MQ=255)
                       ttCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     <  4:1054764/126‑3 (MQ=255)
                        tCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATg                                                                         >  1:1023/1‑137 (MQ=255)
                        tCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGGGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATg                                                                         <  2:1023/137‑1 (MQ=255)
                          ccGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatgt                                                                                    >  4:824893/1‑121 (MQ=255)
                          ccGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatgt                                                                                    <  3:824893/124‑4 (MQ=255)
                            gCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaa                                                                              <  2:848213/128‑1 (MQ=255)
                            gCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaa                                                                              >  1:848213/1‑128 (MQ=255)
                                    gCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaa                                                                                   <  3:234985/115‑1 (MQ=255)
                                    gCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaa                                                                                   >  4:234985/1‑115 (MQ=255)
                                     cTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGt                                                 <  4:392887/148‑1 (MQ=255)
                                            gTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGt                                                      <  1:99874/136‑1 (MQ=255)
                                            gTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGt                                                      >  2:99874/1‑136 (MQ=255)
                                                  cTGAAAAACGTTCTTAAGGAGAGCTGATGCGCCATTGTGGGTGGTTACTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGAc                                                                                            <  2:1040450/92‑1 (MQ=255)
                                                  cTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGAc                                                                                            >  1:1040450/1‑92 (MQ=255)
                                                   tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaa                                                                                   <  4:26724/100‑1 (MQ=255)
                                                   tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGaa                                                                                   >  3:26724/1‑100 (MQ=255)
                                                   tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             <  2:1122685/138‑1 (MQ=255)
                                                   tGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             >  1:1122685/1‑138 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGAc                                                                                            >  3:850931/1‑90 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGAc                                                                                            <  4:850931/90‑1 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCa                                                                                       >  1:774635/1‑95 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCa                                                                                       <  2:774635/95‑1 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGt                                                 >  4:534767/1‑133 (MQ=255)
                                                    gAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGt                                                 <  3:534767/133‑1 (MQ=255)
                                                     aaaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcgg                                 <  4:821468/148‑1 (MQ=255)
                                                     aaaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcgg                                 >  3:821468/1‑148 (MQ=255)
                                                      aaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             <  3:316722/135‑1 (MQ=255)
                                                      aaaaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             >  4:316722/1‑135 (MQ=255)
                                                        aaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGTTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAgg                                                                      <  4:1135783/108‑1 (MQ=255)
                                                        aaCGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGTTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAgg                                                                      >  3:1135783/1‑108 (MQ=255)
                                                          cGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     >  2:981049/1‑89 (MQ=255)
                                                          cGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     <  1:981049/91‑3 (MQ=255)
                                                             tCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAggg                                                              >  1:196215/1‑111 (MQ=255)
                                                             tCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAggg                                                              <  2:196215/111‑1 (MQ=255)
                                                               ttAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCa                       >  1:812037/1‑148 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaaaaa                                                                           >  1:1015205/1‑95 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTaaaaa                                                                           <  2:1015205/95‑1 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAggg                                                              <  4:716417/108‑1 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAggg                                                              >  3:716417/1‑108 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCaa                      <  2:812037/148‑1 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     <  1:144051/85‑3 (MQ=255)
                                                                tAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTCACTGGCatg                                                                                     >  2:144051/1‑83 (MQ=255)
                                                                  aGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCaaa                                                                 <  4:597901/103‑1 (MQ=255)
                                                                  aGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCaaa                                                                 >  3:597901/1‑103 (MQ=255)
                                                                    gagagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             >  4:923530/1‑121 (MQ=255)
                                                                    gagagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             <  3:923530/121‑1 (MQ=255)
                                                                     agagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             <  1:1086227/120‑1 (MQ=255)
                                                                     agagaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTaa                                             >  2:1086227/1‑120 (MQ=255)
                                                                       agaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCaaa                                                                 >  2:851007/1‑98 (MQ=255)
                                                                       agaTGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCaaa                                                                 <  1:851007/98‑1 (MQ=255)
                                                                                 ccATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcg                                  >  4:176342/1‑119 (MQ=255)
                                                                                 ccATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGggcg                                  <  3:176342/119‑1 (MQ=255)
                                                                                             gtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                  <  3:1181/123‑1 (MQ=255)
                                                                                             gtTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                  >  4:1181/1‑123 (MQ=255)
                                                                                              ttGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTgg                                                  >  3:803740/1‑90 (MQ=255)
                                                                                              ttGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTgg                                                  <  4:803740/90‑1 (MQ=255)
                                                                                                   gggCTTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATa                            >  4:322465/1‑107 (MQ=255)
                                                                                                   gggATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATa                            <  3:322465/107‑1 (MQ=255)
                                                                                                    ggATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTa                                              <  3:629472/88‑1 (MQ=255)
                                                                                                    ggATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTa                                              >  4:629472/1‑88 (MQ=255)
                                                                                                       ttGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGAc     <  2:984846/126‑1 (MQ=255)
                                                                                                       ttGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGAc     >  1:984846/1‑126 (MQ=255)
                                                                                                              cGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAAc                 >  1:515643/1‑107 (MQ=255)
                                                                                                              cGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAAc                 <  2:515643/107‑1 (MQ=255)
                                                                                                               gCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATAcc                          <  4:119161/97‑1 (MQ=255)
                                                                                                               gCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATAcc                          >  3:119161/1‑97 (MQ=255)
                                                                                                                         ctGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGa      >  1:422132/1‑107 (MQ=255)
                                                                                                                         ctGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGa      <  2:422132/107‑1 (MQ=255)
                                                                                                                           ggCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                  >  1:394905/1‑93 (MQ=255)
                                                                                                                           ggCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTaa                  <  2:394905/93‑1 (MQ=255)
                                                                                                                               acacATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCa                       <  1:782297/84‑1 (MQ=255)
                                                                                                                               acacATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCa                       >  2:782297/1‑84 (MQ=255)
                                                                                                                               acacATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTgg  >  4:1135107/1‑105 (MQ=255)
                                                                                                                               acacATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTgg  <  3:1135107/105‑1 (MQ=255)
                                                                                                                                           |                                                                                            
CCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGATGATGCGCCATTGTGGGTGGTTGCTGGGATTGTTATCGCTGTTTTCTCTGGCAACACATGCCAGTGACTGGCAAGAAATTAAAAATGAGGCCAAAGGGCAAACCGTCTGGTTTAACGCCTGGGGCGGCGATACCGCAATTAACCGCTATCTCGACTGG  >  NC_000913/1835997‑1836228

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