Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,757,688 T→A 6.4% intergenic (+31/+33) lpp → / ← ldtE murein lipoprotein/murein L,D‑transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,6880TA 6.4% 314.0 / 8.2 94intergenic (+31/+33)lpp/ldtEmurein lipoprotein/murein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base A (4/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGATAC  >  NC_000913/1757552‑1757835
                                                                                                                                        |                                                                                                                                                   
gAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttt                                                                                                                                         <  1:533693/149‑1 (MQ=255)
gAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttt                                                                                                                                         <  2:477019/149‑1 (MQ=255)
    aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtttttca                                                                                                                                     <  3:592008/149‑3 (MQ=255)
    aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtt                                                                                                                                          <  4:22155/144‑1 (MQ=255)
    aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtt                                                                                                                                          >  3:22155/1‑144 (MQ=255)
     acgacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttttt                                                                                                                                    >  4:468139/1‑149 (MQ=255)
     acgacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttttt                                                                                                                                    >  2:663681/1‑149 (MQ=255)
       gacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGt                                                                                                                                  <  1:1043366/149‑1 (MQ=255)
       gacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGt                                                                                                                                  <  4:912880/149‑1 (MQ=255)
               cGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGtt                                                                                                                          <  3:34238/149‑1 (MQ=255)
                     gCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGc                                                                                                                    <  4:637083/149‑1 (MQ=255)
                             aCGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGt                                                                                                            >  3:344102/1‑149 (MQ=255)
                               gTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCa                                                                                                          >  4:957554/1‑149 (MQ=255)
                                ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCAc                                                                                                         <  1:424216/149‑1 (MQ=255)
                                        gctAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACa                                                                                                 >  3:674038/1‑149 (MQ=255)
                                              gATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTcc                                                                                           >  3:261307/1‑149 (MQ=255)
                                              gATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTcc                                                                                           >  2:602793/1‑149 (MQ=255)
                                              gATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTcc                                                                                           <  3:938232/149‑1 (MQ=255)
                                                tGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCt                                                                                         <  2:841562/149‑1 (MQ=255)
                                                tGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCt                                                                                         <  3:957554/149‑1 (MQ=255)
                                                 gACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCtt                                                                                        <  4:86780/149‑1 (MQ=255)
                                                 gACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCtt                                                                                        >  2:802863/1‑149 (MQ=255)
                                                 gACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCtt                                                                                        <  3:727991/149‑1 (MQ=255)
                                                    gNAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGct                                                                                     >  3:261461/1‑149 (MQ=255)
                                                    gCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGct                                                                                     >  1:417686/1‑149 (MQ=255)
                                                     cAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGctc                                                                                    <  3:426447/149‑1 (MQ=255)
                                                       gCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTg                                                                                  >  3:274173/1‑149 (MQ=255)
                                                        cTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTgg                                                                                 <  1:663681/149‑1 (MQ=255)
                                                         tCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGt                                                                                >  3:182839/1‑149 (MQ=255)
                                                         tCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGt                                                                                >  2:523573/1‑149 (MQ=255)
                                                           gTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTc                                                                              >  1:137101/1‑149 (MQ=255)
                                                             gCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCAc                                                                            <  1:510182/149‑1 (MQ=255)
                                                               tAACCAGCGTCTGGACAAAATGGCTACTAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCAt                                                                         >  2:480016/1‑149 (MQ=255)
                                                                  ccAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTc                                                                       >  1:894486/1‑149 (MQ=255)
                                                                   cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCt                                                                      >  3:7710/1‑149 (MQ=255)
                                                                   cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCt                                                                      <  4:646354/149‑1 (MQ=255)
                                                                    aGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTg                                                                     >  4:199861/1‑149 (MQ=255)
                                                                     gCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTgc                                                                    >  3:244625/1‑149 (MQ=255)
                                                                     gCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTgc                                                                    <  4:4450/149‑1 (MQ=255)
                                                                           ggACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAc                                                              <  2:913954/149‑1 (MQ=255)
                                                                           ggACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAc                                                              >  2:4077/1‑149 (MQ=255)
                                                                               aacaTGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTa                                                          <  3:1029637/149‑1 (MQ=255)
                                                                               aacaTGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTa                                                          <  3:199172/149‑1 (MQ=255)
                                                                                   tGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACAGTTCCTGCGCCACGCGTAACATATTCCCTTGCTCTGGTTCCCCATTCTGCGCTGGCTCTAcgaa                                                      >  1:942819/1‑146 (MQ=255)
                                                                                   tGGCTACTAAATACCGCAAGTAATAGTACCAGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCATAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGa                                                      <  1:690843/149‑1 (MQ=255)
                                                                                     gCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAg                                                    <  2:800940/149‑1 (MQ=255)
                                                                                     gCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAg                                                    <  3:87975/149‑1 (MQ=255)
                                                                                       tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAGAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTACCATATTCCCTTGCTCTGGTTCACCATTCTGCTCTGACTCTACTGAAGgc                                                  >  4:884955/1‑149 (MQ=255)
                                                                                       tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                                  >  1:53951/1‑149 (MQ=255)
                                                                                         ctaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACatat                                                                                              >  2:671734/1‑103 (MQ=255)
                                                                                         ctaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACatat                                                                                              <  1:671734/103‑1 (MQ=255)
                                                                                         ctaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgcgc                                                >  3:363650/1‑149 (MQ=255)
                                                                                         ctaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgcgc                                                <  1:478087/149‑1 (MQ=255)
                                                                                          taAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCa                                               >  1:678465/1‑149 (MQ=255)
                                                                                          taAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCa                                               >  2:725216/1‑149 (MQ=255)
                                                                                           aaaTACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCAt                                              <  4:756537/149‑1 (MQ=255)
                                                                                               aCCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTgct                                          <  1:350856/149‑1 (MQ=255)
                                                                                                  gCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACAGAAGGCGCATTgctggc                                       >  1:262915/1‑149 (MQ=255)
                                                                                                       tAATAGTACTTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTGCTGAAGGCGCATTGCTGGCTGCgg                                  <  3:884955/149‑1 (MQ=255)
                                                                                                       tAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCgg                                  <  1:631309/149‑1 (MQ=255)
                                                                                                        aaTAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCggg                                 <  4:658537/149‑1 (MQ=255)
                                                                                                         aTAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGTTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGa                                >  1:775538/1‑149 (MQ=255)
                                                                                                         aTAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGa                                <  3:199861/149‑1 (MQ=255)
                                                                                                           aGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGt                              <  1:75949/149‑1 (MQ=255)
                                                                                                             tACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTg                            >  1:668595/1‑149 (MQ=255)
                                                                                                               ccTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCt                          <  3:468139/149‑1 (MQ=255)
                                                                                                                cTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCAtt                                             <  1:267562/129‑1 (MQ=255)
                                                                                                                cTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCAtt                                             >  2:267562/1‑129 (MQ=255)
                                                                                                                  gtgaagtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCa                       >  1:318662/1‑149 (MQ=255)
                                                                                                                     aagtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTg                    >  4:932505/1‑149 (MQ=255)
                                                                                                                       gtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCt                  <  2:62541/149‑1 (MQ=255)
                                                                                                                        tgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTc                 >  1:1025292/1‑149 (MQ=255)
                                                                                                                        tgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTc                 >  2:422377/1‑149 (MQ=255)
                                                                                                                         gaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTc                         <  2:85220/140‑1 (MQ=255)
                                                                                                                         gaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTc                         >  1:85220/1‑140 (MQ=255)
                                                                                                                         aaaaaaTGGCGCACAATGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCa                >  3:625037/2‑149 (MQ=255)
                                                                                                                          aaaaaTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGTTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCAc               <  2:775538/149‑1 (MQ=255)
                                                                                                                          aaaaaTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCAc               >  4:136662/1‑149 (MQ=255)
                                                                                                                            aaaTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGc                           >  4:13420/1‑135 (MQ=255)
                                                                                                                            aaaTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGc                           <  3:13420/135‑1 (MQ=255)
                                                                                                                             aaTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGa            >  4:62714/1‑149 (MQ=255)
                                                                                                                             aaTGGCGCACAATGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGa            >  1:824668/1‑149 (MQ=255)
                                                                                                                              aTGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGaa           >  2:499159/1‑149 (MQ=255)
                                                                                                                               tGGCGCACAATGTGCGCCATTTTTTTTGACTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGCCTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGaaa          >  4:645315/1‑149 (MQ=255)
                                                                                                                                ggCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAAc         <  3:259401/149‑1 (MQ=255)
                                                                                                                                  cgcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCg       <  4:182839/149‑1 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGActc                                                            >  2:97710/1‑95 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGActc                                                            <  1:97710/95‑1 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg      >  2:416707/1‑149 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg      >  3:290328/1‑149 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg      >  3:121496/1‑149 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg      >  3:73051/1‑149 (MQ=255)
                                                                                                                                   gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg      >  3:49715/1‑149 (MQ=255)
                                                                                                                                       aTTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGATAc  <  2:481192/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
GAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGATAC  >  NC_000913/1757552‑1757835

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: