Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,358 T→A 5.8% L17Q (CTG→CAG)  ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,3580TA 5.8% 785.8 / 26.3 220L17Q (CTG→CAG) ybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base T (92/115);  new base A (7/6);  total (99/121)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 223 total aligned reads displayed.

CCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGC  >  NC_000913/774108‑774631
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
ccTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgc                                                                                                                                                                                                                                                           >  4:240848/1‑275 (MQ=255)
           tGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAg                                                                                                                                                                                                                                                <  4:18152/275‑1 (MQ=255)
              aaaaTGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTCACATCGTGATGGCGCt                                                                                                                                                                                                                                                        >  2:231877/1‑264 (MQ=255)
                                    aagaagATAGTGACCGTAACCCCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGAGGCATTTCGCAGGGATTTTGGGAACGCTGCTTGCCTGTTCGTTTGGGGTAATCCCCGCGCTGGCGCTTGAACCCGTCGAATCAGGCAAAGCCGGTCAAGAAAACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtc                                                                                                                                                                                                                       <  4:270511/275‑1 (MQ=255)
                                                 aCGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgcg                                                                                                                                                                                                          <  2:240848/275‑1 (MQ=255)
                                                    tAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACGTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGGGATGGCGCTTCTGTTAGCCGGTTGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAg                                                                                                                                                                                                       >  1:12965/1‑275 (MQ=255)
                                                     aaCACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCt                                                                                                                                                                                                                                                        <  2:69935/225‑1 (MQ=255)
                                                                              ggAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGGCGCGAAGACCAGTGCACTGGAAATCTGGCTTgg                                                                                                                                                                             >  1:74881/1‑275 (MQ=255)
                                                                                                      aTGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCt                                                                                                                                                                                                                                                                                    >  6:1118895/1‑149 (MQ=255)
                                                                                                           ttCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCggg                                                                                                                                                                                                                                                                               >  6:1006489/1‑149 (MQ=255)
                                                                                                              tGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGgcgc                                                                                                                                                                                                                                                                            <  7:249569/149‑1 (MQ=255)
                                                                                                              tGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGgcgc                                                                                                                                                                                                                                                                            <  7:3322/149‑1 (MQ=255)
                                                                                                               gggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGTTTACCGCGACTTTGTATGCGGTAATGGACAAGGGCCCCCTGCGGGCGCt                                                                                                                                                                                                                                                                           >  5:465222/1‑149 (MQ=255)
                                                                                                                ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGCGTACGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTCTTGCTGATGTGGGCCGTCTGTGCTGgtgt                                                                                                                                            >  2:306795/1‑274 (MQ=255)
                                                                                                                ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCtt                                                                                                                                                                                                                                                                          >  7:302806/1‑149 (MQ=255)
                                                                                                                   aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTcc                                                                                                                                                                                                                                                                      >  8:1360764/1‑149 (MQ=255)
                                                                                                                   aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTcc                                                                                                                                                                                                                                                                      <  8:692016/149‑1 (MQ=255)
                                                                                                                     gCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCtt                                                                                                                                                                                                                                                                    >  8:535498/1‑149 (MQ=255)
                                                                                                                      cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTc                                                                                                                                                                                                                                                                   >  8:224336/1‑149 (MQ=255)
                                                                                                                      cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTc                                                                                                                                                                                                                                                                   <  6:515153/149‑1 (MQ=255)
                                                                                                                      cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTc                                                                                                                                                                                                                                                                   <  6:30454/149‑1 (MQ=255)
                                                                                                                       ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCg                                                                                                                                                                                                                                                                  >  8:825692/1‑149 (MQ=255)
                                                                                                                        tcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACCCGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                                                                                                                                                 <  8:961264/149‑1 (MQ=255)
                                                                                                                        tcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                                                                                                                                               <  4:358396/167‑1 (MQ=255)
                                                                                                                          ttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGt                                                                                                                                                                                                                                                   >  4:139924/1‑161 (MQ=255)
                                                                                                                              cTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgc                                                                                                                                                                                                                                                           >  8:172181/1‑149 (MQ=255)
                                                                                                                              cTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                                                                                                                                               <  2:139924/161‑1 (MQ=255)
                                                                                                                                gTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgcgc                                                                                                                                                                                                                                                         >  6:251376/1‑149 (MQ=255)
                                                                                                                                 ttCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATGAAGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTAGGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCGTTCGTGATGGCGCt                                                                                                                                                                                                                                                        >  8:898104/1‑149 (MQ=255)
                                                                                                                                            tAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAg                                                                                                                                                                                                                                             <  6:1251738/149‑1 (MQ=255)
                                                                                                                                            tAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATCAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAg                                                                                                                                                                                                                                             >  5:330650/1‑149 (MQ=255)
                                                                                                                                               tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                                                                                                                                                          <  5:13554/149‑1 (MQ=255)
                                                                                                                                               tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                                                                                                                                                          <  6:305443/149‑1 (MQ=255)
                                                                                                                                               tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                                                                                                                                                          <  8:150423/149‑1 (MQ=255)
                                                                                                                                               tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                                                                                                                                                          <  6:596445/149‑1 (MQ=255)
                                                                                                                                                        tGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGtttt                                                                                                                                                                                                                                 <  5:1238477/149‑1 (MQ=255)
                                                                                                                                                        tGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGtttt                                                                                                                                                                                                                                 <  8:302806/149‑1 (MQ=255)
                                                                                                                                                         ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgc                                                                                                                                                                                                                                                           <  7:466172/122‑1 (MQ=255)
                                                                                                                                                         ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgc                                                                                                                                                                                                                                                           >  8:466172/1‑122 (MQ=255)
                                                                                                                                                         ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                                                                                                                                                <  5:291561/149‑1 (MQ=255)
                                                                                                                                                         ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                                                                                                                                                <  7:1044550/149‑1 (MQ=255)
                                                                                                                                                         ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                                                                                                                                                <  7:282907/149‑1 (MQ=255)
                                                                                                                                                            gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAACAGGATGTATGTTTTGGGa                                                                                                                                                                                                                             <  8:1238052/149‑1 (MQ=255)
                                                                                                                                                             cTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                                                                                                                                                            <  7:521977/149‑1 (MQ=255)
                                                                                                                                                               tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                                                                                                                                                          <  6:883293/149‑1 (MQ=255)
                                                                                                                                                                 aaCACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAt                                                                                                                                                                                                                        <  5:1428823/149‑1 (MQ=255)
                                                                                                                                                                  acacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtc                                                                                                                                                                                                                       >  7:1224474/1‑149 (MQ=255)
                                                                                                                                                                  acacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtc                                                                                                                                                                                                                       >  5:1287822/1‑149 (MQ=255)
                                                                                                                                                                   cacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtct                                                                                                                                                                                                                      >  6:1214254/1‑149 (MQ=255)
                                                                                                                                                                    acGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCt                                                                                                                                                                                                                                                                     <  6:947163/101‑1 (MQ=255)
                                                                                                                                                                    acGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtctc                                                                                                                                                                                                                     >  8:490657/1‑149 (MQ=255)
                                                                                                                                                                      gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGt                                                                                                                                                                                                                   >  5:108227/1‑149 (MQ=255)
                                                                                                                                                                         gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtttt                                                                                                                                                                                                                <  8:1279094/149‑1 (MQ=255)
                                                                                                                                                                         gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtttt                                                                                                                                                                                                                <  5:1038535/149‑1 (MQ=255)
                                                                                                                                                                               ggCGAAGCCGGTCAAGAAGACATCTGAGGCGTAAGATTATCGCGACATTGCACGCGGTCATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCTTGATGGCGCTTCTGTTAGCAGGATGTCTGTTTTGGGACCTATCTCGGTTTGCcggg                                                                                                                                                                                                          >  8:562186/1‑147 (MQ=255)
                                                                                                                                                                               ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgcg                                                                                                                                                                                                          >  5:688834/1‑149 (MQ=255)
                                                                                                                                                                                 cAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                                                                                                                                                        <  5:261820/149‑1 (MQ=255)
                                                                                                                                                                                 cAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                                                                                                                                                        <  5:673415/149‑1 (MQ=255)
                                                                                                                                                                                  aaaGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAg                                                                                                                                                                                                       <  7:57937/149‑1 (MQ=255)
                                                                                                                                                                                  aaaGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAg                                                                                                                                                                                                       >  7:377472/1‑149 (MQ=255)
                                                                                                                                                                                    aGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGACATCTGGCATGGGTTATTGCTGATGTGGGCCGACGGGGTTGGTGTGATTCCCGGCGTGTGCTTTCGTCCGCGGCAGGTCTTTTGGCAAGGGTTTTTTTGCCCAtttcttgcact                                                                       >  4:17268/1‑272 (MQ=255)
                                                                                                                                                                                      ccGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGat                                                                      <  2:144592/274‑1 (MQ=255)
                                                                                                                                                                                       cGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAg                                                                                                                                                                                                  >  5:236102/1‑149 (MQ=255)
                                                                                                                                                                                        ggTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGt                                                                                                                                                                                                 <  8:1127771/149‑1 (MQ=255)
                                                                                                                                                                                          tCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGa                                                                                                                                                                                               <  6:108227/149‑1 (MQ=255)
                                                                                                                                                                                             agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                                                                                                                                                            >  5:536131/1‑149 (MQ=255)
                                                                                                                                                                                             agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                                                                                                                                                            <  5:1403838/149‑1 (MQ=255)
                                                                                                                                                                                               aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                                                                                                                                                          >  7:15076/1‑149 (MQ=255)
                                                                                                                                                                                                 gaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaa                                                                                                                                                                                        <  5:1006489/149‑1 (MQ=255)
                                                                                                                                                                                                      cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTg                                                                                                                                                                                   <  7:1299766/149‑1 (MQ=255)
                                                                                                                                                                                                           gaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATg                                                                                                                                                                              <  5:367615/149‑1 (MQ=255)
                                                                                                                                                                                                             gTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                                                                                                                                                            <  6:536131/149‑1 (MQ=255)
                                                                                                                                                                                                             gTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                                                                                                                                                            <  6:403152/149‑1 (MQ=255)
                                                                                                                                                                                                              tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                                                                                                                                                >  6:1405279/1‑96 (MQ=255)
                                                                                                                                                                                                              tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                                                                                                                                                           >  5:447297/1‑149 (MQ=255)
                                                                                                                                                                                                              tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                                                                                                                                                           >  6:1374394/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTg                                                                                                                                                                      <  7:222125/149‑1 (MQ=255)
                                                                                                                                                                                                                      tCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTTTCATTGTCGGGCTGATTttcttcttc                                     >  2:113727/1‑275 (MQ=255)
                                                                                                                                                                                                                       cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGa                                                                                                                                                                  <  5:488180/149‑1 (MQ=255)
                                                                                                                                                                                                                        gcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                                                                                                                                                          <  7:968134/124‑1 (MQ=255)
                                                                                                                                                                                                                        gcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAt                                                                                                                                                                 <  8:25384/149‑1 (MQ=255)
                                                                                                                                                                                                                            cTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                                                                                                                                                             <  5:800775/149‑1 (MQ=255)
                                                                                                                                                                                                                                          aTGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgcg                                                                                                                                                                                                          <  4:101638/90‑1 (MQ=255)
                                                                                                                                                                                                                                            ggACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtg                                                                                                                                             >  5:216712/1‑149 (MQ=255)
                                                                                                                                                                                                                                             gACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt                                                                                                                                            <  6:586142/149‑1 (MQ=255)
                                                                                                                                                                                                                                              gaaaGCGCCCGCAGGGGGCGCTTGTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt                                                                                                                                            >  6:1261777/3‑149 (MQ=255)
                                                                                                                                                                                                                                              gaaaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                                                                                                                                                            >  8:1340901/3‑116 (MQ=255)
                                                                                                                                                                                                                                              gaaaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                                                                                                                                                            <  7:1340901/114‑1 (MQ=255)
                                                                                                                                                                                                                                                aaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                                                                                                                                                          >  6:113321/1‑68 (MQ=38)
                                                                                                                                                                                                                                                aaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                                                                                                                                                          <  5:113321/68‑1 (MQ=38)
                                                                                                                                                                                                                                                aaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtgt                                                                                                                                                               <  6:1037321/127‑1 (MQ=255)
                                                                                                                                                                                                                                                aaGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTGTGTCAGCAGGATGTATGTTTTGGGACCCATCTCGTTGTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAGTGCTGATGTTGGCCGTATGTGCGGGTGTGAt                                                                                                                                         >  7:901622/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aaGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgt                                                                                                                                                    <  5:385587/138‑1 (MQ=255)
                                                                                                                                                                                                                                                  gcgcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTc                                                                                                                                       <  7:1358832/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  gcgcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTGTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTc                                                                                                                                       >  6:352560/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCAc                                                                                                                                     <  7:820286/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ccccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                                                                                                                                                >  6:1179617/1‑56 (MQ=255)
                                                                                                                                                                                                                                                       cccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGc                                                                                                                                  >  6:791748/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         ttGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGAACCCATTTCGTTTTGCCGCGAAGACCTGTGAAATGGAAATGTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGc  <  2:201058/274‑1 (MQ=255)
                                                                                                                                                                                                                                                         gtGCGGGGGCTTTCCTTGGGGAGGGCGCTTGTGTTAGCAGGATGTATGTTTTGGGACCCATTTCGTTTGGCCGCGAAGACCAGTGAACTGGAAAGCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGGGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGc  <  4:231825/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
CCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCGCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGC  >  NC_000913/774108‑774631

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