Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,177,646 T→G 8.3% N109K (AAT→AAG lolE → lipoprotein‑releasing system transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,177,6460TG 8.3% 273.1 / 13.5 82N109K (AAT→AAGlolElipoprotein‑releasing system transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base G (2/5);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGC  >  NC_000913/1177498‑1177789
                                                                                                                                                    |                                                                                                                                               
tGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAt                                                                                                                                                 <  3:664315/149‑1 (MQ=255)
tGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAt                                                                                                                                                 >  1:515706/1‑149 (MQ=255)
    ggtggtGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTc                                                                                                                                             <  1:775744/149‑1 (MQ=255)
    ggtggtGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTc                                                                                                                                             <  1:25581/149‑1 (MQ=255)
        gtgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgc                                                                                                                                         <  1:859656/149‑1 (MQ=255)
         tgCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgcc                                                                                                                                        <  2:816504/149‑2 (MQ=255)
          gCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCaa                                                                                                                                       >  1:791996/1‑149 (MQ=255)
           ccGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCAAGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCAGGTGGAAAGTGGCGCGAAGATTcgcaccac                                                                                                                                      <  4:112514/149‑6 (MQ=255)
             gCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATcc                                                                                                                                    >  3:730286/1‑149 (MQ=255)
             gCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGACTCgcgccactt                                                                                                                                    <  3:871939/149‑6 (MQ=255)
               aTGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAg                                                                                                                                  <  2:833223/149‑1 (MQ=255)
                tGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAgg                                                                                                                                 <  3:672209/149‑1 (MQ=255)
                tGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAgg                                                                                                                                 >  1:292026/1‑149 (MQ=255)
                   tGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGa                                                                                                                              >  1:96358/1‑149 (MQ=255)
                       aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAggg                                                                                                                          >  1:408095/1‑149 (MQ=255)
                        tCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGc                                                                                                                         >  3:362493/1‑149 (MQ=255)
                         cGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCg                                                                                                                        <  1:102184/149‑1 (MQ=255)
                              cGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAc                                                                                                                   <  1:710958/149‑1 (MQ=255)
                              cGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAc                                                                                                                   <  2:475919/149‑1 (MQ=255)
                              cGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAc                                                                                                                   <  1:262783/149‑1 (MQ=255)
                              cGGTGGAGTAGCCGTGGACTGCCTGGGAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTGACCAGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAc                                                                                                                   <  2:868285/149‑1 (MQ=255)
                                    aTCCGCGGTGGGCTAACTGGAAGGAAGGGCTGGATCACGTGCAAAAGGTGCCAGGTATTGCCGCTGCCGCGCCGTATCTCAATTTCACCCGGTTGGTGAAAAGTGGCGCGAATCGCCGCGCAATCCAGGTGATGGGCGTTAAACCGaaa                                                                                                             >  2:312509/1‑149 (MQ=255)
                                    aTCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCaa                                                                                                             >  4:23690/1‑149 (MQ=255)
                                              gACTAACTGGCAGGAAGCACTGGATCACGTCCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGAGAAGGGCGTTAACCCGCAACAGGAACAGc                                                                                                   <  3:717530/149‑1 (MQ=255)
                                                    cTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGa                                                                                             >  3:321116/1‑149 (MQ=255)
                                                    cTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGa                                                                                             <  4:362493/149‑1 (MQ=255)
                                                     tGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAg                                                                                            >  2:136144/1‑149 (MQ=255)
                                                         aGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCa                                                                                        >  1:594471/1‑149 (MQ=255)
                                                         aGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCa                                                                                        >  2:675472/1‑149 (MQ=255)
                                                         aGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCa                                                                                        >  2:287358/1‑149 (MQ=255)
                                                          ggAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCAt                                                                                       >  4:184101/1‑149 (MQ=255)
                                                              gCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAcc                                                                                   <  3:743800/149‑1 (MQ=255)
                                                               cACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                          >  4:579100/1‑125 (MQ=255)
                                                               cACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAg                                                                                                          <  3:579100/125‑1 (MQ=255)
                                                                aCTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCt                                                                                 >  1:32462/1‑149 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  <  3:408624/115‑1 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  >  4:408624/1‑115 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                <  4:730286/149‑1 (MQ=255)
                                                                 cTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTCCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                >  3:823842/1‑149 (MQ=255)
                                                                  tGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCg                                                                               >  4:144731/1‑149 (MQ=255)
                                                                   ggATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGt                                                                              <  4:179863/149‑1 (MQ=255)
                                                                        aCGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGaa                                                                                                       >  1:378355/1‑119 (MQ=255)
                                                                        aCGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGaa                                                                                                       <  2:378355/119‑1 (MQ=255)
                                                                          gTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCa                                                                       >  2:441030/1‑149 (MQ=255)
                                                                           tGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCNg                                                                      <  3:23690/149‑1 (MQ=255)
                                                                             cAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAggg                                                                    >  2:80263/1‑149 (MQ=255)
                                                                             cAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgcca                                                                                                                                       <  1:634851/82‑3 (MQ=255)
                                                                             cAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAAGATTCgcgcca                                                                                                                                       >  2:634851/1‑80 (MQ=255)
                                                                              aaaaaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGc                                                                   >  2:846368/1‑149 (MQ=255)
                                                                                 aaGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGAt                                                                >  1:712173/1‑149 (MQ=255)
                                                                                   gTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGc                                                              <  1:820839/149‑1 (MQ=255)
                                                                                      ccGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTg                                                           >  3:595233/1‑149 (MQ=255)
                                                                                      ccGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTg                                                           <  4:595233/149‑1 (MQ=255)
                                                                                       cGGGGATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGGGGAAAGTGGGGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAACCGCAACAGGAACAGCGTCTGAGCGAATTACCCTCGTTCGTCAAGGGCGATGCGTgg                                                          <  4:441712/149‑1 (MQ=255)
                                                                                            aTTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCaa                                                     <  2:292026/149‑1 (MQ=255)
                                                                                             ttGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  <  2:96305/87‑1 (MQ=255)
                                                                                             ttGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  >  1:816189/1‑87 (MQ=255)
                                                                                             ttGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  <  2:816189/87‑1 (MQ=255)
                                                                                             ttGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAAcc                                                                                                                  >  1:96305/1‑87 (MQ=255)
                                                                                                 cgctgccgcGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTa                                                >  2:211970/1‑149 (MQ=255)
                                                                                                       cgcgCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCgg                                          <  4:724623/149‑1 (MQ=255)
                                                                                                        gcgcCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCggg                                         >  1:848563/1‑149 (MQ=255)
                                                                                                           ccGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                <  1:348152/107‑1 (MQ=255)
                                                                                                           ccGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTc                                                                                >  2:348152/1‑107 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGGAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                     <  2:851396/149‑1 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                     >  1:758734/1‑149 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                     <  2:426964/149‑1 (MQ=255)
                                                                                                            cGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTAACCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                     <  3:602352/149‑1 (MQ=255)
                                                                                                             gTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAAc                                    <  3:154762/149‑1 (MQ=255)
                                                                                                                  tCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaa                               >  4:417782/1‑149 (MQ=255)
                                                                                                                   cAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaaa                              <  1:136144/149‑1 (MQ=255)
                                                                                                                   cAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaaa                              <  2:338284/149‑1 (MQ=255)
                                                                                                                     aTTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAtt                            >  2:823744/1‑149 (MQ=255)
                                                                                                                        tCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatc                         <  4:188687/149‑1 (MQ=255)
                                                                                                                           ccGGGCGGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcatc                      <  4:8337/149‑1 (MQ=255)
                                                                                                                                 tggtggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCaaa                >  4:841745/1‑149 (MQ=255)
                                                                                                                                  ggtggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAg               <  1:243915/149‑1 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCg            >  4:452944/1‑149 (MQ=255)
                                                                                                                                     ggAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCg            >  1:662060/1‑149 (MQ=255)
                                                                                                                                         gattgcgcGAAGATTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGc        >  2:686603/5‑149 (MQ=255)
                                                                                                                                         aGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGc        <  1:80263/149‑1 (MQ=255)
                                                                                                                                         aGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGc        <  1:268/149‑1 (MQ=255)
                                                                                                                                              cgcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATg   <  3:585594/149‑1 (MQ=255)
                                                                                                                                               gcgAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATgc  >  3:278658/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGGTATTGCCGCTGCCGCGCCGTATATCAATTTCACCGGGCTGGTGGAAAGTGGCGCGAATCTTCGCGCAATCCAGGTGAAGGGCGTTAACCCGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGC  >  NC_000913/1177498‑1177789

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